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ABSTRACT 


A  computer  program  has  been  developed  to  perform  a  first  order 
elastic  analysis  of  a  three-dimensional  structure.  The  structure  is 
assumed  to  be  free  to  rotate  about  its  longitudinal  axis,  and  to  trans¬ 
late  in  its  principal  directions. 

The  structure  is  considered  to  be  composed  of  a  series  of  planar 
bents.  In  each  bent  the  frame  and  shear-wall  elements  are  lumped  into 
an  analytical  model  composed  of  a  single  equivalent  frame  and  single 
equivalent  shear-wall.  In  each  individual  bent,  frame  and  shear-wall 
elements  may  be  coupled  or  uncoupled. 

The  structure  is  analyzed  under  lateral  loads  applied  at  each  floor 
level.  A  stiffness  approach  is  used  to  develop  the  equilibrium  equa¬ 
tions.  These  equations  are  then  solved  by  an  iterative  procedure,  which 
begins  by  assuming  locations  for  the  centres  of  rotation  of  the  structure. 
These  locations  are  then  adjusted  until  the  structure  is  in  equilibrium. 

The  validity  of  the  method  is  checked  against  other  techniques. 

The  effect  of  torsional  action  on  the  shear  distribution  within  the  struc¬ 
ture  is  illustrated  by  a  series  of  analyses  of  a  ten  story  structure. 

The  results  are  tabulated,  and  discussed  in  detail.  The  limitations  of 
the  analysis  are  indicated  and  possible  extensions  are  discussed. 
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CHAPTER  I 


INTRODUCTION 

Due  to  the  complexity  of  many  multi -story  structures,  it  has  been 
common  practice  to  analyze  the  three-dimensional  system  as  a  series  of 
braced  and  unbraced  planar  bents.  Within  each  bent  members  are  designed 
to  carry  the  directly  applied  loads,  with  some  modification  for  the  inter¬ 
action  between  adjacent  bents. 

The  actual  structure,  however,  deforms  in  a  three  dimensional  manner 
in  order  to  develop  internal  forces  which  are  in  equilibrium  with  the 
applied  loads.  The  manner  in  which  the  structure  reacts  depends  on  the 
arrangement  of  the  various  stiffening  elements,  their  relative  stiffnesses, 
the  interaction  between  adjacent  elements,  as  well  as  the  type  and  point 
of  application  of  the  loading.  Within  the  individual  bents,  the  normal 
shear- type  of  frame  deformation  is  modified  by,  and  must  be  compatible 
with,  the  cantilever  deformation  occurring  in  the  shear-walls. 

Computer  solutions  have  been  presented  for  determining  the  shear 
distribution  between  frame  and  shear-wall  elements  of  planar  bents  (1,  2). 
These  solutions  are  satisfactory  when  the  arrangement  of  the  stiffening 
elements  within  a  structure  is  symmetric  and  the  applied  lateral  load  at 
each  floor  level  passes  through  the  centre  of  resistance.  Under  the 
above  conditions  the  floor  diaphragms  do  not  rotate  and  the  deflections 
of  the  stiffening  elements  are  equal  at  each  floor  level. 

If  a  torque  is  applied  about  the  longitudinal  axis  of  the  structure; 


1 


I  P37-WO 


i  au  tc  \/  -'t?-  rd  r  JiH  v.ri6n  o  y  xyiqT'.v.  of  t  o)  ujG 

9«^lT9b  o'i  •  1  d  r!.  69  n  rr:1  H  .  sdrrsd  h  so  go  dm  bflfi  boosud 

-rt9dnr  art*  io'l  »  (ribo^rbom  ms  riJyv  teb6of  bsKqqfi  yldoanib  ertj  ynso  od 

-  n  - 

.sbfeor  bsrfqqs 

dnroq  bni5  oqy}  9fii  \t-  f  f9w  ae  ,  afro*  ifi  met  ft  n  s  nordoBT.Mr  ';• 

. 

. rsv9r  -roof*  rtoes  J6  feups  90.  aJ  9:’  f9  onffis^  riii  9rtJ  ^0 
;yiuJo:/  \  i  /e  f  nrbudrpnor  9fi7  d .<ods  be r f qq/..  .  y  ?od  b  ^1 

'' 


2 


or  if  the  layout  of  the  stiffening  elements  is  asymmetric,  the  structure 
will  rotate  as  well  as  translate  (3,  4,  5).  The  interaction  between 
shear-wall  and  frame  elements  in  each  bent  is  altered  considerably  by 
the  torsional  action  of  the  structure  (6).  The  shear  distribution  among 
individual  bents,  as  well  as  between  shear-wall  and  frame  elements  in 
each  bent,  may  be  quite  different  from  that  predicted  when  rotation  of 
the  structure  is  ignored. 

The  three  dimensional  aspect  of  structural  behavior  forms  the  basis 
of  the  present  investigation.  A  first-order  elastic  analysis  is  developed 
which  is  used  to  investigate  the  behavior  of  structures  free  to  rotate 
about  a  longitudinal  axis,  as  well  as  translate  under  the  action  of  lat¬ 
eral  loads.  The  analysis  is  used  to  examine  the  shear  distribution  in 
various  example  structures  for  several  arrangements  of  stiffening  elements, 
and  the  results  are  presented  in  a  systematic  manner.  Comparison  is  made 
between  this  and  other  methods  of  analysis.  The  limitations  of  the  present 
method  are  discussed  and  modifications  suggested  to  extend  the  application 
into  the  inelastic  range. 
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CHAPTER  II 


PREVIOUS  INVESTIGATIONS 

It  has  been  common  practice  to  simplify  the  analysis  of  a  multi¬ 
story  structure  by  dividing  it  into  a  series  of  two-dimensional  plane 
frames.  The  lateral  loads  were  distributed  to  the  bents  in  the  structure 
on  the  basis  of  their  tributary  areas.  The  effect  of  the  torsional  de¬ 
formations  of  the  structure  under  the  action  of  these  lateral  loads  was 
neglected  and  only  translational  movements  (in  both  principal  directions) 
were  considered. 

If  the  arrangement  of  bents  within  a  structure  is  not  symmetrical, 
the  above  assumptions  are  erroneous.  In  this  case  the  structure  will 
deform  so  that  certain  bents  are  subjected  to  larger  deflections  than 
are  others  and  must  carry  a  correspondingly  greater  share  of  the  lateral 
load. 

In  the  last  decade  high  intensity  earthquakes  have  occurred  in 
areas  containing  contemporary  multi-story  structures  (7,  8,  9).  Many 
of  the  structures  failed  in  a  torsional  mode.  The  damage  to  these  build¬ 
ings  has  emphasized  the  necessity  of  considering  the  effect  of  torsion 
in  design. 

A  completely  general  method  of  analysis  is  not  available  for  asym¬ 
metric  frame  shear-wall  structures.  In  Canada  (10),  however,  torsional 
effects  must  be  considered  when  designing  for  earthquake  loading.  Because 
of  the  uncertainties  in  the  analysis,  the  eccentricities  computed  for 
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asymmetric  structures  are  increased  by  50  per  cent  and,  in  fact,  addi¬ 
tional  insurance  factors  added. 

In  Californian  earthquake  codes  (11)  provision  is  also  made  for  an 
increase  in  shear  due  to  torsion.  The  torque  is  computed  as  the  product 
of  the  lateral  earthquake  loading  and  the  distance  separating  the  centre 
of  mass  and  the  centre  of  rigidity.  It  is  recommended  that  negative  tor¬ 
sional  shears  be  neglected  since  the  earthquake  motion  may  have  any  dir¬ 
ection  of  propagation.  It  is  further  stipulated  that  if  the  vertical 
elements  depend  on  the  diaphragm  action  of  the  floor  system  for  shear 
distribution,  the  minimum  torsional  moment  shall  be  the  product  of  the 
story  shear  and  an  eccentricity  equal  to  5  per  cent  of  the  maximum  build¬ 
ing  width  at  that  level. 

A  design  method  proposed  by  the  Portland  Cement  Association  (12)  pro¬ 
vides  for  the  possibility  of  torsional  loading  even  when  the  structure 
is  symmetrical.  In  this  case  a  minimum  design  torque  is  specified.  For 
asymmetric  structures,  the  centre  of  rotation  for  each  floor  is  assumed 
to  be  located  at  the  centre  of  rigidity,  A  torsional  moment  is  applied 
to  the  structure,  which  is  equal  to  the  product  of  the  lateral  load  and 
the  distance  separating  the  centre  of  mass  and  the  centre  of  rigidity. 

The  above  design  codes  specify  the  loadings  to  be  used  but  do  not 
provide  methods  to  determine  the  distribution  of  the  resisting  shears  to 
the  bents  in  a  structure.  The  conservative  nature  of  these  provisions 
emphasizes  the  uncertainty  surrounding  the  three-dimensional  aspects  of 
structural  behavior. 

Analyses  of  unsymmetric  single  story  shear  wall  structures  have 
been  presented  (4,  13).  These  methods  take  into  consideration  the  shear 
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contribution  to  the  overall  deflection  of  the  structure.  The  analyses 
assume  that  the  centre  of  rotation  is  at  the  centre  of  resistance  of  the 
structural  elements.  A  torsional  load  is  applied  which  is  equal  to  the 
product  of  the  lateral  load  and  its  eccentricity.  The  total  load  on  any 
particular  element  is  then  obtained  by  superimposing  the  lateral  and 
torsional  load  contributions.  This  approach  is  satisfactory  for  single 
story  structures,  where  shear  walls  are  not  coupled  to  frames.  Extension 
of  the  method  to  multi -story  structures  would  be  difficult  because  of  the 
interaction  of  the  various  elements  and  the  carry-over  effects  between 
adjacent  stories. 

Wilbur  (3)  presented  equations,  for  determining  the  shear  distribu¬ 
tion  to  the  individual  bents  of  a  structure.  These  equations  allow  for 
the  torsional  deformations  which  arise  from  the  lack  of  coincidence  of 
the  point  of  load  application  and  the  centre  of  rigidity  of  the  bents. 

The  equations  are  based  on  the  assumption  that  points  of  inflection  occur 
at  the  mid  points  of  all  members  and  that  all  joint  rotations  are  equal 
at  a  particular  floor.  The  method  is  further  limited  to  bents  containing 
only  uncoupled  frames. 

Blume,  Newmark  and  Corning  (12)  suggest  a  method  which  is  basically 
the  same  as  that  presented  by  Wilbur.  In  both  cases  the  rotation  of  the 
beam  to  column  joints  are  considered  but  no  allowance  is  made  for  the 
interaction  between  frame  and  shear  wall  elements. 

Khan  and  Sbarounis  (14)  have  presented  an  analysis  of  planar  bents 
which  considers  the  interaction  between  frame  and  shear-wall  elements. 

The  method  lumps  the  structure  into  an  equivalent  bent  which  is  composed 
of  two  parts:  a  frame  element  and  a  shear  wall  element.  The  entire 
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lateral  load  is  applied  to  the  wall  element  and  its  deflected  position 
calculated.  The  loads  necessary  to  make  the  frame  element  deform  to  the 
same  position  are  then  calculated.  The  load  on  the  wall  is  modified  until 
the  frame  and  shear  wall  are  compatible  and  in  equilibrium  under  the 
applied  lateral  load.  If  the  distribution  of  frame  and  shear  wall  elements 
is  not  symmetrical  in  the  structure  a  torsional  analysis  is  suggested.  In 
this  analysis  the  centre  of  rotation  is  to  be  located  on  the  basis  of  the 
shear  resisted  by  each  element  and  its  location  in  the  structure.  Since 
the  analysis  considers  interaction  effects,  the  centre  of  rotation,  computed 
on  this  basis  no  longer  coincides  with  the  centre  of  rigidity  computed 
simply  on  the  basis  of  El. 

Winokur  and  Gluck  (6)  have  formulated  the  analysis  of  asymmetric 
structures,  containing  frame  and  shear-wall  elements,  in  matrix  form. 

Three  degrees  of  freedom  are  considered  for  each  floor;  translation  in  the 
X  and  Y  direction  and  rotation  about  the  vertical  axis  of  the  structure, 

0.  The  analysis  consists  of  the  following  three  steps: 

1.  The  overal  lateral  stiffness  matrix  for  the  structure  is  first 
evaluated.  This  is  the  sum  of  the  lateral  stiffness  matrices 
for  each  of  the  stiffening  elements.  The  overall  lateral  stiff¬ 
ness  matrix  for  the  structure,  |K|,  is  a  3n  by  3n  square  matrix, 
where  n  is  the  number  of  floors  in  the  structure. 


|  Kxx  | 

|Kxy| 

|  Kx0 1 

|Kyx| 

|Kyy| 

|Ky0| 

|K0x| 

1  Key  | 

|  Kee  | 
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where  each  of  the  submatrices  represent  the  stiffness  of  the 
structure  in  the  direction  of  the  first  subscript,  due  to  a 
displacement  in  the  direction  of  the  second  subscript. 

2.  Next,  the  three  equations  of  equilibrium  are  formulated  at 
each  floor  level.  The  solution  of  these  equations  involves 
the  inversion  of  a  3n  by  3n  matrix  and  determines  the  deforma¬ 
tions  of  the  structure. 

3.  The  member  forces  are  obtained  by  multiplying  the  lateral 
stiffness  submatrices  for  each  element  by  the  deformations  ob¬ 
tained  in  step  2. 

A  method  is  not  specified  for  calculating  the  stiffness  matrix  of 
the  individual  elements  in  the  first  step.  Reference  is  made  to  the  work 
of  Khan  and  Sbarounis  (14),  and  Clough,  King  and  Wilson  (2). 

Weaver  and  Nelson  (5)  have  developed  an  analysis  which  accounts 
for  rotation  of  the  floors  as  well  as  translation  in  the  principal  dir¬ 
ections  of  the  structure.  The  method  is  limited  to  framed  structures 
and  shear  wall  elements  are  not  considered.  The  method  consists  of  first 
determining  member  stiffness  matrices.  As  both  axial  shortening  and  tor¬ 
sion  are  considered  the  member  stiffness  matrices  consist  of  6  by  6 
arrays  for  beams  and  12  by  12  arrays  for  columns.  The  geometric  rela¬ 
tions  of  the  members  within  the  structure  are  then  used  to  transform  the 
member  stiffness  matrices  to  a  floor  stiffness  matrix.  The  equilibrium 
equations  for  the  structure  are  formulated  for  each  pair  of  adjacent 
floors.  The  equations  contain  joint  rotations  for  both  floors  and  floor 
translations  and  rotations  for  all  floors  in  the  structure,  as  unknowns. 
These  unknowns  are  reduced  by  a  forward  elimination  process  starting  at 
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the  top  of  the  structure.  Joint  rotations  and  floor  displacements  are 
zero  at  the  base,  thus  the  first  floor  displacements  are  obtained  by 
back  substitution.  Member  forces  are  determined  by  multiplying  the  joint 
displacements  by  the  member  stiffness  matrices. 

As  a  summary  of  the  above  literature  survey;  it  is  apparent  that, 
there  is  no  completely  general  method  for  analysing  multi-story  frame 
shear-wall  structures.  Approximate  design  methods  for  framed  structures 
have  been  developed  based  on  simplifying  assumptions  to  reduce  the  number 
of  unknowns.  These  methods  are  approximations  at  best  and  do  not  make 
any  allowance  for  the  interaction  between  individual  elements. 

Complex  exact  analyses  are  available  for  analysing  framed  structures 
but  these  are  not  applicable  to  the  large  number  of  structures  that  in¬ 
corporate  shear  walls.  An  analysis  has  been  presented  for  considering 
planar  bents  incorporating  frame  and  shear-wall  elements,  but  this  is  not 
suitable  for  three  dimensional  structures  subject  to  torsional  loads.  Even 
for  symmetric  structures  not  subjected  to  torsion,  the  techniques  of  lump¬ 
ing  the  structure  into  a  planar  bent  introduces  errors  of  unknown  magnitude. 
A  method  is  proposed  for  a  three  dimensional  analysis  of  framed  shear  wall 
structures  but  this  is  not  applicable  to  the  case  where  the  frame  and 
shear-walls  are  coupled. 

There  is  an  immediate  need  for  a  method  of  analysis  that  can  be 
applied  to  structures  containing  coupled  frame  shear-wall  bents  as  well  as 
uncoupled  frame  and  shear-wall  elements.  This  method  must  be  applicable 
to  large  multi -story  structures  of  varying  geometric  layouts. 
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CHAPTER  III 


ANALYTICAL  MODEL 

In  the  method  of  analysis  discussed  in  this  and  the  succeeding 
section,  it  will  be  assumed  th*at  the  structure  is  composed  of  a  linear 
elastic  material  and  that  the  response  of  each  member  can  be  predicted 
by  a  linear  analysis.  The  loads  are  lateral  loads,  concentrated  at  each 
floor  level,  and  assumed  to  be  applied  statically.  In  addition,  the  dia¬ 
phragms  at  each  floor  level  are  assumed  to  be  infinitely  rigid  in  their 
own  plane. 

The  motion  of  a  typical  floor  diaphragm  is  shown  in  FIGURE  3.1. 

The  lateral  loads  cause  translation  as  well  as  a  rotation  of  the  floor 
diaphragm,  so  that  general  point  A  moves  to  a  new  position  A'.  If  it 
is  assumed  that  the  deformations  are  small  relative  to  the  overall  dim¬ 
ensions  of  the  structure,  then  bent  k  ,  which  spans  in  the  Y  direction, 

A 

will  undergo  a  translation,  v^,  in  the  Y  direction,  where 

vk  =  va  +  RXk  x  6  (3-l) 

In  Equation  (3.1),  vk  is  the  total  deflection  of  bent  kx  in  the  Y  dir¬ 
ection,  vA  is  the  translation  of  the  general  point  A  in  the  Y  direction, 

RX.  is  the  co-ordinate  of  bent  k  with  respect  to  point  A,  measured  per- 

K  X 

pendicular  to  the  bent,  and  6  is  the  rotation  of  the  floor  diaphragm 
about  the  general  point  A.  From  Equation  (3.1)  the  difference  in  deflec¬ 
tion  between  adjacent  bents  is  equal  to  the  product  of  the  angle  of 
rotation  and  the  perpendicular  distance  between  the  bents.  Similarly 
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for  bent  k^,  which  spans  in  the  X  direction,  the  deflection  in  the  X 
direction  is  given  by 


uk  =  uA  -  RYk  x  0  (3.2) 

where  uk  is  the  total  deflection  of  bent  in  the  X  direction,  is 
the  translation  of  the  general  point  A  in  the  X  direction,  RYk  is  the  co¬ 
ordinate  of  bent  kx  with  respect  to  point  A,  measured  perpendi cul ar  to 
the  bent.  In  Equations  (3.1)  and  (3.2),  the  rotation  is  assumed  to  be 
positive  if  clockwise  and  the  distances  from  the  general  point  A  to  the 
X  and  Y  bents,  RX  and  RY,  are  assumed  to  be  positive  as  shown  in  FIGURE 
3.1.  Similar  expressions  may  be  developed  for  each  bent,  k,  at  each 
floor  level ,  i . 

If  the  torsional  stiffness  of  the  beams  linking  adjacent  bents  is 
neglected,  the  behavior  of  the  structure  will  depend  on  that  of  the  in¬ 
dividual  planar  bents  as  shown  in  FIGURE  3.2.  The  deformation  of  a  parti¬ 
cular  bent  at  each  floor  level  differs  from  the  deformation  of  adjacent 
bents  due  to  their  different  locations  relative  to  the  general  point  A. 

Each  bent  in  the  system  may  contain  frames  coupled  to  one  or  more  shear 
walls. 

To  obtain  the  analytical  model  used,  the  individual  walls  and  frames 
must  be  lumped  within  each  bent.  The  lumping  technique  has  been  used 
previously  for  the  analysis  of  symmetrical  structures  consisting  of  shear 
walls  and  frames  (1).  In  each  story  of  the  lumped  model,  each  joint  has 
one  rotational  degree  of  freedom  and  the  story  has  one  translational  de¬ 
gree  of  freedom.  Each  degree  of  freedom  requires  one  equilibrium  equation 
in  the  solution.  It  is  assumed  that  within  a  particular  bent  the  beam 
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to  column  connections  undergo  equal  joint  rotations  and  the  beam  to  wall 
connections  undergo  equal  joint  rotations  (but  different  from  those  at 
the  beam  to  column  connections);  at  a  particular  floor  level.  This  re¬ 
duces  the  number  of  unknowns  to  three  per  story  for  the  case  of  a  coupled 
frame  shear-wall  bent,  and  to  two  for  an  uncoupled  frame  bent.  A  typical 
coupled  frame  shear-wall  bent  is  shown  in  FIGURE  3.3A  and  the  equivalent 
model  in  FIGURE  3.3B  .  FIGURES  3.4A  and  3.4B  depict  an  uncoupled 

frame  bent,  and  the  equivalent  model.  In  the  model  the  wall  stiffness, 

KW,  is  equal  to  the  sum  of  stiffnesses  of  all  the  walls  in  a  given  story. 

The  column  stiffness,  KC,  in  the  model,  is  equal  to  the  sum  of  stiffnesses 
of  all  the  columns  of  the  story.  The  beams  are  divided  into  two  groups; 
those  linking  the  shear  walls  to  the  columns  are  referred  to  as  wall  beams 
and  those  linking  individual  columns  are  referred  to  as  frame  beams.  In 
the  analytical  model,  the  wall  beam  has  a  stiffness  equal  to  the  sum  of 
the  wall  beam  stiffnesses  and  the  frame  beam  has  a  stiffness  double  that 
of  the  sum  of  the  frame  beam  stiffnesses.  One  end  of  the  frame  beam  is 
placed  on  a  roller  and  the  joint  rotation  of  the  roller  end  is  considered 

to  be  equal  to  that  at  the  column  end.  The  distinction  between  wall  beam 

and  frame  beam  is  necessary  because  of  the  finite  width  of  the  wall.  This 
causes  a  vertical  deflection  of  the  wall  beam,  at  the  wall  joint,  due  to 
the  wall  rotation.  This  effect  increases  the  bending  moment  in  the  wall 
beam  for  a  given  wall  joint  rotation  and  is  similar  to  an  increased  rota¬ 
tional  restraint  on  the  wall  at  each  floor  level. 

The  final  stage  in  developing  the  analytical  model  is  to  reduce  the 
number  of  degrees  of  freedom  of  each  floor  diaphragm,  from  three  to  one. 
This  is  achieved  by  assuming  that  each  floor  diaphragm  rotates  about  a 
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point  in  its  plane  as  shown  in  FIGURE  3.5.  This  point,  referred  to  as 
the  centre  of  rotation,  has  a  unique  location,  relative  to  the  assumed 
general  point  A  in  FIGURE  3.1,  such  that  the  following  geometric  condi¬ 
tions  are  satisfied, 


vA  =  (X0  -  XA)  X  9 


(3.3) 


UA  ■  (Yo  *  V  x  0 


(3.4) 


where  v^,  u^  and  6  have  been  previously  defined,  and  are  the  co¬ 
ordinates  of  point  A  and  X  and  Y  are  the  co-ordinates  of  the  centre  of 

oo 

rotati on . 

To  initiate  the  analysis,  the  position  of  the  centre  of  rotation  must 

be  assumed,  relative  to  a  reference  point.  The  two  unknowns,  X  and  Y  , 

o  o 

are  thus  eliminated.  This  assumed  location  is  varied  until  the  angle  of 
rotation  selected  results  in  a  balance  between  the  applied  loads  and  the 
resisting  forces,  in  both  the  X  and  Y  directions;  and  a  torque  which  is 
equal  to  zero  about  any  point  in  space.  This  type  of  iterative  solution 
requires  a  relatively  small  amount  of  computer  storage  space  and  thus 
more  complex  structures  may  be  analyzed.  The  iteration  process  itself, 
however,  is  a  complex  portion  of  the  program  and  requires  a  significant 
amount  of  computer  calculation  time. 
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FIGURE  3.1  TYPICAL  FLOOR  DIAPHRAGM 


FIGURE  3.2  TWO  DIRECTIONAL  SYSTEM  OF  PLANAR  BENTS 
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FIGURE  3, 3A  COUPLED  BENT 


FIGURE  3.3B  LUMPED  COUPLED  BENT 
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KC  =  (KC1  +  KC2  +  KC3) 


FIGURE  3.4B  LUMPED  UNCOUPLED  BENT 
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FIGURE  3.5  FLOOR  DIAPHRAGM  ROTATION 


CHAPTER  IV 


METHOD  OF  ANALYSIS 

The  present  analysis  is  based  on  the  simplified  model  described  in 
the  previous  chapter  and  involves  two  distinct  steps;  first  a  lateral 
stiffness  matrix  is  developed,  in  terms  of  the  story  sway  deformations, 
for  each  of  the  bents  in  the  structure  (in  both  the  X  and  Y  directions); 
secondly  a  rotational  stiffness  matrix  is  computed  for  the  entire  struc¬ 
ture.  This  is  possible  because  the  story  sway  vectors  of  each  bent  are 
related,  since  the  floor  diaphragms  are  assumed  to  be  infinitely  rigid  in 
their  own  planes.  The  rotational  stiffness  matrix  is  then  inverted  and 
multiplied  by  the  applied  torque  to  obtain  the  floor  diaphragm  rotation 
vector.  By  back  substitution  the  resisting  floor  forces  are  calculated 
(in  both  the  X  and  Y  directions).  If  the  resisting  forces  and  applied 
loads  are  not  in  balance,  an  iterative  process  is  used  to  relocate  the 
assumed  positions  of  the  centres  of  rotation.  This  iterative  procedure 
is  repeated  until  the  resisting  floor  forces  are  in  balance  with  the  applied 
loads  and  the  moment  about  any  point  in  space  is  equal  to  zero,  for  all 
floors  in  the  structure. 

The  procedure  used  to  determine  the  lateral  stiffness  matrix  for 
each  bent  depends  on  whether  the  bent  consists  of  a  coupled  frame  and 
shear  wall  or  an  uncoupled  frame.  The  case  of  a  coupled  frame  and  shear 
wall  will  be  treated  in  detail  and  the  modifications  necessary  for  an 
uncoupled  frame  will  be  outlined  briefly. 

Consider  the  single  bent  shown  in  FIGURE  4.1.  As  lateral  loads  are 
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applied,  the  bent  deforms  as  shown  in  FIGURE  4.2.  In  each  story  the 
sway  displacements  of  the  frame  and  shear  wall  are  equal.  Because  of  the 
finite  width  of  the  shear  wall,  the  wall  beam  undergoes  a  vertical  dis¬ 
placement  at  the  connection  to  the  wall.  The  end  rotations  of  each  frame 
beam  are  assumed  to  be  equal. 

For  each  floor  level  there  are  three  unknowns;  A,  0  and  0  .  The 

c  w 

lateral  stiffness  matrix  is  developed  by  writing  the  equilibrium  equations 
in  terms  of  the  slope  deflection  coefficients  and  expressing  the  joint 
rotations  in  terms  of  the  story  sway  deformations.  The  elimination  tech¬ 
nique  used  to  express  joint  rotations  as  story  sway  displacements  is  de¬ 
scribed  in  APPENDIX  A.  The  resulting  floor  force  and  joint  rotation  vec¬ 
tors  are  as  follows: 


V 

0 

0 


'k 

I c 
!k 

iW 

■k 


k  X  Nk 

(4.1a) 

kc*  Nk 

(4.1b) 

<3 

X 

(4.1c) 

where  |K|£,  |K|®c  and  |K|^W  are  square  matrices  which  depend  only  on  the 
bent  member  properties  and  their  geometry.  The  story  sway  vector,  | A | ^ , 
contains  unique  values  for  each  bent.  These  must  be  determined  before 
the  joint  rotation  vectors  can  be  obtained  by  back  substitution  in 
Equations  (4.1b)  and  (4.1c). 

The  next  step  in  the  analysis  is  to  determine  the  story  sway  vector 
| A | .  .  The  lateral  displacement  of  a  given  floor  in  a  particular  bent  is 
equal  to  the  product  of  the  floor  diaphragm  rotation  and  the  perpendicular 
distance  between  the  bent  and  the  centre  of  rotation  of  the  floor  diaphragm. 
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For  a  bent  spanning  in  the  Y  direction: 

Sk,i  =  ei  x  (x0,i  ‘  Xk,i)  (4-2a) 

where  .  is  the  displacement  (from  the  undeformed  oosition)  of  bent  k 

at  floor  level  i,  0.  is  the  rotation  of  floor  diaphragm  i  about  its 
centre  of  rotation  and  (X  .  -  X.  .)  is  the  co-ordinate  distance  of  bent 

U  )  I  K  5  I 

k  from  the  centre  of  rotation  of  the  ith  floor  diaphragm.  The  sway  defor- 

t  h 

mation  in  the  iLn  story  of  bent  k  can  then  be  expressed  as. 


Ak,i  Sk,i  "  Sk,i-1 


(4.2b) 


therefore 


,1  Vx  (Xo,i  ‘  Xk,i)  '  ei-l  x  (Xo,i-l  '  xk,1-l^  (4-2c) 


J.  L. 

where  A,  .  is  the  sway  deformation  of  the  l  n  story  of  bent  k.  This  may 

K  ,1 

be  written  in  general  matrix  form  as  follows: 


Nk  =  |D|k  x  1 0 1  (4. 2d) 

where  |A|k  is  the  story  sway  vector  for  bent  k,  |D|k  is  a  two-diagonal 
square  matrix  in  which  each  row,  i,  contains  the  perpendicular  distances 
(X  •  -i  -  X.  .  ,),  (X  .  -  X.  .),  between  the  bent  and  the  centres  of  ro- 

0,1-1  K  ,1  -  I  0,1  K ,  l 

tation  of  two  adjacent  floor  diaphragms,  and  |e|  is  the  floor  diaphragm 
rotation  vector.  The  floor  force  vector  |V|k  may  now  be  written  as, 

|vlk  -  lKlkx  I® 


(4.3a) 
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where 

|K|“  =  |K|£  x  | D | k  (4.3b) 

The  next  step  in  the  analysis  is  to  formulate  the  torsional  equilibrium 
equation.  A  reference  point  is  chosen  at  each  floor  level  so  that  all  re¬ 
ference  points  lie  on  a  common  vertical  axis  and  all  of  the  bents  lie  in 
the  first  quadrant  of  the  co-ordinate  system,  as  shown  in  FIGURE  4.3.  The 
building  is  oriented  so  that  the  centre  of  rotation  lies  on  the  positive 
X  side  of  the  centre  of  resistance.  This  will  result  in  positive  rotations 
of  the  floor  diaphragms  and  resisting  floor  forces  which  are  in  opposition 
to  lateral  loads  applied  in  the  positive  X  and  Y  directions.  The  resisting 
torques,  about  the  reference  point,  developed  by  the  bents  spanning  in  both 
the  X  and  Y  directions  are  computed  individually  and  then  summated.  Written 
in  general  matrix  form, 

NX  n 

|T|X  =  l  (|K|°  x  |X|k)  (4.4a) 

k“  1 

NY 

| T | v  =  l  ( | K | P  x  | Y | . )  (4.4b) 

where  |T|^  is  a  torsional  stiffness  matrix  which  includes  all  the  bents 
in  the  X  direction,  and  | T | y  is  the  torsional  stiffness  matrix  for  all 
the  bents  in  the  Y  direction.  |X|k  is  the  vector  representing  the  dis¬ 
tance  of  bent  k  (in  the  X  direction)  from  the  reference  point,  |Yjk  is 
the  vector  representing  the  distance  of  bent  k  (in  the  Y  direction)  from 
the  reference  point,  NX  is  the  number  of  bents  spanning  in  the  X  direction 
and  NY  represents  the  number  of  bents  spanning  in  the  Y  direction.  The 
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total  torsional  stiffness  matrix,  |T|R,  for  the  entire  structure  may  now 
be  written  as, 

ItIr  =  lTlx  +  I T I Y  +  I T I SV 

where  |T|<.y  is  the  St.  Venant  torsional  stiffness  matrix  for  the  structure. 
This  represents  the  restraint  offered  to  the  floor  diaphragm  rotation  by 
the  St.  Venant  stiffness  of  the  vertical  stiffening  elements.  The  torque, 
computed  about  the  reference  point,  due  to  the  applied  lateral  loads  is 
given  by, 

| T | p  =  | F | x  x  |X|H  -  | F | y  x  |Y|h  (4.6) 

where  |T|p  is  the  torque  vector,  |F|X  and  |F|y  are  the  applied  load  vectors 
in  the  positive  Y  and  X  directions  respectively,  |X|H  and  |Y|H  are  vectors 
representing  the  distance  between  the  reference  point  and  the  point  of 
application  of  the  applied  load  vectors  | F | x  and  |F|y  respectively.  As 
the  resultant  torque  about  the  reference  point  must  be  equal  to  zero, 

| T | F  +  |T|r  x  | e |  =0 

and  therefore 

M  -  -  ItIr1  x  |t|f 

The  assumed  positions  of  the  centres  of  rotation  are  now  checked  by 
substituting  for  |e|  in  Equation  (4.3a).  This  step  will  determine  the 
resisting  floor  forces  for  each  bent.  The  total  resisting  forces  in  the 
X  and  Y  directions  are  obtained  by  summing  the  appropriate  forces  for  the 


(4.7a) 

(4.7b) 
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individual  bents.  If  the  resisting  forces  in  the  X  and  Y  directions  are 
equal  to  the  applied  loads  then  the  initial  assumed  locations  of  the 
centres  of  rotation  are  correct.  If  the  resisting  forces  do  not  balance 
the  applied  loads  then  a  correction  technique  is  applied  which  shifts  the 
centres  of  rotation  in  both  the  X  and  Y  directions  until  the  structure  is 
in  equilibrium. 

The  correction  is  applied  by  first  determining  an  influence  coef¬ 
ficient  matrix  for  the  shift  of  the  centres  of  rotation  in  the  X  direction. 
The  Y  co-ordinates  of  the  centres  of  rotation  are  held  constant.  The  in¬ 
fluence  coefficient  matrix  is  obtained  by  increasing  the  distance  from  the 
reference  point  to  the  centre  of  rotation  for  the  first  floor  by  1  inch, 
while  maintaining  the  position  of  the  centres  of  rotation  for  all  other 
floors  at  the  initial  values.  With  this  new  arrangement  of  the  centres 
of  rotation,  the  torsional  analysis  described  above  is  repeated  and  the 
resultant  floor  forces  are  calculated.  The  influence  coefficient  for  the 
first  floor  is  equal  to  the  difference  in  total  floor  force  between  the 
initial  value  and  that  obtained  for  the  shifted  centre  of  rotation.  This 
procedure  is  repeated  for  all  floors;  in  each  case  the  centres  of  rotation 
are  held  at  the  initial  values  except  for  the  floor  under  consideration. 

The  physical  interpretation  of  the  process  is  that  if  the  centre  of  rota¬ 
tion  of  a  particular  floor  is  shifted  by  x  inches,  then  the  resisting  floor 
force  at  any  floor  is  changed  by  the  nroduct  of  x  and  the  influence  coef¬ 
ficient  for  that  floor.  Consequently,  if  the  resisting  floor  forces  for 
the  initial  trial  differ  from  the  applied  loads  it  is  possible  to  deter¬ 
mine  the  corrections  to  the  positions  of  the  centres  of  rotation,  in 
order  for  the  resisting  floor  forces  to  approach  the  applied  lateral  loads. 
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This  may  be  written  in  general  matrix  form  as, 

I  AX  I  =  Id'1  X  (- 1  F|x  -  |V|X)  (4.8) 

where  | AX |  is  the  vector  representing  the  corrections  to  be  applied  to 
the  positions  of  the  centres  of  rotation  (in  the  X  direction),  | C | ^  is  the 
transpose  of  the  influence  coefficient  matrix  (for  a  shift  of  the  centres 
of  rotation  in  the  X  direction),  and  |V|X  is  the  resisting  floor  force 
vector  in  the  Y  direction,  for  the  initial  position  of  the  centres  of 
rotation.  Due  to  the  manner  in  which  the  iterative  solution  is  formulated, 
the  influence  coefficients  are  non-linear.  The  correction  vector,  | AX | ,  in 
Equation  4.8  does  not  immediately  lead  to  the  correct  position  for  the 
centres  of  rotation  and  several  iteration  cycles  may  be  required.  After 
each  cycle  the  corrected  position  of  the  centres  of  rotation  become  the 
initial  values  for  the  next  cycle.  When  the  floor  forces  in  the  Y  direc¬ 
tion  balance  the  applied  loads  within  a  specified  limit,  a  similar  correc¬ 
tion  technique  is  applied  to  shift  the  centres  of  rotation  in  the  Y  dir¬ 
ection.  In  matrix  form, 

| AY |  =  ICl'1  x  (-| F|y  -  |V | y)  (4.9) 

where  | AY |  is  the  vector  representing  the  corrections  to  be  applied  to 
the  position  of  the  centres  of  rotation  (in  the  Y  direction),  |C|y  is  the 
transpose  of  the  influence  coefficient  matrix  (for  a  shift  of  the  centres 
of  rotation  in  the  Y  direction)  and  |V|y  is  the  vector  of  the  resisting 
floor  forces  in  the  X  direction  (for  the  initial  position  of  the  centres 
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of  rotation).  The  number  of  iteration  cycles  required  to  determine  the 
correct  position  of  the  centres  of  rotation  will  depend  on  how  correct 
the  initial  choices  were.  If  the  building  is  symmetrical  about  a  centroi- 
dal  axis  and  there  are  no  applied  loads  in  this  direction  then  the  centre 
of  rotation  will  lie  on  this  axis.  This  factor  facilitates  the  choice  of 
the  initial  position  of  the  centres  of  rotation. 

A  computer  program  was  written  in  Fortran  IV  language  for  use  on  an 
IBM  360/67  system  in  order  to  facilitate  the  analysis.  A  print-out  of 
the  program  is  included  as  APPENDIX  C.  In  this  particular  program  the 
structure  size  was  limited  to  twenty  stories,  with  a  maximum  of  four  bents 
in  the  X  direction  and  three  bents  in  the  Y  direction.  The  size  of  struc¬ 
ture  that  may  be  analyzed,  however,  is  limited  only  by  the  computer  stor¬ 
age  capacity. 

The  computer  program  has  three  distinct  divisions.  COMPAK  reads  in 
the  member  properties  and  computes  the  bent  stiffness  matrices.  TORPAK 
performs  the  torsional  analysis  until  an  equilibrium  condition  is  reached. 
CALPAK  uses  the  story  sway  and  joint  rotation  vectors  calculated  by 
TORPAK  to  compute  the  member  forces  for  each  bent.  The  procedure  followed 
by  the  program  is  described  in  greater  detail  below. 

The  COMPAK  stage  of  the  program  first  determines  whether  the  bent 
includes  a  shear  wall.  It  then  reads  in  the  member  properties  and  bent 
geometry.  From  these  quantities  the  member  stiffness  submatrices  are 
computed.  If  the  bent  consists  solely  of  a  frame,  four  submatrices  are 
required  but  if  the  frame  is  coupled  to  a  shear  wall  the  number  of  sub¬ 
matrices  increases  to  nine.  A  subroutine  then  calculates  the  bent  stiff¬ 
ness  matrix,  |K|^  ,  and  matrices  |K|®C  and  |K|®W  in  Equations  (4.1b)  and 
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(4.1c) 


The  torsional  stiffness  matrix,  |T|R  is  calculated  in  the  TORPAK 
stage.  The  location  of  bents  within  the  structure,  the  magnitude  of  the 
applied  lateral  loads,  their  points  of  application,  the  location  of  the 
centre  of  resistance  at  each  floor  level,  and  the  St.  Venant  torsional 
stiffness  for  each  floor  level  are  read  in  and  used  by  TORPAK.  TORPAK 
may  next  perform  a  torsional  analysis  or  an  analysis  with  rotation  of  the 
floor  diaphragms  suppressed. 

If  a  torsional  analysis  is  performed  the  torque  about  the  reference 
point,  due  to  the  applied  loads,  is  calculated.  If  the  centre  of  resist¬ 
ance  of  the  structure  is  within  5  per  cent  of  the  overall  width  of  the 
structure,  from  the  geometric  centre,  then  a  minimum  torque  is  applied. 
TORPAK  first  computes  |T|SV,  |T|X,  and  | T | Y ;  the  latter  two  will  differ 
for  each  change  in  the  position  of  the  centres  of  rotation.  The  initial 
position  of  the  centres  of  rotation  are  assumed  to  be  at, 


(4.10a) 


(4.10b) 


if  there  is  an  applied  lateral  load  in  both  the  X  and  Y  directions.  If 
there  is  no  applied  load  in  one  direction  the  initial  positions  of  the 
corresponding  centres  of  rotation  are  modified  to, 
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In  the  above  four  equations  X  .  is  the  distance  between  the  reference 

0,1 

point  and  the  centre  of  rotation  (in  the  X  direction),  Y  .  is  the  dis- 

0,1 

tance  between  the  reference  point  and  the  centre  of  rotation  (in  the  Y 
direction),  Xc  ^  is  the  distance  between  the  reference  point  and  the  centre 
of  resistance  (in  the  X  direction),  Y  .  is  the  distance  between  the  ref- 

C  ,  I 

erence  point  and  the  centre  of  resistance  (in  the  Y  direction),  all  taken 
t  h 

at  the  i  floor  level.  The  distance  between  the  reference  point  and  the 
X  bent  furthest  away  from  the  reference  point  is  denoted  as  X^x  ..  and 
Yny  is  the  distance  between  the  reference  point  and  the  Y  bent  furthest 
away  from  the  reference  point. 

The  floor  diaphragm  rotation  vector  |e|  is  calculated  by  inverting 

the  torsional  stiffness  matrix  |T|^„  The  resisting  floor  forces  in  both 

the  X  and  Y  directions  are  obtained  by  back  substitution  and  summation. 

A  check  is  performed  to  determine  whether  the  resisting  floor  forces  are 

within  1  per  cent  of  the  applied  lateral  loads.  If  so,  TORPAK  returns 

control  to  COMPAK,  after  calculating  the  story  sway  and  joint  rotation 

i  c  w 

vectors  | A | ^ ,  | 0 | ^  and  | 0 | ^  for  each  bent. 

If  the  resisting  forces  are  not  within  1  per  cent  of  the  applied  loads, 
TORPAK  calls  up  the  subroutine  which  calculates  the  influence  coefficients 
for  correcting  the  centres  of  rotation  in  the  X  direction.  The  correction 
cycle  is  repeated  until  the  resisting  forces  and  applied  loads  balance 
in  the  Y  direction.  Then  the  subroutine  for  correcting  the  position  of 
the  centres  of  rotation  in  the  Y  direction  takes  over.  After  each  cor¬ 
rection  cycle,  if  the  forces  and  loads  in  the  Y  direction  are  not  in  bal¬ 
ance,  these  are  balanced  first  before  correction  in  the  X  direction  con¬ 
tinues.  The  correction  procedure  is  continued  until  either  equilibrium 
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is  reached  or  the  number  of  cycles  exceeds  the  program  limit. 

COMPAK  then  calls  up  CALPAK  which  calculates  all  of  the  member  forces 
for  each  bent.  These  member  forces  are  based  on  the  last  values  of  I A I .  , 

1 6 1 ^  ar>d  |e|£  calculated  by  TORPAK.  Unless  the  structure  is  in  equili¬ 
brium,  therefore,  the  member  forces  calculated  by  CALPAK  are  not  neces¬ 
sarily  the  correct  values.  If  the  correction  procedure  is  terminated  by 
the  program  limit,  the  unbalance  of  forces  and  applied  loads  is  indicated 
by  a  print  statement  in  the  output, 

CALPAK  also  includes  the  write  subroutine,  which  prints  out  the 
member  forces,  before  returning  to  COMPAK  and  stopping,  A  flow  diagram 
of  the  program  is  presented  in  FIGURES  4,4,  4,5  and  4,6,  In  APPENDIX  B 
a  list  of  all  subroutines  is  presented  together  with  an  explanation  of 
the  tasks  that  each  subroutine  performs.  Also  included  in  APPENDIX  B 
is  a  list  of  the  input  data  cards  and  their  formats. 

Two  additional  variations  are  possible  within  the  TORPAK  stage  of 
the  analysis;  a  minimum  torque  may  be  applied  to  the  structure  or  an  in 
plane  analysis  may  be  performed  on  the  structure.  In  the  former  case, 
if  the  structure  is  approximately  symmetrical,  the  rotation  of  the  floor 
diaphragms  become  very  small  and  the  centres  of  rotation  are  a  great 
distance  from  the  reference  point.  The  solution  requires  many  iteration 
cycles  and  results  in  floor  deformations  that  are  almost  equal  at  all 
bents.  In  this  case,  it  may  be  required  by  design  codes  that  a  minimum 
torque  be  applied  to  the  structure.  This  is  done  automatically  by  TORPAK 
when  the  centre  of  resistance  of  the  structure  is  within  5  per  cent  of 
the  overall  width  of  the  structure  from  the  geometric  centre.  TORPAK 
applies  a  minimum  torque  equal  to  the  product  of  the  applied  loads  and 
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an  eccentricity  equal  to  10  per  cent  of  the  overall  width  of  the  structure. 
The  second  variation  is  initiated  by  the  input  data.  If  requested, 
TORPAK  performs  a  symmetrical  analysis  on  the  structure,  in  which  all 
bents  in  any  one  direction  translate  equal  amounts.  This  is  achieved 
by  summing  the  lateral  stiffness  matrices  of  all  bents  in  one  direction. 

The  total  lateral  stiffness  matrix  is  then  inverted  and  multiplied  by  the 
applied  lateral  loads  (in  that  direction)  to  yield  the  story  sway  vectors 
(which  will  be  the  same  for  all  bents  in  that  direction).  The  subroutine 
that  carries  out  the  symmetrical  analysis  is  called  SYMPAK  and  after  re¬ 
turning  control  back  to  COMPAK  all  the  torsional  analysis  subroutines  are 
bypassed  except  SWAY  which  calls  up  CALPAK  as  before. 
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FIGURE  4.1  SECTION  OF  TYPICAL  BENT 


FIGURE  4,2  DEFORMED  SECTION  OF  TYPICAL  BENT 
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FIGURE  4s 3  CO-ORDINATE  AXES  AND  BUILDING  ORIENTATION 
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FIGURE  4.4  MAIN  PROGRAM 
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IfiURE  4.4  MAIN  PROGRAM  (CONTINUED) 
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FIGURE  4.5  SUBROUTINE  TORPAK 
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FIGURE  4.5  SUBROUTINE  TORPAK  (CONTINUED) 
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FIGURE  4.5  SUBROUTINE  TORPAK  (CONTINUED) 
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FIGURE  4,5  SUBROUTINE  TORPAK  (CONTINUED) 
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FIGURE  4,5  SUBROUTINE  TORPAK  (CONTINUED) 
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FIGURE  4.6  SUBROUTINE  CALPAK 
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CHAPTER  V 


RESULTS  AND  DISCUSSION 


As  a  check  on  the  validity  of  the  assumptions  used  in  the  analysis, 
several  structures  were  analyzed  and  the  results  compared  with  those  ob¬ 
tained  by  other  methods „  Although  a  method  does  not  exist  for  performing 
a  general  analysis  of  a  three  dimensional  frame  shear-wall  structure,  it 
was  possible  to  check  the  various  sections  of  the  analysis  indi vi dual ly. 

The  calculation  of  the  stiffness  matrix  for  an  uncoupled  frame  was 
checked  by  analyzing  the  twenty  story  structure  shown  in  elevation  in 
FIGURE  5 . 1 A .  A  structure  containing  this  bent  was  used  by  Weaver  and 
Nelson  (5)  to  illustrate  the  rigorous  analysis  of  three-dimensional  framed 
structures,  The  members  were  of  structural  steel  having  a  modulus  of 
elasticity  of  30,000  k,s.i,  and  were  assumed  to  be  coated  with  one  inch 
of  concrete  fireproofing.  The  moments  of  inertia  were  computed  on  the 
basis  of  the  transformed  section  using  a  modular  ratio  of  nine.  The  lat¬ 
eral  load  on  the  structure  consists  of  a  wind  pressure  of  20  p,Sof.,  con¬ 
centrated  at  each  floor  level  and  applied  to  the  south  face  as  shown  on 
the  plan  of  the  structure,  FIGURE  5JB0 

The  first  stage  of  the  comparison  consisted  of  performing  an  in¬ 
plane  analysis  using  the  present  method.  Then,  the  one-bay  bent,  XI, 
was  isolated  as  shown  in  FIGURE  5,2  and  subjected  to  the  lateral  loads 
resisted  by  the  bent  in  the  in-plane  analysis  of  the  entire  structure. 

This  analysis  was  performed  using  STRUDL  (15),  The  column  moments,  joint 
rotations  and  floor  deflections  obtained  by  the  two  analyses  show  excellent 


39 


o  :  :  b:  "i  r>3 

v 


... 

* 

■ 


. 


40 


agreement  and  are  tabulated  in  TABLE  5.1.  The  discrepancies  occur  only 
in  the  third  decimal  place.  The  bents  used  for  the  comparison  were  chosen 
to  eliminate  errors  due  to  the  lumping  technique  employed  in  the  analysis. 

The  calculation  of  the  stiffness  matrix  for  a  coupled  frame  shear- 
wall  bent  was  checked  by  analyzing  the  ten-story  structure  shown  in 
FIGURE  5.3.  The  moment  of  inertia  of  the  shear-wall  was  taken  as  443,000 

4 

inches  and  that  of  the  frame  members  was  based  on  the  CISC  handbook  prop¬ 
erties  for  the  sections  used  (16).  The  lateral  load  consists  of  a  wind 
pressure  of  16  p.s.f.  on  the  long  side  of  the  structure.  This  was  applied 
as  an  equivalent  concentrated  load  at  each  floor  level.  The  columns  were 
of  constant  cross  section  for  the  bottom  five  stories  then  changed  at  the 
sixth  story  and  were  again  constant  for  the  top  five  stories.  All  column 
sections  in  a  given  story  are  the  same.  The  beam  sections  were  the  same 
for  all  floors  in  a  particular  bent.  The  shear  wall  dimensions  were 
8'-0"  x  0‘-6"  and  its  modulus  of  elasticity  was  3,160  k.s.i.  The  modulus 
of  elasticity  of  the  steel  was  taken  as  29,000  k.s.i. 

An  in-plane  analysis  was  performed  on  this  structure  using  the  pre¬ 
sent  method.  Bent  XI  was  then  isolated  as  shown  in  FIGURE  5.4  and  sub¬ 
jected  to  the  lateral  loads  resisted  by  this  bent  in  the  in-plane  analy¬ 
sis  of  the  entire  structure.  This  planar  bent  was  then  analyzed  using  an 
approximate  procedure,  developed  by  Guhamajumdar  et  al  (1),  for  two-dim¬ 
ensional  frame  shear-wall  structures.  The  column  shears,  wall  shears  and 
floor  deflections  obtained  by  the  two  methods  are  presented  in  TABLE  5.2. 
The  agreement  was  again  excellent.  If  the  bents,  in  the  above  two  com¬ 
parisons,  had  consisted  of  more  than  one  bay,  the  lumping  technique  would 
introduce  errors  in  the  stiffness  matrix  calculation. 
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To  provide  a  check  on  the  torsional  portion  of  the  analysis  the 
example  structure  in  FIGURE  5.1  was  again  used.  A  torsional  analysis 
was  performed  using  the  present  method  and  the  results  were  compared 
with  a  "Space  Frame"  analysis  performed  using  STRUDL.  For  the  latter 
analysis  the  cross-sectional  areas  of  all  members  were  set  at  an  artifi¬ 
cially  large  value  (99999.0  square  inches)  so  as  to  suppress  the  effect 
of  axial  deformations.  Pin-ended  diagonal  members  were  introduced  be¬ 
tween  adjacent  bents,  in  a  horizontal  plane,  to  simulate  the  action  of 
floor  diaphragms  having  large  rigidities  in  their  own  planes.  The  floor 
deflections  of  bents  XI  and  X4,  and  the  floor  diaphragm  rotations,  as 
computed  by  the  two  analyses,  have  been  presented  in  TABLE  5.3.  The 
values  in  all  cases  agree  remarkably  well .  This  structure  contains  sev¬ 
eral  multi -bay  bents  which  will  result  in  errors  in  the  stiffness  matrix 
calculation,  for  the  present  analysis.  The  maximum  discrepancy,  however, 
is  less  than  4  per  cent  which  is  acceptable  for  a  structure  of  this  size. 
The  maximum  floor  deflections  and  floor  diaphragm  rotations  obtained 
agree  very  well  with  results  quoted  by  Weaver  and  Nelson  (5). 

Several  structures  have  been  analyzed  by  the  present  method  to  de¬ 
termine  the  effect  of  varying  degrees  of  asymmetry  on  the  structure. 

The  results  of  the  analysis  of  one  such  structure  will  be  presented  in 
detai 1 c 

The  ten  story  structure  in  FIGURE  5.3  was  analyzed  for  four  different 
positions  of  bent  X3.  For  simplicity,  bent  X3  is  considered  to  have  zero 
stiffness  in  the  X  direction.  The  stiffnesses  of  bents  Y1  and  Y2  are 
computed  by  ignoring  the  contribution  of  the  columns  of  bent  X3  and 
thus  have  the  same  properties  for  all  positions  of  X3.  The  shear-wall  has 
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a  constant  section  for  the  full  height  of  the  structure.  Since  the  col¬ 
umn  sizes  change  at  the  sixth  story,  the  ratio  of  rigidity  of  shear-wall 
to  rigidity  of  all  the  columns  is  83  for  the  first  five  stories  and  172 
for  the  top  five  stories „ 

The  position  of  bent  X3  is  defined  by  the  dimension  X,  and  analyses 
were  performed  for  values  of;  0,  90,  180,  270  and  360  inches,  TABLE  5,4 
contains  the  story  shears  developed  by  the  columns  in  bent  XI  and  the 
shear-wall  in  bent  X3;  for  the  various  values  of  X„  In  TABLE  5,5  the 
story  shears,  V^,  for  X  =  180  and  360  inches  have  been  expressed  as  ratios 
of  the  corresponding .shears ,  V^y^,  obtained  from  the  in-plane  analysis 
(X  =  0).  To  show  the  effect  of  the  increasing  asymmetry,  FIGURE  5,5  plots 
the  deflection  at  the  top  floor  of  each  bent  for  values  of  X  from  0  to 
360  inches. 

The  change  in  the  interaction  between  shear-wall  and  frame  elements, 
for  increasing  asymmetry  of  the  shear-wall  bent,  is  clearly  shown  in 
TABLES  5,4  and  5,5,  For  a  symmetric  structure,  the  interaction  between 
the  can ti lever- type  of  deformation  of  the  wall  and  the  shear- type  defor¬ 
mation  of  the  frame,  results  in  the  shear-wall  carrying  almost  all  of  the 
lateral  load  at  the  base  of  the  structure.  Due  to  the  large  deflections 
of  the  wall  in  the  upper  stories  of  the  structure  the  frame  must  act  as 
a  restraining  member.  Consequently  the  frame  carries  a  larger  part  of 
the  overall  lateral  load  in  the  upper  stories  and,  as  shown  in  TABLE  5,4, 
may  actually  carry  more  than  the  externally  applied  lateral  load,  due 
to  the  negative  load  imposed  on  the  frame  by  the  wall. 

In  TABLE  5,4  as  the  value  of  X  increases,  more  and  more  of  the  load 
formerly  carried  by  the  wall  has  to  be  resisted  by  the  frame  of  bent  XI „ 
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As  the  shear  ratios  ^ ^x/ ^SYM ^  s^ow>  change  may  be  very  large.  The 
change  is  greatest  for  the  lower  stories  of  the  structure  and  decreases 
towards  the  top.  The  discontinuity  at  the  sixth  floor  is  due  to  the 
abrupt  change  in  the  column  moment  of  inertia  at  this  level.  It  is  a 
general  observation  that  the  effect  of  any  abrupt  change  in  member  prop¬ 
erties  serves  to  attract  load  to  the  immediate  area, 

FIGURE  5,5  shows  the  same  effect  for  bent  X2,  but  on  a  reduced  scale 
of  magnitude.  As  the  eccentricity  increases,  the  increase  in  deflection 
at  the  top  of  bent  XI  is  of  the  order  of  300  per  cent  while  the  increase 
for  bent  X2  is  less  than  200  per  cent.  The  deflection  characteristic 
of  bent  X4  is  explained  by  the  fact  that  the  torsional  behavior  of  a 
bent  is  influenced  by  its  location  in  relation  to  the  centre  of  rotation, 
as  well  as  its  stiffness.  Up  to  an  X  value  of  approximately  200  inches 
the  centre  of  rotation  of  the  structure  lies  inside  the  structure,  so 
that  as  the  floor  diaphragms  rotate  the  net  deflection  of  bent  X4  is  re¬ 
duced,  Beyond  this  point  the  deflections  tend  to  increase  once  again. 

To  observe  the  effect  of  a  reduction  in  shear  wall  rigidity,  relative 
to  the  sum  of  the  column  rigidities,  in  a  structure,  the  shear-wall  mom¬ 
ent  of  inertia  in  the  structure  of  FIGURE  5,3  was  reduced  from  443,000 
to  80,000  inches^,  The  results  obtained  from  a  series  of  torsional  analy¬ 
ses  for  various  positions  of  the  shear  wall  bent  indicated  the  same  trends 
as  those  of  TABLES  5,4,  5,5  and  FIGURE  5,5,  The  only  difference  was  in 
the  proportion  of  total  lateral  load  carried  by  the  shear-wall  at  any 
one  time. 

These  trends  were  confirmed  by  the  results  of  a  series  of  analyses 
on  a  three  story  structure  with  a  plan  similar  to  that  of  the  ten  story 
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structures  mentioned  above. 

The  symmetrical  structure,  of  FIGURE  5.3,  with  bent  X3  at  X  =  0 
inches,  was  also  analyzed  under  a  torque  equal  to  the  product  of  the 
applied  lateral  load  and  an  eccentricity  of  10  per  cent  of  the  overall 
width  of  the  structure.  As  bents  X2  and  X3  are  at  the  centre  of  resist¬ 
ance  of  the  structure,  the  floor  deflections  and  story  shears  for  these 
two  bents  are  the  same  as  those  obtained  from  an  in-plane  analysis.  For 
this  reason  TABLE  5,6  contains  the  story  shears,  floor  deflections  and 
the  ratios  of  the  story  shears  to  the  shears  obtained  from  an  in-plane 
analysis;  for  bents  XI  and  X4  only. 

The  results  in  TABLE  5,6  illustrate  the  large  change  in  shears  that 
can  occur  in  certain  bents  of  a  symmetric  structure  if  the  structure  is 
subjected  to  a  torque  about  its  longitudinal  axis.  The  increase  in 
shear  in  bent  XI  is  greatest  for  the  bottom  stories  and  gradually  reduces 
with  height.  Because  of  symmetry  the  decrease  in  shear  in  bent  X4  is 
equal  to  the  increase  in  shear  in  bent  XI.  This  decrease  in  shear  was 
sufficient  to  cause  a  sign  reversal  for  the  bottom  story.  Because  bents 
X2  and  X3  are  located  at  the  centre  of  rotation  of  the  structure  there 
is  no  net  change  in  floor  deflections  or  story  shears  for  these  bents. 
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FIGURE  5.1A  EXAMPLE  STRUCTURE  ELEVATION  (TWENTY  STORIES) 
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POINT  OF  LOAD  APPLICATION 


FIGURE  5. IB  EXAMPLE  STRUCTURE  PLAN  (TWENTY  STORIES) 
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FIGURE  5.3  EXAMPLE  STRUCTURE  (TEN  STORIES) 
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50 


FIGURE  5.5 


BENT  DEFLECTIONS  AT  TOP  FLOOR 
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Story 

or 

Floor 

No. 

PRESENT  ANALYSIS 

STRUDL 

Column  Moments 

Joint 

Rotation 

X10~3Rad. 

Floor 

Deflection 

Inches 

Column  Moments 

Joint 

Rotation 

X  10~ 3  Rad 

Floor 

Deflection 

Inches 

Bottom 
Kip  Inch 

Top 

Kip  Inch 

Bottom 

Kip  Inch 

Top 

Kip  Inch 

1 

4172 

-62 

3.4031 

0.324 

4172 

-62 

3.4030 

0.324 

2 

2015 

992 

4.2252 

0.932 

2015 

992 

4.2251 

0.931 

3 

1432  1329 

4.3167 

1.605 

1432 

1329 

4.3165 

1.605 

4 

1147  1313 

4.1691 

2.269 

1147 

1313 

4.1688 

2.269 

5 

1078  1305 

3.9448 

2.910 

1078 

1305 

3.9444 

2.910 

6 

957  1196 

3.7096 

3.512 

957 

1195 

3.7090 

3.512 

7 

932  1167 

3.4522 

4.083 

932 

1166 

3.4513 

4.082 

8 

813 

982 

3.2670 

4.614 

813 

982 

3.2660 

4.613 

9 

891  1011 

3.0563 

5.150 

891 

1011 

3.0551 

5.149 

10 

742 

910 

2.7593 

5.638 

741 

910 

2.7583 

5.637 

11 

672 

808 

2.5206 

6.081 

672 

807 

2.5196 

6.080 

12 

638 

758 

2.2633 

6.497 

637 

757 

2.2623 

6.496 

13 

540 

678 

1.9681 

6.864 

540 

678 

1.9671 

6.863 

14 

450 

539 

1.7793 

7.185 

450 

538 

1.7782 

7.183 

15 

481 

533 

1.5548 

7.531 

481 

533 

1.5536 

7.529 

16 

358 

434 

1.2262 

7.815 

358 

434 

1.2252 

7.812 

17 

269 

322 

0.9965 

8.037 

269 

322 

0.9958 

8.034 

18 

250 

273 

0.7146 

8.308 

249 

273 

0.7140 

8.306 

19 

136 

167 

0.3551 

8.471 

136 

167 

0.3548 

8,469 

20 

37 

58 

0.1016 

8.531 

37 

58 

0.1015 

8.528 

TABLE  5.1  STIFFNESS  MATRIX  CHECK  (UNCOUPLED  BENT) 
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. 


Story 

or 

Floor 

No. 

PRESENT  ANALYSIS 

GUHAMAJUMDAR 

Column 

Shear 

Kips 

Wall 

Shear 

Kips 

Floor 

Deflection 

Inches 

Column 

Shear 

Kips 

Wall 

Shear 

Kips 

Floor 

Deflection 

Inches 

1 

7.63 

85.21 

0.110 

7.62 

85.22 

0.110 

2 

13.27 

59.01 

0.347 

13.27 

59.01 

0.347 

3 

14.60 

44.33 

0.626 

14.60 

44.32 

0.626 

4 

13.88 

35.06 

0.903 

13.88 

35.06 

0.903 

3 

14.11 

23.49 

1.160 

14.11 

23.49 

1.160 

6 

8.72 

27.45 

1.389 

8.72 

27.45 

1.389 

7 

8.45 

16.85 

1.579 

8.45 

16.85 

1.579 

8 

6.35 

11.17 

1.727 

6.34 

11.18 

1.726 

9 

4.44 

4.93 

1.836 

4.43 

4.95 

1.834 

10 

4.25 

-5.60 

1.916 

4.25 

-5.61 

1.913 

TABLE  5 , 2 


STIFFNESS  MATRIX  CHECK  (COUPLED  BENT) 
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PRESENT  ANALYSIS 

STRUDL 

Floor 

Number 

Floor  Deflection 
Inches 

Floor 

Rotation 

ofl 

Floor  Deflection 
Inches 

Floor 

Rotation 

Rad.  XIO"2 

Bent  XI 

Bent  X4 

Rad.  X10  | 

Bent  XI 

|  Bent  X4 

1 

0.421 

0.252 

0.020 

0.422 

0.256 

0.019 

2 

1.226 

0.712 

0.059 

1.229 

0.726 

0.058 

3 

2.129 

1.216 

0.106 

2.136 

1.241 

0.103 

4 

3.023 

1.708 

0.152 

3.035 

1.746 

0.149 

5 

3.891 

2.182 

0.198 

3.905 

2.232 

0.194 

6 

4.708 

2.625 

0.241 

4.725 

2.689 

0.236 

7 

5.483 

3.046 

0.282 

5.502 

3.123 

0.275 

8 

6.203 

3.438 

0,320 

6.226 

3.527 

0.312 

9 

6.925 

3.835 

0.358 

6.953 

3.939 

0.348 

10 

7.584 

4.198 

0.392 

7.616 

4.317 

0.382 

11 

8.182 

4.526 

0.423 

8.217 

4.659 

0.412 

12 

8.742 

4.837 

0.452 

8.780 

4.983 

0.439 

13 

9.236 

5.110 

0.478 

9.277 

5.270 

0.464 

14 

9.671 

5.347 

0.500 

9.713 

5.521 

0.485 

13 

10.130 

5.611 

0.523 

10. 175 

5.799 

0.506 

16 

10.504 

5.826 

0.541 

10.553 

6.027 

0.524 

17 

10.796 

5.997 

0.555 

10.851 

6.207 

0.537 

18 

11:148 

6.209 

0.572 

11.202 

6.433 

0.552 

19 

11.357 

6.338 

0.581 

11.414 

6.571 

0.561 

20 

11.443 

6.379 

0.586 

11.518 

6.618 

0.567 

TABLE  5.3  TORSIONAL  ANALYSIS  CHECK  (FRAMED  STRUCTURE) 
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Floor 

Number 

EFFECT  OF  INCREASING  ASYMMETRY 

Story  Shear  in  Bent  XI 

Kips 

Story  Shear  in  Shear  Wall  Kips 

X 

0 

90 

180 

270 

360 

0 

90 

180 

270 

360 

1 

5.53 

16.68 

27.35 

35.31 

40.22 

85.21 

79.88 

67.01 

52.94 

40.94 

2 

8.55 

16.31 

24.12 

30.35 

34.54 

59.01 

56.09 

48.74 

39.98 

31.93 

3 

9.16 

15.43 

21.62 

26.76 

30.28 

44.33 

42.40 

37.29 

31.08 

25.17 

4 

8.65 

13.80 

18.97 

23.24 

26.21 

35.06 

33.59 

29.74 

24.97 

20.36 

5 

8.59 

12.30 

16.05 

19.35 

21.84 

23.49 

22.89 

20.98 

18.23 

15.27 

6 

5.22 

9.76 

14.04 

17.34 

19.32 

27.45 

25.83 

21.84 

17.45 

13.67 

7 

5.01 

8.05 

10.99 

13.28 

14.87 

16.85 

16.02 

13.92 

11.42 

9.13 

8 

3.76 

5.86 

7.94 

9.56 

10.60 

11.17 

10.60 

9.25 

7.62 

6.11 

9 

2.63 

3.78 

4.84 

5.74 

6.33 

4.93 

4.76 

4.21 

3.55 

2.91 

10 

2.37 

2.12 

1.92 

1.89 

1.97 

-5.61 

-5.10 

-3.96 

-2.80 

-1.94 

TABLE  5.4  SHEAR  TABULATION  FOR  BENT  XI  AND  SHEAR-WALL 
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Floor 

Number 

SHEAR  RATIOS  VX/VsyM 

Bent  XI 

Shear  Wall 

X 

180 

360 

180 

360 

1 

4.95 

7.28 

0.79 

0.48 

2 

2.82 

4.04 

0.82 

0.54 

3 

2.36 

3.31 

0.84 

0.57 

4 

2.19 

3.03 

0.85 

0.58 

5 

1.87 

2.55 

0.89 

0.65 

6 

2.69 

3.70 

0.80 

0.50 

7 

2.20 

2.97 

0.83 

0.54 

8 

2.11 

2.82 

0.83 

0.55 

9 

1.84 

2.41 

0.85 

0.59 

10 

0.81 

0.83 

0.71 

0.35 

TABLE  5.5  SHEAR  RATIOS  BENT  XI  AND  SHEAR-WALL 
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BENT  X 1 

BENT  X4 

Floor 

Shear 

Deflection 

Ratio 

Shear 

Deflection 

Ratio 

Min.  Torque 

of 

Min.  Torque 

of 

Number 

Kips 

Inches 

Shears 

Kips 

Inches 

Shears 

1 

14.38 

0.241 

2.60 

-3.32 

-0.021 

-0.60 

2 

15.60 

0.683 

1.83 

1.49 

0.011 

0.17 

3 

15.34 

1.158 

1.68 

2.97 

0.093 

0.32 

4 

13.92 

1.609 

1.61 

3.35 

0.198 

0.39 

5 

12.88 

2.015 

1.50 

4.29 

0.304 

0.50 

6 

9.62 

2.408 

1.84 

0.83 

0.369 

0.16 

7 

8.34 

2.729 

1.67 

1.67 

0.429 

0.33 

8 

6.14 

2.970 

1.63 

1.37 

0.483 

0.36 

9 

4.09 

3.135 

1.56 

1.18 

0.535 

0.45 

10 

2.83 

3.239 

1.19 

1.92 

0.592 

0.81 

■i 

TABLE  5.6  COMPARISON  OF  IN-PLANE  AND 


MINIMUM  TORQUE  ANALYSES 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

A  method  of  analysis  has  been  presented  which  determines  the  distri¬ 
bution  of  resisting  shears  among  the  various  bents  in  asymmetric  three- 
dimensional  structures.  The  analysis  neglects  second  order  effects  and 
assumes  that  the  material  is  ideally  elastic.  The  bents  of  the  structure 
may  contain  either  coupled  or  independent  shear-wall  and  frame  elements. 

The  actual  structure  has  been  simplified  to  a  system  of  two-dimen¬ 
sional  bents.  A  lumping  technique  has  been  used  to  reduce  the  individual 
bents  to  single  equivalent  frame  and  shear-wall  elements.  Each  floor 
diaphragm  is  free  to  rotate  about  a  vertical  axis,  and  to  translate  in 
the  X  and  Y  directions. 

The  computer  program  developed  for  the  analysis,  performs  a  torsional 
analysis  if  the  structure  is  asymmetric  but  when  the  degree  of  asymmetry 
is  minor  a  specified  torque  is  applied  to  the  structure.  The  program 
also  performs  an  in-plane  analysis  if  required.  A  list  of  the  subrou¬ 
tines  used  in  the  computer  program,  and  the  required  input  data  is  pre¬ 
sented  in  APPENDIX  B.  A  print-out  of  the  computer  program  is  included 
in  APPENDIX  C. 

The  various  steps  in  the  analysis  were  verified  by  performing  analy¬ 
ses  on  segments  of  the  structure  for  which  other  methods  were  available. 
The  present  method  was  then  used  to  analyze  several  structures  containing 
coupled  frame  and  shear-wall  bents.  The  results  of  these  analyses  clearly 
show  the  important  effect  that  the  rotation  of  the  structure  has  on  the 
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shear  distribution,,  With  increasing  asymmetry  of  the  structure  the  in¬ 
teraction  between  shear-wall  and  frame  elements  is  changed  dramatically. 
Even  for  small  eccentricities  the  increased  shears  imposed  on  certain 
bents  near  the  perimeter  of  the  structure  are  significant  and  must  be 
considered. 

For  the  series  of  ten  story  structures  used  to  illustrate  the  ef¬ 
fect  of  torsion,  the  number  of  iteration  cycles  required  for  convergence 
varied  between  ten  and  twenty.  The  corresponding  computer  time  amounted 
to  approximately  0,8  minutes.  For  the  twenty  story  structure  convergence 
was  considerably  slower.  This  was  because  the  structure  does  not  have 
an  axis  of  symmetry,  so  that  both  the  X  and  Y  co-ordinates  of  the  centre 
of  rotation  must  be  determined  by  iteration.  Each  iterative  cycle  re¬ 
quires  the  inversion  of  more  than  twenty  matrices.  To  alleviate  this 
problem  it  is  suggested  that  the  program  be  modified  to  a  more  rigorous 
solution,  in  which  the  three  equilibrium  equations,  for  each  floor,  will 
be  solved  at  the  one  time.  This  approach  will  lead  to  reduced  computer 
time,  at  the  expense  of  increased  computer  storage  space. 

The  next  stage  to  be  considered  is  the  inelastic  response  of  the 
structure.  The  greatest  increase  in  shear  due  to  torsional  action  occurs 
in  the  bents  located  furthest  from  the  centre  of  rotation.  The  formation 
of  plastic  hinges  in  these  bents  reduces  the  bent  stiffness  and  would 
result  in  a  shift  of  the  centre  of  resistance  as  the  structure  is  loaded 
into  the  inelastic  range.  In  future  extensions  of  the  program  the  struc¬ 
ture  will  be  analyzed  as  described  until  a  plastic  hinge  forms.  The  stiff¬ 
ness  matrix  for  the  particular  bent  involved  will  then  be  modified  by 
the  insertion  of  a  plastic  hinge  in  the  member.  The  structure  will  be 
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analyzed  with  the  modified  stiffness  matrix  until  a  second  hinge  forms. 

The  process  is  then  repeated.  At  each  load  increment,  iteration  will  be 
used  to  obtain  a  balance  between  the  applied  loads  and  the  resisting  floor 
forces  using  the  actual  stiffness  matrix  for  each  bent. 

Additional  effects  that  must  be  investigated  are  primarily  due  to 
the  axial  loads;  the  influence  of  axial  shortening,  the  P-A  effect  caused 
by  the  vertical  loads  and  the  change  in  carry-over  and  stiffness  coeffi¬ 
cients  due  to  the  presence  of  axial  load. 
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NOMENCLATURE 


u 

Displacement  in  the  X  di recti  on . 

V 

Displacement  in  the  Y  direction. 

X 

V* 

-< 

General  co-ordi nates e 

<c 

>- 

0% 

<c 

X 

Co-ordinates  of  general  point  A. 

X  ,  Y 

0 9  0 

Co-ordinates  of  the  centre  of  rotation. 

X 

o 

ss 

-< 

o 

Co-ordinates  of  the  centre  of  resistance,, 

NX,  NY 

Number  of  bents. 

F 

Applied  lateral  load. 

V 

Resisting  floor  force. 

A 

Story  sway  deformation. 

S 

Displacement  from  the  undeformed  position. 

6 

Floor  diaphragm  rotation. 

D 

Two-diagonal  matrix. 

T 

Torsional  stiffness  matrix. 

C 

Transpose  of  influence  coefficient  matrix. 

AX 

Centre  of  rotation  correction  in  the  X  direction 

AY 

Centre  of  rotation  correction  in  the  Y  direction 

E 

Modulus  of  elasticity  of  frame. 

EW 

Modulus  of  elasticity  of  wall. 

h 

Story  height. 

Cl 

Moment  of  inertia  of  equivalent  column. 

WI 

Moment  of  inertia  of  equivalent  wall. 

0C 

Beam  to  column  joint  rotation. 
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ew 

MC 

MW 

FBI 

FBL 

WBI 

WBL 


Beam  to  wall  joint  rotation,, 

Resultant  moment  at  beam  to  column  joint. 
Resultant  moment  at  beam  to  wall  joint. 
Moment  of  inertia  of  equivalent  frame  beam, 
Length  of  equivalent  frame  beam. 

Moment  of  inertia  of  equivalent  wall  beam. 
Length  of  equivalent  wall  beam. 


V 

- 
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NOMENCLATURE  FOR  FORTRAN  IV  PROGRAM 


In  the  following  nomenclature,  X  and  Y  in  any  of  the  terms  refer 
to  the  direction  of  the  perpendicular  to  the  plane  of  the  bent, 

KJ  Matrix  relating  column  joint  rotation  to  the  story  sway  for 

an  uncoupled  frame,, 

KP  Matrix  relating  column  joint  rotation  to  the  story  sway  for 

coupled  frame, 

KS  Matrix  relating  wall  joint  rotation  to  the  story  sway  for  a 

coupled  frame  and  shear-wall  bent, 

KX  Bent  stiffness  matrix  in  terms  of  story  sway, 

KY  Bent  stiffness  matrix  in  terms  of  story  sway, 

KV  St,  Venant  torsional  stiffness  matrix, 

K  Bent  stiffness  matrix  term  used  in  subroutines  MACAL1  and 

MACAL2, 

NS  Number  of  stories, 

NX,  NY  Number  of  bents, 

ICX  Dummy  term  representing  the  number  of  iteration  cycles, 

ICY  Dummy  term  representing  the  number  of  iteration  cycles, 

XO  Distance  between  reference  point  and  centre  of  rotation, 

YO  Distance  between  reference  point  and  centre  of  rotation, 

KXROTC  Matrix  relating  column  joint  rotation  to  the  story  sway, 
KXROTW  Matrix  relating  wall  joint  rotation  to  the  story  sway, 
KYROTC  Matrix  relating  column  joint  rotation  to  the  story  sway, 
KYROTW  Matrix  relating  wall  joint  rotation  to  the  story  sway. 
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H 

Story  height* 

Cl 

Equivalent  column  moment  of  inertia* 

WI 

Equivalent  wall  moment  of  inertia* 

WBL 

Wall  beam  length* 

U 

Wall  width* 

FBL 

Frame  beam  length* 

FBI 

Frame  beam  moment  of  inertia* 

WBI 

Wall  beam  moment  of  inertia* 

E 

Modulus  of  elasticity  of  frame* 

EW 

Modulus  of  elasticity  of  wall. 

EGC 

Rotational  stiffness  of  foundation  under  column* 

EGW 

Rotational  stiffness  of  foundation  under  wall* 

FX 

Applied  lateral  load  in  the  Y  direction* 

FY 

Applied  lateral  load  in  the  X  direction* 

XX 

Distance  between  reference  point  and  particular 

bent* 

YY 

Distance  between  reference  point  and  particular 

bent* 

cx 

Distance  between  reference  point  and  centre  of  resistance  (El) 

CY 

Distance  between  reference  point  and  centre  of  resistance  (El) 

RX 

Distance  between  a  particular  bent  and  the  centre  of  rotation* 

DX 

Relative  distance  between  top  and  bottom  of  each 

story  and 

their  respective  centres  of  rotation* 

KDX 

Stiffness  matrix  in  terms  of  unit  rotation  about 

the  centres 

of  rotation* 

TTX 

Torsional  stiffness  contribution,  computed  about  the  reference 

point* 

' 
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ROT 

VX 

TVX 

RY 

DY 

KDY 

TTY 

VY 

TVY 

MT 

TVXA 

TVYA 

XH 

YH 

INDICX 

INDICY 

HX 

CIX 

WIX 

WBLX 


Floor  diaphragm  rotation  vector. 

Floor  forces  developed  by  a  particular  bent. 

Sums  of  the  floor  forces. 

Distance  between  a  particular  bent  and  the  centres  of  rotation. 
Relative  distance  between  top  and  bottom  of  each  story  and  their 
respective  centres  of  rotation. 

Stiffness  matrix  in  terms  of  unit  rotation  about  the  centres 
of  rotation. 

Torsional  stiffness  contribution,  computed  about  the  reference 
point. 

Floor  forces  developed  by  a  particular  bent. 

Sum  of  the  floor  forces. 

Torsional  moment,  about  the  reference  point,  due  to  the 
applied  lateral  loads. 

Sum  of  the  floor  forces  for  the  initially  assumed  position  of 
the  centres  of  rotation. 

Sum  of  the  floor  forces  for  the  initially  assumed  position  of 
the  centres  of  rotation. 

Distance  between  reference  point  and  FX. 

Distance  between  reference  point  and  FY, 

Dummy  term  indicating  whether  frame  is  coupled  or  not. 

Dummy  term  indicating  whether  frame  is  coupled  or  not. 

Story  height. 

Equivalent  column  moment  of  inertia. 

Equivalent  wall  moment  of  inertia. 

Wall  beam  length. 
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wx 

Wall 

width. 

FBLX 

Frame  beam  length. 

WBIX 

Wall 

beam  moment  of  inertia 

FBIX 

Frame  beam  moment  of  inerti 

HY 

Story  height. 

CIY 

Equivalent  column  moment  of 

WIY 

Equivalent  wall  moment  of  ii 

WBLY 

Wall 

beam  length. 

WY 

Wall 

wi dth . 

FBLY 

Frame  beam  length. 

WBIY 

Wall 

beam  moment  of  inertia 

FBIY 

Frame  beam  moment  of  inerti, 

KA 

Bent 

stiffness  submatrix. 

KB 

Bent 

stiffness  submatrix. 

KC 

Bent 

stiffness  submatrix. 

KD 

Bent 

stiffness  submatrix. 

KE 

Bent 

stiffness  submatrix. 

KF 

Bent 

stiffness  submatrix. 

KG 

Bent 

stiffness  submatrix. 

KH 

Bent 

stiffness  submatrix. 

KI 

Bent 

stiffness  submatrix. 

KK 

Work 

matrix. 

KL 

Work 

matrix. 

KM 

Work 

matrix. 

KN 

Work 

matrix. 

KO 

Work 

matrix. 
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KQ 

KR 

KT 

KW 

TX 

TY 

RT 

SX 

SY 

ROTXC 

ROTXW 

ROTYC 

ROTYW 

TVXC 

TTVXC 

TVYC 

TTVYC 

CMXB 

CMXT 

WMXB 

WMXT 

FBMX 

WBMXC 

WBMXW 


Work  matrix. 

Work  matrix. 

Work  matrix. 

Work  matrix. 

Torsional  stiffness  contribution,  computed  about  the  reference 
point. 

Torsional  stiffness  contribution,  computed  about  the  reference 
point. 

Torsional  stiffness  matrix  about  the  reference  point. 

Story  sway  vector. 

Story  sway  vector. 

Column  joint  rotation  vector. 

Wall  joint  rotation  vector. 

Column  joint  rotation  vector. 

Wall  joint  rotation  vector. 

Influence  coefficient  matrix. 

Transpose  of  TVXC. 

Influence  coefficient  matrix. 

Transpose  of  TVYC. 

Moment  at  bottom  of  column. 

Moment  at  top  of  column. 

Moment  at  bottom  of  wall. 

Moment  at  top  of  wall . 

Moment  at  the  ends  of  frame  beam. 

Moment  at  column  end  of  wall  beam. 

Moment  at  wall  end  of  wall  beam. 
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CMYB 

CMYT 

WMYB 

WMYT 

FBMY 

WBMYC 

WBMYW 

CSX 

WSX 

FBSX 

WBSX 

CSY 

WSY 

FBSY 

WBSY 

CFX 

WFX 

FBFX 

WBFX 

CFY 

WFY 

FBFY 

WBFY 

DEX 

DEY 

INSYM 


Moment  at  bottom  of  column. 

Moment  at  top  of  column. 

Moment  at  bottom  of  wall. 

Moment  at  top  of  wall . 

Moment  at  the  end  of  frame  beam. 

Moment  at  column  end  of  wall  beam. 

Moment  at  wall  end  of  wall  beam. 

Column  shear. 

Wall  shear. 

Frame  beam  shear. 

Wall  beam  shear. 

Column  shear. 

Wall  shear. 

Frame  beam  shear. 

Wall  beam  shear. 

Column  axial  load. 

Wall  axial  load. 

Frame  beam  axial  load. 

Wall  beam  axial  load. 

Column  axial  load. 

Wall  axial  load. 

Frame  beam  axial  load. 

Wall  beam  axial  load. 

Floor  deflection. 

Floor  deflection. 

Dummy  term  indicating  whether  a  rotational  analysis  is  to 
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ITERX 

ITERY 

IX 

IY 


be  performed. 

Dummy  term  that  limits  the  number  of  iteration  cycles 
Dummy  term  that  limits  the  number  of  iteration  cycles 
Dummy  term  indicating  the  number  of  resets. 

Dummy  term  indicating  the  number  of  resets. 
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APPENDIX  A 


DERIVATION  OF  BENT  STIFFNESS  MATRIX 

For  the  section  of  the  bent  shown  in  FIGURE  4.1,  three  equilibrium 
equation  may  be  written  in  terms  of  slopes  and  deflections. 


V ( i )  =  A(i)  x 


-12  x  E  x  Cl (i )  12  x  EW  x  Wl(i ) 


(h(i ) ) 


(h(i ) ) 


A(i+1)  x 


12  x  E  x  Cl (i  +  1 )  +  12  x  EW  x  WI (i  +  1 ) 


0C(i-l)  x 


(h(i+l))3 

(h(i+l)}3 

6  x  E  x  C I  ( i ) 

+  0C(i )  x 

6  x  E  x  CI(i)  6 

(h(i))2 

( h  ( i ) ) 2 

6  x  E  x  Cl 


(i  +  1) 


(h(i+l ) )' 


+  ec(i+i)  x 


-6  x  E  x  CI(.+1) 

-1-  Ql»l  .  .  V 

6  x  EW  x  WI (i ) 

( h { i + 1 ) )2 

0W(i-l) 

(h(i))2 

0W(i) 


6  X  EW  X  WI(i)  _  6  X  EW  X  WI (i+1 ) 

(h(i))2  (h(i+i))2 


ew( i+i )  x 


-6  x  EW  x  WI (i+1 ) 
(h(i+l))2 


(A.  1 ) 


MC(i) 


A  ( i )  x 


ec(i-l) 


-6  x  E  x  Cl (i ) 

+  A  ( i + 1 )  x 

-6  x  E  x  Cl (i+1 ) 

(h(i))2 

(h(i+l))Z 

2  X  hU)Ci('^'  +  6C(1)  x 

4  x  E  x  CI(i) 

— urn — 

+ 
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4  x  E  x  Cl  (i+1 )  x  J2_x  E  x  FBI  (i )  ,  4  x  E  x  WBI(i) 

h(i  +  D  FBL(i )  WBL(i) 


0C(i+l)  x 


2  x  E  x  Cl (i+1 ) 

hXT+T) 


+ 


ew(i)  x 


2  x  E  x  WBI(i) 

WBL(i) 


+  3  x  w(i)  x  E  x  WBI(i) 

(WBL(i))2 


(A. 2) 


MW(i )  =  A(i )  x 

-6  x  EW  x  WI (i ) 

+  A( i +1 )  x 

-6  x  EW  x  WI (i+1 ) 

o 

(h(i))2 

(h  (i+1  )r 

0  C  ( i )  x 


3  x  E  x  WB 1(1)  x  W(i) 
(WBL(i))2 


,  2  x  E  x  WBI(i) 
WBL(i) 


+ 


ew(i-i)  x 


2  X  EW  X  W I  ( i ) 

— m — 


+  ew(i )  x 


4  x  EW  x  W I  ( i ) 

- FTTH - 


4  x  EW  x  WI (i+1 ) 

hH+T) 


+  6  x  E  x  WBI(i)  x  W(i )  +  3  x  E  x  WBI(i)  x  (W(i))2 

(WBL(i))2  (WBL(i))3 


4  x  E  x  WBI(i) 
+  "  WBL(i) 


+  0W( i+1 )  x 


2  x  EW  x  WI (i+1 ) 

- RTTm - 


(A. 3) 


Where  V(i)  is  the  force  exerted  on  the  bent  by  floor  diaphragm  i,  MC(i) 
is  the  resultant  moment  at  the  beam  to  column  joint  of  floor  level  i  and 
MW(i)  is  the  resultant  moment  at  the  shear  wall  to  beam  joint  of  floor 
level  i .  The  remaining  terms  are  defined  in  the  nomenclature.  The 
above  three  equations  may  be  written  for  each  of  the  n  floors  in  the  bent. 
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This  procedure  results  in  3  x  n  equations  for  3  x  n  unknowns  which  may 
be  expressed  in  matrix  form  as  follows. 


|  V 1 

|  KA  | 

|  KB  | 

|  KC  | 

|  A 1 

|M|C 

= 

|  KD  | 

|  KE  | 

|  KF  | 

X 

|  e  | c 

|  M  | w 

|  KG  | 

|  KH  | 

|  KI  | 

|e|w 

where  KA,  KB,  KC,  KD,  KE,  KF,  KG,  KH  and  KI  are  square  submatrices  con¬ 
taining  n  x  n  elements,  most  of  which  are  zero  terms.  The  elements  of 
the  submatrices  consist  of  the  member  stiffness  terms  expressed  in 

J.  L. 

Equations  (A.l)  to  (A. 3)  above.  As  an  example  the  l  row  of  submatrix 

t  h 

KE  expresses  the  resultant  moment  at  the  beam  to  column  joint  on  the  i tr 
floor  due  to  unit  rotation  of  6C(i-l),  0C(i)  and  6C(i+l).  All  columns 
of  row  i  would  contain  zeros  except  for  the  following; 


column  (i-1) 


2  x  E  x  Cl  (i ) 

- hTTJ - 


column  (i) 


4  x  E  x  Cl (i )  ,  4  x  E  x  Cl (i+1 )  .  12  x  E  x  FBI (i ) 

- RTT} -  hTT+T)  FBLTD 


4  X  E  X  WBI(i) 
WBL(i) 


column  (i+1) 


2  x  E  x  Cl (i+1 ) 

— rtt+d — 


The  matrix  Equation  (A. 4)  may  be  expanded  into  the  following; 
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|V|  =  | KA |  x  | A |  +  | KB |  x  |e|c+  | KC |  x  | 0 |w  (A. 5a) 

| M | C  =  | KD |  x  | A |  +  | KE |  x  | e | c  +  |KF|  x  |e|w  (A. 5b) 

|  M  |W  =  |  KC  |  x  j  A |  +  |  KH  |  x  |  0  |  C  +  |  KI  |  x  |  0  |w  (A, 5c) 

In  order  to  satisfy  equilibrium  the  resultant  moments  at  the  beam  to 
column  joints  and  at  the  beam  to  wall  joints  must  be  equal  to  zero.  It 
is  therefore  possible  to  eliminate  the  joint  rotation  vectors  |e|  and 
| 0 |w  from  Equation  (A. 5a).  Transposing  Equations  (A. 5b)  and  (A. 5c)  yields, 

| e | C  =  -  [ KE r 1  X  ( | KD |  X  | A |  +  | KF|  x  1 6 |W)  (A. 6) 

|  e  |w  =  -  I  Kir1  x  ( [  kg  [  x  |  &  [  +  |  kh  |  x  |  e  5 c)  (a. 7) 

and  substituting  for  | 0 |w  in  Equation  (A. 6), 

|  e  I  C  =  [  ( |  KF  |  X  |  Kir1  x  |  KH  | )  -  |  KE  |  ]_1  x  [  |  KD  |  -  |KF| 

x  i ki r 1  x  ikgi]  x  i*i  (a. 6a) 

or  |0|C  =  |K|0C  x  |A|  (A.6b) 

Substituting  for  j 6 1 c  in  Equation  (A. 7), 

| 0 |w  =  -  | KI|-’  x  [|KGI  x  1*1  +  lKHl  x  lKpl  x  lAH  (A. 7a) 

or  | 6 |w  =  |K|0W  x  | A ]  (A, 7b) 


and  substituting  for  1 0 1 c  and  |0|w  in  Equation  (A. 5a)  yields, 
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V 


KA I  x  I A I  +  I  KB 


K 


0C 


KC 


K 


0w 


x  I A |  (A. 5d) 

or  | V |  =  | K | A  x  | a |  (A.5e) 

where  | K |  is  the  bent  stiffness  matrix  in  terms  of  the  story  sway  defor¬ 
mations  . 

When  the  bent  is  composed  of  an  uncoupled  frame  all  terms  containing 
WI(i)  are  eliminated  in  Equation  (A.l);  in  Equation  (A. 2),  all  terms  con¬ 
taining  WBI(i)  are  eliminated  and  Equation  (A.3)  is  eliminated  entirely. 
Equation  (A. 4)  is  then  written  as. 


|V| 

|  KA 1  | 

1  KB  | 

|A| 

|M|C 

|  KD  | 

LU 

X 

|e|c 

(A. 8) 


where  the  submatrices  have  the  same  significance  as  in  Equation  (A. 4) 
except  for  (KA'|  and  |KE'|,  which  have  been  modified  as  described  above. 
Furthermore,  for  an  uncoupled  frame  Equations  (A. 5c),  (A. 7),  (A. 7a)  and 
(A. 7b)  are  eliminated  and  Equation  (A.5d)  now  reads, 

|  v  |  =  1  KA 1  |  x  |  A |  +  |  KB  |  x  |  K |  0 C  x  |  A  |  (A. 9) 

where  IKA'I  and  I K ' | 0C  have  been  modified  due  to  the  absence  of  the  shear 
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APPENDIX  B 


SUBROUTINE  LISTING  AND  DATA  CARDS 


In  the  following  subroutines,  X  and  Y  (where  they  occur)  denote  the 

directions  of  a  line  constructed  perpendicular  to  the  plane  of  a  bent. 

1.  COMPAK:  -  Calls  up  subroutines  to  read  in  the  member  properties  and 

calculate  the  bent  stiffness  matrix. 

2.  READ1 :  -  Reads  in  the  member  properties  for  an  uncoupled  bent. 

3.  READ2:  -  Reads  in  the  member  properties  for  a  coupled  bent. 

4.  MACALl :  -  Calculates  the  bent  stiffness  matrix  using  the  submatrices 
calculated  by  READ! . 

5.  MACAL2:  -  Calculates  the  bent  stiffness  matrix  using  the  submatrices 
calculated  by  READ2. 

6.  TOKYO:  -  Matrix  inversion  subroutine, 

7.  TORPAK:  -  Calls  up  subroutines  to  develop  the  rotational  stiffness 
matrix  for  the  structure,  and  also  tests  for  equilibrium  of  resisting 
floor  forces  and  applied  loads.  If  there  is  an  unbalance,  subroutines 
are  called  which  carry  out  the  correction  technique, 

8.  SYMPAK:  -  When  called  up  by  TORPAK  this  performs  an  in-plane  analy¬ 
sis  on  the  structure, 

9.  WRITE1 :  -  Prints  out  member  properties,  loads,  number  of  stories  and 
bent  location  within  the  structure. 

10.  HXODX:  -  Computes  the  two-diagonal  geometric  matrix,  which  relates 
the  bent  location  to  the  centres  of  rotation  of  successive  floors. 

11.  DXOTTX:  -  Computes  the  torsional  stiffness  matrix  contribution. 
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12.  VXOTVX:  -  After  the  floor  rotations  have  been  calculated  this  sub¬ 
routine  calculates  the  resisting  floor  forces. 

13.  HYODY:  -  Computes  the  two-diagonal  geometric  matrix,  which  relates 
the  bent  location  to  the  centres  of  rotation  of  successive  floors, 

14.  DYOTTY:  -  Computes  the  torsional  stiffness  matrix  contribution, 

15.  VYOTVY:  -  After  the  floor  rotations  have  been  calculated  this  sub¬ 
routine  calculates  the  resisting  floor  forces, 

16.  RTOROT:  -  Adds  up  the  three  contributions  to  the  rotational  stiff¬ 
ness  matrix  then  inverts  the  total  and  calculates  the  rotation  vector. 

17.  SWAY:  -  When  the  story  sway  vectors  have  been  calculated,  this  sub¬ 
routine  calculates  the  joint  rotation  vectors  for  all  bents. 

18.  WRITE3:  -  Prints  out  the  resisting  floor  forces,  location  of  the 
centres  of  rotation,  rotation  vector,  the  number  of  iteration  cycles 
in  the  X  and  Y  directions,  and  whether  there  is  an  unbalance  in  the 
resisting  floor  forces  and  applied  loads  after  the  last  iteration 
cycle. 

19.  ICXO:  -  This  subroutine  corrects  the  positions  of  the  assumed  centres 
of  rotation,  in  the  X  direction. 

20.  ICYO:  -  This  subroutine  corrects  the  positions  of  the  assumed  centres 
of  rotation,  in  the  Y  direction. 

21.  CALPAK:  -  Calls  up  subroutines  to  calculate  the  member  forces, 

22.  CALC1 :  -  Calculates  the  member  forces  and  bent  deformations  for  the 
bents  located  in  the  X  plane. 

CALC2 :  -  Calculates  the  member  forces  and  bent  deformations  for  the 
bents  located  in  the  Y  plane. 


23. 


„ 

f rU  ,t  :  ubfjs  a  se  u*  )9\  yswa  yioy  9  \  .*»' 

6  «*cft  2T'i  99V  no  1 0"<  in  0 c  >rti  aeifiTuofeD  •TiiubY 

'  0  '.go  :  '.r  -  tug  ainrV?  -  ri 5  ;i 


1 

m 


M,.  "  f:  '  a  Hi  enoiii  ’O' 


. 


76 


24.  WRITE4:  -  Prints  out  the  member  forces  and  bent  deformations  for 
all  the  bents  in  the  structure. 

Data  is  transmitted  to  the  program  by  means  of  punched  cards.  The 
contents  of  each  card  is  explained  below  and  the  order  of  cards  must  be 
the  same  as  the  numeric  order  in  which  they  are  listed. 

1.  NS,  NX,  NY;  Format  (3110).  Where  NS  is  the  number  of  stories,  NX 
the  number  of  bents  in  the  X  plane  and  NY  the  number  of  bents  in 
the  Y  plane.  (A  bent  in  the  X  plane  has  a  perpendicular  which  is 
parallel  to  the  X  axis. ) 

2.  E,  EW;  Format  (2F10.2).  Where  E  is  the  modulus  of  elasticity  for 
the  frame  and  EW  is  the  modulus  of  elasticity  for  the  shear  wall. 

The  units  of  both  constants  are  kips  per  square  inch. 

3.  INDICX  or  INDICY;  Format  (110).  Where  INDIC  is  the  indicator  that 

tells  whether  the  bent  consists  of  a  coupled  frame  and  shear  wall 
or  an  uncoupled  frame.  The  numeral,  1,  in  column  10  indicates  a 
coupled  frame  shear  wall  bent  while  the  numeral,  0,  in  column  10 
indicates  an  uncoupled  frame.  This  particular  card  must  be  repeated 
for  each  bent  and  precedes  each  bent  member  property  card  group  as 
defined  below. 

4.  H,  Cl,  WI,  W,  WBL,  FBL,  WBI ,  FBI;  Format  (F8.2,  2F12.2,  F8.2,  4F10.2), 

Where  for  each  story  or  floor  level  H  is  the  story  height,  Cl  is  the 

lumped  column  moment  of  inertia,  WI  is  the  lumped  shear  wall  moment 
of  inertia,  W  is  the  shear  wall  width  (in  the  plane  of  the  bent), 

WBL  is  the  wall  beam  length,  FBL  is  the  frame  beam  length,  WBI  is  the 
wall  beam  moment  of  inertia  and  FBI  is  the  frame  beam  moment  of  in¬ 
ertia.  The  units  for  all  of  the  above  are  inches.  For  each  bent 
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there  must  be  firstly  an  INDIC  card  followed  by  a  member  property 
card  for  each  floor  in  the  bent.  Then  an  INDIC  card  for  the  next 
bent  and  so  on.  For  the  case  where  the  bent  consists  of  an  uncoup¬ 
led  frame;  WI  and  WBI  are  left  blank  but  W  and  WBL  must  still  have 
non  zero  values.  This  is  necessary  so  that  the  computer  will  not 
divide  by  zero  which  would  cause  a  program  interrupt. 

5.  INSYM;  Format  (110).  INSYM  indicates  whether  a  torsional  analysis 
is  required.  If  the  numeral,  1,  is  placed  in  column  10  then  TORPAK 
calls  up  SYMPAK  and  the  torsional  analysis  is  suppressed.  When  a 
torsional  analysis  is  required  the  card  is  left  blank. 

6.  FX;  Format  (10F8.2).  FX  is  the  applied  lateral  load  in  the  Y 
direction.  The  loads  are  read  in  starting  with  the  bottom  story 
and  moving  to  the  top.  If  the  number  of  stories  is  greater  than  10, 
additional  cards  must  be  used. 

7.  FY;  Format  (10F8.2).  FY  is  the  applied  lateral  load  in  the  X 
direction.  Both  FX  and  FY  are  in  kips. 

8.  XH;  Format  (10F8.2),  XH  is  the  distance  from  the  reference  point 
to  load  FX, 

9.  YH;  Format  (10F8.2).  YH  is  the  distance  from  the  reference  point 
to  load  FY. 

10.  ITERX ,  ITERY;  Format  (2110).  ITER  is  the  number  of  iteration 
cycles  to  which  the  correction  technique  is  limited.  X  and  Y  re¬ 
fer  to  the  direction  in  which  the  correction  technique  is  to  be 
applied.  After  the  number  of  iterations  in  either  direction  has 
reached  the  limit  then  TORPAK  calls  up  CALPAK  which  calculates  the 
member  forces  on  the  basis  of  the  last  analysis. 
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11.  GKT;  Format  (5F16.6).  GKT  is  the  sum,  for  each  story,  of  the 
St.  Venant  torsion  constants.  The  units  are  kip  square  inches. 

12.  XX;  Format  (10F8.2).  XX  is  the  perpendicular  distance  between 
each  bent  and  the  reference  point.  In  this  and  the  following  card 
names  the  X  and  Y  refer  to  the  direction  in  which  the  measurement 
is  taken. 

13.  YY ;  Format  (10F8.2).  YY  is  the  perpendicular  distance  between  each 
bent  and  the  reference  point. 

14.  CX;  Format  (10F8.2).  CX  is  the  distance  between  the  reference  point 
and  the  centre  of  resistance  (El)  for  each  story. 

15.  CY ;  Format  (10F8.2).  CY  is  the  distance  between  the  reference  point 
and  the  centre  of  resistance  (El)  for  each  story. 


t  ■)  *  •>  r  .lc  ■  ’  .  t  -  . . 

k 


•  f  •  '  t.c  ■'  0  !  2.  VJ  ;{  r  )f  j6  !  O 


•£  s  tf-s  (I  i  s.  :  V?  inso  :  .tv 


’ 


. 


■  It  I ■ . 


non 


APPENDIX  C 


79 


PRINT-OUT  OF  FORTRAN  IV  PROGRAM 


c 

(..  CGMPAK  MAIN  PROGRAM 

c 

c 

r. 

COMMON  KJ(  2f>,2C  )  ,KP  (2  )  ,KS(  2'"  ,20)  ,KX  (4 ,2r',  20  ,KY<  3,20,20)  ,KV(  2C 

1 .20)  ,M  20  ,  PC  ),  NS  ,NX  ,  NY,  I  CX,  I  CY  ,  X0(  21"  )  ,Y0(  2\>)  ,KXROTW  (  4, 29, 20)  ,K  YKuT 
2C(3,20,2C),KYR0Tw<3,20,2'O,Hl20!,CI(2P),Wl(20),„BLI20),W(20),FBL(2 
30)  ,WBI(2?I  »E,LW,EGC»EGW,FX(20)  ,  F  Y  (  2  0  )  ,XX(4)  ,  YY  (  3 )  ,CX(2n),CY(2u)  ,KX 
4(4 ,20) ,UX( 4,20,20) ,KDX (4,20,20 ,TTX( 20 ,20 ) , ROT (20) , VX ( 4 , 20 ) , T V X ( 20 
5)  ,RY(  3,20)  ,UY(  3,20 ,2C  I  ,Kl)Y(  3,2°  ,2>'.  )  ,TTY  (2C  ,20  )  ,VY(  3,20)  ,  TVY(  20)  ,MT 
6(20  ,  TVXA(  20  ,TVYA(  20  >  ,XH(  20  >  ,  YH(  20)  ,  FBI  (  20  ,KXROTC  (4,20,20)  ,  INDIC 
7X(4)  ,  IN0ICY(3)  ,  H  X  ( 4  ,20  )  ,CIX(4, 20)  ,w»IX(4,2r)  *WBLX(4,2G)  ,  w  X  (  4 , 2  C  )  ,Fb 
BLX(4,20)»tof3IX(4»2O  ,FBIX(4,20),HY(3,20),CIY(3,2C),WIY(3,20),wBLY(3 
9,29) ,WY( 3,PO  ,  F  B  L  Y  (  3,2'“')  ,WBIY(3,20)  ,  F  B I  Y  (  3,20) 

COMMON  K  A  (  2  ,  2  "  )  ,KB(2  0,20),KD(2r  ,  20),KE(2r,2C),KC(20,20),KF(2C,20) 
1»KG( 20 ,20 ) ,KH(20, 20 ,Kl ( 20 ,29)  ,KK( 20,  20) ,KL(20,20)  ,KM (20,20) ,KN<  20 

2.20)  ,KfJ(  20,20)  ,K«  (2f  ,20)  ,  KR(2C  ,  20  )  ,KT(  2  0  20  )  ,  Kw  (20, 2C  )  ,  T  X  (  4, 20, 20  ) 

3*TY(  3,20*20  ,RT(  29,/CI  ,SX(4, 20  ,SY(  3, 20  )  ,K(JTXC(4,20)  ,R0TXh<4,20)  , 
4R0  T  YC ( 3 ,20 ,ROTYW( 3,2  ) ) , T T VXC ( 20 , 20 ) , T VXC ( 2 0 , 2C ) , T T V YC ( 20 , 20 ) , T VYC 
5<  20, 20)  ,CMXB(  4 ,2  3  )  ,CMXT(  4,20  )  ,  WMXB(  4, 20  )  ,  WMXT  (4 ,20)  ,FBMX(  4,  2C  )  ,  wtJM 
6XC ( 4,20 ), WBMXw( 4 ,20 ) ,CMY3( 3,20 ) ,CMYT( 3,20 ), WMYB ( 3,20 ,WMYT ( 3,20 ), F 
7BMY ( 3,20 ) , WBMYCI 3,2  )  ,  W6MYW  (  3 , 20  )  ,CSX  <  4 , 20  ,  WSX  (  4 , 20  )  ,  F  BSX  (  4 , 20  )  ,  W 

8BSX ( 4,2  0  )  ,CSY( 3 ,20  ,  wSY(  3 ,20) ,FBSY{ 3, 2C  )  , WBSY ( 3,20 ) ,CFX(4,20 I  , WFX( 

94 .20)  ,  FBFX(  4,20  )  ,„B)  X  (  4,2  o  )  ,CF  Y(  3 ,20'  (  ,  WFY  (  3,20)  ,FBF  Y(  3, 20  ) 

COMMON  wBFY ( 3,20 ) , OCX (4, 20) ,DEY ( 3,20) , INSYM, I TFRX, I TERY, I  X, I Y 
RFAL  KJ,KP,KS,KX,KY,KV,k, KXRU T C , KXRUT W , KY RO TC , NYR OT W , KOX.KOY ,MT 
REAL  KA , KB »  KD  »  Kfc  , KC , KF , KG , KH , K I , KK , KL , KM , KN , KU , KU, Kk , KT , KW 

C 

THIS  SECTION  R.FADS  IN  THE  BENT  MEMBER  PROPERTIES  AND  THEN 
CALCULATES  THE  BENT  STIFFNESS  MATRICES  (IN  TERMS  OE  STORY  SWAY) 

READ! 5, 100  INS, NX, NY 
130  FORMAT (3110) 

READ(5»99)E,EW,EGC,EG« 

99  F0RMAT(2FlO2,2E12.5) 

DO  101  N  = 1 »  N X 
READ(5, 117)1 NDICX(N) 

117  FORMAT  (  IK) 

I  F  C  INOICX(N)  .EU.l  )GU  TO  116 

115  CALL  RE  AD  1 
DO  120  1=1, NS 
DO  121  J  =  1 , N  S 
KXRCTC ( N,  1 , J ) =KJ ( 1 , J  ) 

KXROT  W ( N , I  ,  J)=C.C 

121  CONTINUE 
120  CONTINUE 

102  CONTINUE 

I  F  (  I  N  0  I  C  X'(  N(  .  E  U  •  0  )  G  U  TO  103 

116  CALL  READ  2 
DO  122  1=1, NS 
DO  123  J= 1 , NS 
KXRCTC ( N, 1  ,  J ) =  KP ( I , J ) 

K  X  ROT  w I N » I  ,  J ) =KS ( I , J ) 

123  CONTINUE 

122  CONI  1 NUE 

103  CONTINUE 
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DO  104  1=1, NS 
DO  105  J= 1 , NS 
KXIN, l » J)=K( I , J) 

105  CONTINUE 
104  CONTINUE 

DO  130  1=1, NS 
HX( N  * IJ  =H( 1 1 
c  txi n, i»=ci 1 1 1 

W I  X  (  N ,  I  )  =  W I  {  I  ) 

WBLXIN, I  I  =  WBL ( I ) 

WXIN.I )=W(  I  I 
FBLXI N,  I  ) =FBL ( I ) 

MB  I  X ( N , I ) =WB III) 

FBIXIN, I )  =  FB I ( I ) 

130  CONTINUE 
101  CONTINUE 

DO  106  N=1,NY 

READ15, 1171 INDICYINI 

IF( INDICYI N) .EQ. 1 IGO  TO  119 

118  CALL  RE ADI 

DO  124  1  =  1, NS 
DO  125  J= 1 , NS 
KYROTCI N, I ,J)=KJI I ,J ) 

KY  ROTW l N , I ,J»=0.0 

125  CONTINUE 
124  CONTINUE 

107  CONTINUE 

IF( INDICYINI .EQ. 0)00  TO  108 

119  CALL  READ  2 
DO  126  1=1, NS 
DO  127  J= 1 , NS 

KY  ROT  C I N , I , J I =KP (  I, J) 

KYROTW ( N, I  ,  J)=KS( I, J) 

127  CONTINUE 

126  CONTINUF 

108  CONTINUE 

DO  109  1=1, NS 
DO  110  J=1 , NS 
KYIN, I , J)=K( I , J) 

110  CONTINUE 

109  CONTINUE 

DO  131  1=1, NS 
H Y ( N , I  I =H I  I) 

Cl  YIN, I  )=CI I n 
W I Y I N  ,  I  I  =WI I  I ) 

WBLYIN, I )  =  WBL I  I  I 
WYI N » I ) =M I  1 1 
FBLYIN, I ) =FBL I  I  I 
WBIYI N, I )=WBI I  I  I 
FBIYIN,  I ) =  FB I  I  I  ) 

131  CONTINUE 
106  CONTINUE 

READ!  5,  1 1*7  I  INSYM 
CALL  TORPAK 

132  CONTINUE 
STOP 
END 

SUBROUTINE  READ1 

COMMON  KJI 20,20)  , KP I  20 ,20 ) ,KS<  Z0,20 ) ,KXI4 ,20,20) ,KY<  3,20,20)  ,KV( 20 
1  ,20)  ,KI 20,20) ,NS,NX,NY,ICX,ICY,X0(20) , Y 01 20 > , KXRO T W I  4 , 20 , 20 >  .KYROT 
2C I  3,20,20 ), KYROTW I  3,20 ,2 0 ) ,H I  20 ) , C I  I  20 ) , W I  I  20  I , WBLI 20 ) , W I  20 ) ,F BL I  2 
30 ) ,WBI I  20)  , E,EW, EGC ,EGW,FX(20)  ,FY(20),XXI4),YY(3),CXI20),CY(20),RX 
414,20)  ,DX I  4,20,20) , KUX 1 4 , 2 0 , 20 ) ,T TXI 20, 20 )  ,  ROT l 20) , VX I  4 ,20 ) , TVXI 20 
5) ,RY(3,20) ,DY(3»20,20I  ,KDY(3,20,2C)  ,TTY(20,20)  »VY( 3,20)  *  TVY  I 20) ,MT 
6(20) , TVXAI  20 ) ,TVYAI 20)  ,XHl 20) , YH( 20 ) ,FBI I  20 ) , KXROTC I  4 , 20 , 20 )  , INDIC 


,  I  c- 


- 
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■ 


. 


* 
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7X  (  4  )  ,  I NDICYI  3)  *HX ( 4  *20 )  ,  L  I  X(  4 , 2°  )  *  W  1  X  (  4. 20  |  ,  wBL  X  (  4  ,  ?C  J  , MX ( 4 , 20  J  *  F  B 
8LX(4,20),WBIX(4,20),FB1X(4,20),HY(3,2C),CIY<3,20),W1Y(3,20),WBLY(3 

9.20)  t MY ( 3,20  I ,  Hi  LY<  3,2°)  *  MbI Y( 3*20) , F HI Y( 3, 20  I 

COMMON  KA( 20 ,20  I  , K3 (2  0,20 ) ,KU( 2F ,20 ) , KE  <20, 20 ) ,KC ( 2U, 20) ,KF { 20, 20 ) 
1 ,KG( 20 , 20) ,  M(( 20  ,20) ,K t ( 20,20) ,KK ( 20,  20 ) ,KL ( 20, 20) ,KM( 20, 20) ,KN( 20 

2. 20 ) , K0( 20, 20), KU<20, 20), KR (2r  , 20 ),KT(20,20),KW(20,cC),TX(4, 20, 20) 
3, TY( 3,20  ,20) *  K  T ( 20,20  ,SX(4,20 )  ,SY(  3 , 20 ) , KOT XC ( 4 , 20 ) , ROTXM ( 4 , 20 ) , 
4R0TYC( 3,20) ,ROtYW( 3,20) , T  TVXC ( 2°, 20) , VV XC (2 0,20) ,TTVYC( 20 ,20) ,TVYC 
5| 20 ,20) , CMXBt  4  *20 ) ,CMXT( 4,20 ) ,WMXBl 4, 20 ) , WMXT  <  4,20 )  ,FBMX(4,20 ) , WBM 
6XC (4,20) , WBMXW( 4,20)  ,CMYB(  3 ,20 )  ,LMYT( 3,2C ) , WMYB ( 3,20) , WMYT( 3,20 ) , F 
7BMY( 3,20) , WBMYCt  3,20 ,MBMYW( 3, 20)  , CSX (4,2C ) ,WSX( 4, 20) ,FBSX(4, 20) , W 
8BSX(4,2P) ,CSY( 3,23) ,WSY( 3,20) , F  BS Y ( 3,20) ,WBSY(3,2C) ,CFX(4,20) ,  ME  X  ( 

94.20)  , FbF  X ( 4 , 20 ) , MBE  X ( 4 , 20 ) ,CFY( 3,20)  ,WFY( 3,20  ,FBFY( 3,20) 

COMMON  WBFY(3,20  )  ,DEX(4, 20) ,  ObY  (  3,20  )  ,  IIMSYM,  I  TfcRX,  I  TbRY,  I  X,  I  Y 
REAL  KJ,KP,KS  »KX,KY,KV*K,  KXROTC , KXROTW , KYRDTC  ,  K  YROT  M  ,  KUX  ,  KDY  ,  MT 
REAL  KA ,KB ,KD , KE , KC ,KF ,KG , KH , K I , KK, KL ,KM, KN, KO ,KU ,KR , KT ,KW 

DO  1^1  1=1, NS 

READ(5,102)H(I),CI(l),Wl(I),W(I),WbL(I),FBL(I),WBI(I),FBI(I) 

102  F0RMAT(F8.2,2F12.2,F8.2,4F1C.2) 
loi  CONTINUE 

DO  103  1=1, NS 
DO  104  J= 1 , NS 
K A ( I  ,  J)=0.0 
KB( I , j) =0.0 
KD ( I , J ) =0  - 0 
KE ( I , J) =0.0 

104  CONTINUE 

103  CONTINUE 

DO  105  1=1, NS 

KA(  I, I  )  =  -(  (12»C*E*CI(  I ) )/(H( 1) *H(  I)*H(  I ))  ) 

IF ( I  . EQ  . NS ) GO  TO  105 
N=  1  + 1 

K  A (  I  , N  )  =  (  (  12.C*E*CI <N )  )  / ( H ( N ) *H (  N  I  *H ( N  )  )  ) 

105  CONTINUE 
L=NS-1 

DO  106  1=1, L 
N=  I  + 1 

K  B  (  I  ,  I  )  =  6 • 0*E  *(  (CI(I)/(H(I)*H(I)))-(CI(N)/(H(N)*H(N)))) 

KB(  I ,N  )=-(  <6. C*E*CI (  N  )  ) / ( H ( N ) *H  (  N )  )  ! 

KB ( N , I ) = ( <6.P*E*CI (N) ) / ( H ( N ) *H ( N ) ) ) 

106  CONTINUE 

KB(NS,NS)=( (6.0*E*CI(NS))/(H(NS)*H(NS) ) ) 

DO  107  1=1, NS 

KD(  I , I  ) =-(  (6.0*E*C I (  I  )  1/ (H( I )*H( I ) ) ) 

I F ( I . EQ .NS ) GO  TO  107 
N=  I  + 1 

KD( I , N ) =- ( ( 6 . 0*E*C I (N) )/(H(N)*H(N) ) > 

107  CONTINUE 
L=NS-1 

DO  108  1  =  1, L 
N=  I  +  1 

KE (  I , I  ) =E* ( (4.0*0  I  ( I  ) /H(  I ) )  +<4.^*C I (N )/H( N) )  +  ( 12.0*FBl ( I )/FBL(  I ) )  ) 
KE ( IfN)=(2.0*E*CI(N) )/H(N) 

K E  (  N  ,  I  )  =  (*2 . 0*E*C  I  (  N  )  )/H(N) 

108  CONTINUE 

KE(NS,NS)=E*( (4.0*CI (NS) /H(NS)  )M  12.0*FBI(NS)/F8L(NS) ) ) 

CALL  MACAL1 

RETURN 

END 

SUBROUTINE  READ2 

COMMON  KJ( 20, 20 ) , KP (20,20 ) ,kS( 20,20) , KX ( 4, 20,20 ) ,KYl 3,20,20) ,KV( 20 

1.20)  ,K(  20,20  )  »NS  ,NX,NY,  I  CX.ICY,  XCJ(2C)  ,YU(  20)  ,  KXROTW  (  4 , 2  0 , 20  )  ,KYRUT 
2C(  3,20,20) ,KYROTW( 3,20 ,23)  ,H(20 ) ,C I ( 20 ) ,MI (20 ) ,MBL ( 20) ,M( 20 ) ,FBL(2 
30)  ,  MB  I  ( 20)  ,b,EW,EGC,EGW,FX(2C)  , F Y ( 20 ) , X  X ( 4 )  , Y Y ( 3 ) , CX ( 20 ) , CY ( 20 ) ,RX 
4(4,20)  ,DX( 4, 20,20 ) ,KDX (4,20, 20)  ,TTX( 20,20 ),  ROT (20) ,VX( 4,20) ,TVX( 20 
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5  )  ,R  Y(  3,20  )  ,  U  Y  (  3 ,2)  ,  20  )  ,KI)Y  (  3,  20  ,2'"'  )  ,  TT  Y  (  20  ,  2C  )  ,  VY(  3,2U  I  ,  1  VY(  20)  ,MT 
6(20)*TVXA(?0),TVYA(?C)  ,<U(20),YH(?(.  ),EBI(20)«K  XkOTL  <  4 , 20, 20)  ,  INlJl  C 
7X(4l,lNI>lLYm,tlX(4,20)»CIX(4,2i),WlX(4,2r),WtllX('t,2C).*X(4,20),FU 
8LX  (4,20  )  ,  WBI  X  (  4 ,2r  )  ,  F  H  1  X  (  >  ,20  )  ,  NY  (  3,  ?<"'  )  ,  C  l  Y  (  3  ,  20  I  ,  W  I  Y<  3 , 2  0  I  ,  WBLYI  3 
0,2°)  *  W  Y  (  3,20  ,  Fill  Y<  3,2")  .  WHl  Y  (  3 ,2  0)  ,  F  8  1  Y(  3,20) 

COMMON  K  A  (  2  r'  ,  2 1  )  ,  KBI  2  j  ,  L’l  )  ,  KU  (  2r' ,  20  )  ,  K fc  (  2 r  ♦  20  )  ,  KC  (  20 , 20  »  ,  M-  (  20 , 20  I 
1  »KG(  20 ,20)  ,  KH  (  2  0  » 2  f  )  *  K  I  (  2  "  *  2C  )  ,KK(  2-1,20  )  *  K  L  (2  0,20)  ,  KM ( 2  0  *  20 )  » KN ( 20 
2,20,  MU  20 , 20  )»KU(2r,2-l)  *  KK (  2C  »  21 ) »  KT  (2°,20),KW(20,20),  1X14,20,20) 
3»TY(3,2'',2O),RT(2T,?0>,SX(4,2Ct,SY(3,20>,RUTXU4,20),kOTX*(4,20>, 
4PU  T YC (  3,2°)  *  kOT  Y W (  3,20)  ,  TTVXCI  20*20)  *TVXC(20,20)  ,TTVYC(2C,2C)  »  T  VYL 
b(  20,20)  ,CMXB(4,20  J  ,  CM  XT  (  4 ,27  )  ,  WMXIi(4, 20)  ,  wMXI  (4,20)  ,  FBMX  (  4,20  )  ,  WbM 

6  XC ( 4 , 2r  )  , WBMXW(4  ,2  0 )  ,  CMYBI  3,2" )  ,CMY  T ( 3,20  I , WMYR ( 3 ,20 ) , WMYT ( 3,20) , 7 
7BMY(  3,2")  ,  WBMYC.l  3,21)  ,  WKMYw(  3,;:'0  ,CSX(4 , 20  )  ,WSX(  4,20)  ,FBSX(  4, 20)  ,  W 
PHSX(4,20) ,TSY( 3,20) »  W  S  Y  ( 3*2") ,  F  13  S  Y  ( 3,20) ,WBSY( 3,20) ,CFX(4,20) ,  W  F  X  ( 
04, 2r )  ,)  BF  X(  4  ,  ?(  )  ,  W  0 )  X  (  4 , 2  1 1  ,LFY(3,2r  )  ,wFY(  3,20)  ,FBFY(  3,20) 

COMMON  wB(  Y  l  3,2  0  )  ,  0  [.  X  (  4 , 2  3  )  ,  BE  Y  <  3 , 20  )  ,  l  N S  YM  ,  1  T LRX ,  I  T E R  Y  ,  I  X  ,  I  Y 
P  r  A I  K  J  ,Kl  ,KS  ,  Kx  ,  KY  ,  KV  ,K  ,  K  XRUTC  ,  KXP  NT  w  ,  KYRUTC  ,  KYRUT  W,  KOX,  KDY  ,  MT 
REAL  K  A  ,  K.B  ,  KIJ ,  Kt  ,  K  C. ,  K  F  ,  K(j  ,  KU  ,  K.  I  ,KK,KL  ,KM,KN,Kl),KU,KK,4T  ,KW 
DO  111  1=1,  Nli 

RE  ADI  0,1  "2  )H(  I  )  ,CI(  I  )  ,  Wl  (  I  )  ,W(  I  )  ,WBI  (  I  )  ,  FBL  (  I  )  ,  WBI  (  1  )  ,  FBI  (  I  I 
K  2  FORMAT ( Fh. 2 ,2F 12 .2, E  8.2 , 4F  U  . 2 ) 
lrl  CUNT  [Mill 

[)' '  PI  1  =  1,  '):■> 

1)0  114  j  =  l,ris 
K  A (  I , J )  =i.O 
KB ( 1 , J) =0.0 
KC ( I , J) =0.1 
KD ( 1 , J)=o.O 
K  E (  1 , J )  -I." 

KF( T , J) =0.0 
KG ( I , J) =0.0 
KH(  1 , J)  =0.0 
K  1  (  I  ,  J  )  =0  .  o 
1 r  4  CONTINUE 
103  CONTINUE 

DO  1°  5  1=1,  NS 

KA(  I ,  I  )=-(  ( 12 .0) / l HI  I ) *H<  I  )*H(  I  ) >  ) *<  <  E*C I (I) )  +  ( fcW*WI (I)  ) ) 

IF ( I .EQ.NS ) GU  TO  l"b 
N  =  I  ♦  1 

K A (  I  , N )  =  ( (1 2.0 ) / (H( N) *H( N  )  *H (N  )  )  )  *( ( E*C I l N )  )  +  ( E w*W I (N )  )  ) 

If  5  CONTINUE 
L=NS-1 
DO  l?t>  1  =  1,1 
N=  I  +  1 

KB  (  I  ,  I  )  =  6 . 0*E  *  < 1  Cl  (  I  )  /  (  H  (  I  )  *H(  I  >  )  )■-  t  C  I  (  N  )  /  (  H(  N  )  *H(  N  )  )  )  ) 

KB(  I  »N)=-(  (6.r)/(H(N)*H(N)  )  )  *  (  E  *C  I  (  N  )  ) 

K  B  (  N  ,  1  )  =  (  (  6. 0)  /(  HI  N)  *H(N>  )  )* ( E*L  I  <  N ) ) 
lie  CONTINUE 

KB(NS,NS)-((6.n>/(H(NS)*H(NS> ) ) *( E*C I ( NS ) ) 

1C  1  CONTINUE 
L=NS-1 

DO  1 08  1=1 , L 

N=  I  +  1 

KC(  1 , I  )=6.0*EW*(  <W1  (I)/( H(  I  >*H(  I  I ) )-(WI (N)/IH(N)*H(N) ) )  ) 

KC ( I , N ) =— ( (6.0)/ ( H ( N ) *H ( N ) ) )*(EW*Wl(N) ) 

KC(N,I  )  =  (  (6.0)/(F)(N)  *H  (N)  I  )*(EW*WI(N)  ) 

If  8  CONTINUE 

KC(NS,NS)  =  ( ( 6.0  » /(H(NS)*H(NS) )  ) *( EW*W I (NS ) ) 

DO  110  1=1 , N  S 

KD(  I  ,  I  )=-(  (6.0) /(H(  I ) *H(  l)> )*( E*C  l(  I )  > 

KG  (  I  ,1  )=-(6.0*EW*Wl  (I  )  )/()(<  I  )*H(  1)1 

IF ( 1 .EQ.NS )G0  TO  110 

N=m 

KD(  I,N)=-(  (6.0)/(H(N)*H(N) )  ) * ( fc  *C  1  ( N )  ) 
KG(I,N)=-(6.0*EW*WI(N))/(H(N)*H(N)) 

110  CONTINUE 
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L  =  N  S  -  l 

DO  111  I = 1  * L 
N  =  I  +  1 

K E  (  l  ,  I  )  =E  *  (  (  4 . 0*  t  C  I  (  l  )  /i  I  (  I  ) )  )  ♦  <  4 .  <"  *  ( (  l  (  N  )  /H  (  N )  )  )  ♦  (  1 2 . 0  *  (  F  6 1 1  I  )  /  F  B  L 
1(1  ) ) ) + ( 4 . C  * ( WRI ( ! ) / WBL ( 1  )  )  )  ) 

KF  <  I ,  N  )  =  E*  ( (2  •  0  *  C  I  ( N )  l  /H(  N)  ) 

KC ( N  » I  )  =  E  *( ( 2 . 0*C I <  N )  > /HEN)  ) 

111  CONTINUE 

KE(NS,NS)  =  E*tl4.0*(Ci(N3)/H(NS>>)M12.'*(FBIlNS)/FBL(NS)>)M4.0* 

1  (WBI  (  NS  )  / WBL ( N S )  I  )  ) 

DO  113  1=1, NS  * 

KF(I,I)=((F*WBl<  m/W!3L(  I  ))*(<2.n)+(3.n*w(  I)/WBL(  1))) 

KH(  I  ,  I  )  =  (  (E*WB1  (  I  )  )/WBL(  I  )  )  *  (  (  2 . 0  )  F  (  3 . 0*W (  I  )  / WBL (  I  1  )  ) 

113  CONTINUE 
L  =  NS- 1 

DO  114  1  = 1 , L 
N=  I  + 1 

K  I  (1  ,  I  )  =  (  (  (E*WBI  (  I  )  >/WBL  (  m  *(  (  4.  P)+(  6.C*W(  I  )  /WBL  (  I  )  )  +  (  3.0*W(  I  ) 

1  *W  (  I  )  /  (  Wt3L  (  I  )  *WBL ( I  )  )  )  I  )  *■  (  (4.0  *EW*W  I  (  I)  )  /  H  C  [  ))+<(4.0*Ew*Wl(N)  )  /  H  (  N 
2)  ) 

K I ( I , N ) =EW*( ( 2.0*W I ( N) ) / H ( N ) ) 

KI (N, 1  )  =  EW  * ( ( 2 • 0  *W I (N)  ) / H ( N ) ) 

114  CONTINUE 

KI  (NS  ,NSl  =  (  t  (  E*WBI  (  NS  )  )/WBL  (NS  )  )  *(  <4.C)+(  6.C*W(N's  )/WBL(  NS)  )  +  (  3.0* 
1W(NS)*W  (NS)  /(  WBL  (NS  )  *WBL  (  NS  ))  )  I)  +  (  (  4.  C  *EW  *W  I  (  NS  )  )/H(NS)  I 
CALL  MAC AL  2 
RETURN 
END 

SUBROUTINE  MAC  AL 1 

COMMON  KJ( 20,2" >  ,KP( 20,20 ) ,KS<  29,20  > ,KX( 4, 20,20 ), KYI  3,20, 2G)  ,KV( 20 

1 .20)  , K( 20*20 ) , NS, NX, NY, I CX, ICY ,X0(20 ) ,YOt  20) , KXROTW < 4 , 20 , 20 )  ,KYROT 
2C(  3,20,20)  ,  K  YROT  W(  3 , 2  0  ,2  0  )  ,H(  20  )  ,  C  I  (  2f  )  ,  w  I  ( 20  )  ,  WBL  (  20  )  ,  W  (  20  )  ,  FBL  (  2 
30) ,wBl<  20)  .E,EW,EGC,EGW,FX(20) ,FY(2D  ,XX(4)  ,  Y  Y ( 3 )  , C X ( 2?  )  , C Y ( 20 ) ,KX 
4(4 ,20  ,DX( 4,2C ,20)  ,KDX(4 ,2  0,2C )  ,TTX(20,20) ,kQT( 20) , VX(4,20) , TVX120 
5)  ,  RY  (  3 , 20  )  ,UY(3*2T,29)  ,KDY(3,2Q,2C)  ,  T  T  Y  (  2C‘  ,  20  )  ,VY(3,2  0)  ,TVY(20)  »MT 
6(20)  ,  TVXA(  £§  )  ,TVY4(  20  )  ,XH(  20)  ,  YH(  20  )  ,FBI  (  20  , KXROTC < 4 , 20 , 20 )  ,  INDIC 
7X(4),INDICY(3),HX(4,20),CIX(4,2n),WlX(4,20),WBLX(4,2C),WX(4,20),FB 
8LX(4,20) ,WBIX(4,20) , FBI X( 4,20 , HY<3»20) ,CIY(3,20) ,WIY(3,20) , WBLY(3 
9,2") , W  Y ( 3 , 20 ) ,FbLY(3,20)  ,wBIY(3,2C),FBIY(  3,20) 

COMMON  KA( 20, 2r ), KB (20,70 ) ,KD( 20,20 ) ,KE( 20,20 ) ,KC( 20,20) ,KF( 20,20) 
1,KG( 20,20)  ,KH( 20,2")  ,KI ( 20,20 ) ,KK( 20,20 ) ,KL (20,20) , KM (20,20 ) »KN (20 

2.20 )  ,K0( 20,20 ), KOI  20 ,29) ,KR( 2C ,2r ) ,KT( 20, 20  I ,KW (20,20 ,TX( 4, 20,20 ) 
3,TY( 3,20,20) ,RT( 20,20)  ,SX( 4,20  ,SY( 3, 20 ) ,ROTXC( 4,20) ,ROTXW( 4,20) , 
4R0TYCI  3,20  *ROTYW(  3,21)  ,  TTVXC120 ,20  )  ,  TVXC  (  20' ,  20  )  ,  T  TV  YC  (  20 , 20  )  ,TVYC 
5(  2"  ,20)  ,CMXB(  4,2  0)  ,CM  XT  (4,20  )  ,WMXB<  4,20  )  ,WMXT(  4,20  ,FBMX  (4,20)  ,  WBM 
6XC(4,20) , WBMXW( 4,20 ) ,CMY6( 3 ,20 ) ,CMYT{ 3,20 ) , WMYB<3,20) ,WMYT( 3,20) ,F 
7BMY( 3,2n) ,WBMYC(  3,2'  )  ,  wBM  Y W ( 3,20)  , CSX  <  4 , 2(  ) , W  SX ( 4 , 20 ) ,FBSX(4,20) , W 
8BSX(  4,'2")  ,CSY(  3, 2C  )  ,WSY(  3,20)  ,FBSY(  3,29  )  ,  WBSY(  3,20)  ,CFX(4,20)  ,  WFX( 

94. 20)  ,  FBFX(  4,29  )  ,  WBF  X  (4,20  )  ,CFY(  3 ,20)  ,  WFY(  3 ,20)  ,FBFY(  3,  20) 

COMMON  WBF  Y(  3, 20  )  ,i)tX(4,20  )  ,DEY(  3,20)  ,  INSYM,  ITERX,  ITERY,  I  X,  I  Y 
REAL  KJ,KR,KS,KX,KY,KV,K,KXRUTC, KXROTW, KYROTC,KYROTW,KDX,KDY,MT 
REAL  KA,KB,KD,KE,KC,KF,KG,KH,KI , KK , KL , K M , KN , KO , KU , KR , KT , KW 
CALL  T U K Y (J  (  N S  ,  K L  ) 

1^0  CONTINUE 

DC  101  1=1, NS 
DO  102  J= 1 , NS 
KJ ( I , J)=0.0 
DO  103  L= 1 , NS 

K J ( I , J ) =KJ ( I,J)-KE(I,L)*KD(U,J) 

1 ° 3  CONTINUE 
102  CONTINUE 
ir 1  CONTINUE 

DO  104  1=1, NS 
DO  105  J=1 ,NS 
KK (  I  ,  J  )  =0.0 
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DU  106  L  =  l  *  NS 

KK(  1  »  J ) =KK (  I  *J  I  »  L ) *K  J ( L , J  ) 

106  CONTINUE 

105  CONTINUE 

104  CONTINUE 

DO  107  1=1, NS 
DU  108  J=  1 , NS 
K(  l , J ) =KA (  I  , J)+KK( I , J  ) 

108  CONTINUE 

107  CONTINUE 
RETURN 
END 

SUBROUTINE  MAC  AL2 

COMMON  KJ( 20,20) ,KP( 20,20) ,KSI 20,201 »KX( 4, 20,20 ) .KY( 3,20,201 ,KV( 20 

1 .20 )  ,K( 20,20) , NS  ,  NX, NY, I CX, ICY , X0< 20) ,Y0( 20) ,  KX ROT W ( 4 , 20 , 20 ) , KYRUT 
2C( 3,20,20) »  K  YROT  W (  3,20,20)  ,H(20)  , C  I  C  2  0  )  ,WI  (20)  ,WBL(20)  ,  w  I  20 )  *F  BL  ( 2 
30)  »WBI(20)  *E*EW*EGC»EGW*FX(20) , F  Y l 2  0 ) ,XX(4) ,YY(3) ,CX(20),CYI20) ,RX 
4(4,20) ,DX( 4, 20,20) ,KDX( 4, 20, 201 ,TTX( 20,20) , ROT  I  20 ) , VX ( 4 , 20 ) , T VX ( 20 
5) ,RY( 3,20)  ,DY( 3,20,20  > ,KDY( 3,20,20  I , TTY l 2^,20) , VY( 3,20)  ,TVY( 20) ,MT 
6(20), TVX  A  ( 20  )  ,TVYA( 20  )  ,XH (20) . YH( 20) ,ER t ( 20 > , K XROTC ( 4 , 2 0 , 20 )  , IND1C 
7X(4I , (NDICYJ  3 ) ,HX(4 ,20) ,CI X(4,20 ) , Wl X(4,2C ) ,WBLX( 4, 20 ) , WX(4, 20 ) ,F8 
8LX( 4,20 ) , WB IX ( 4,20)  ,FBI X( 4,2  0) ,HY<  3,20 ) ,C I Y( 3,20) ,M I Y( 3,20) , WBL Y ( 3 

9.20)  ,WY( 3,20) ,FBLY( 3,20 ,WB1Y( 3,20)  , FBI Y( 3, 20  ) 

COMMON  KA(  20,20)  ,  KB ( 2  0,20 ) ,KD<  20,20  > , Kb ( 20, 2C ) ,KC(  20, 20  > ,KF ( 20,20) 
1*KG( 20, 20) ,KH(20,20> ,KI ( 20,20) , KK ( 20*20) ,KL(20,20) , KM (2 0,20) ,KN(20 

2.20) »KO(20,20),KQ(20*20)»KRC20t20),KT(20*20)»KW(20*2C),TX(4,20»20) 
3, TY( 3,20,20) ,RT( 20,20} ,SX( 4*20) ,SY( 3,20 ) ,R0TXC(4,20) , R0TXW(4  » 20 ) , 
4R0TYC( 3,20 ) ,ROTYW( 3,20 ) , TTVXC( 20 ,20 ) , T V XC ( 20 , 20 ) , T T V YC ( 20  ,20 ) , T VYC 
5(20,20) »CMXB( 4,20) ,CMXT ( 4,20) »WMXB(4*29) »WMXT (4,20) ,FBMX(4,20) ,WBM 
6XC(  4,20 ) ,WB MX W( 4,20 ) , CMY B ( 3 , 20 » , C MY T ( 3 , 20 ) ,WMYB( 3,20) , WMYT ( 3 , 20 ) , F 
7BMY!  3,20)  ,  WBMYCI  3,20),V*8MYW(3,20),CSX(4,20),WSX(4,20),FBSX(4,20),W 
8BSX(4,20) ,CSY( 3,20) ,WSY<  3,20) ,FBSY( 3, 20) ,WBSY( 3,20) ,CFX( 4,20) ,WFX( 

94.20)  , FBFXi 4,20) ,WBFX (4,20) ,CFY<  3,20) ,WFY( 3,20) ,FBFY( 3, 20) 

COMMON  WBFY(3,20) ,DEX( 4,20 ,DEY( 3,20) , I  NS YM , I TERX , I TE R Y , I  X , I Y 
REAL  KJ  ,KP,KS,KX,KY,KV,K , KXRUTC , KXROT W , KY RO TC , K YR OT W , KDX, KDY , MT 
REAL  KA»KB»KD*KE,KC,KF*KG*KH,K1  ,KK,KL»KM,KN*KO»KQ»KR*KT  *  KW 
CALL  TOKYO  (NS,KI) 

100  CONTINUE 

DO  101  1=1, NS 
DU  102  J= 1 , NS 
K J ( I , J ) =0 . 0 
DO  103  L= 1 , NS 

KJ( I , J)=KJ( I,J)+KI ( I ,L)*KG(L, J) 

103  CONTINUE 
102  CONTINUE 

101  CONTINUE 

DO  104  1=1, NS 
DO  105  J=1 *NS 
KK ( I , J ) =0 .0 
DO  106  L= 1 »  NS 

KK( I, J)=KK( I,J)+KF( I,L)*KJ(L,J) 

106  CONTINUE 

105  CONTINUE • 

104  CONTINUE 

DO  107  1=1, NS 
DO  108  J= 1 , NS 
KL ( I , J)=KD( I , J ) -KK ( I, J) 

108  CONTINUE 

107  CONTINUE 

DO  109  1=1, NS 
DO  110  J= 1 , NS 
KM ( I , J ) =0 . 0 
DO  111  L= 1  *  NS 

KM(I,J)=KM(I,J)+KI(I,L)*KH(L,J) 

111  CONTINUE 
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1C  9  CONTINUE 

DO  1 1  2  1=1  , NS 
DO  1 1 3  J  =  1 ,  NS 
KNl I  ,  J) =0.0 
DU  114  L  =  1  *  NS 

KNU  ,  J)=KN  C  I  »  J  )  +  KF  (  1  ,  L)*KM(L  ,  J  ) 

114  CONTINUE 
113  CONTINUE 
112  CONTINUF 

DO  115  1=1, NS 
DU  116  J= 1 , NS 
KOI  I  ,  J  )  =KN ( 1 , J )-KE(  I , J ) 

116  CONTINUE 

115  CONTINUE 

CALL  TOKYO  (NS,KO) 

DO  117  1=1, NS 
DO  118  J  =  1  ,NS 
KP(  I  ,  J  1=0.0 
00  119  L= 1 , NS 

KP(I,J)=KP(I,J)+KG(I,L)*KL(L,J) 

119  CONT iNUt 
118  CONTINUF 

117  CONTINUE 

DO  120  1=1, NS 
DO  121  J  =  1  , NS 
KQ( I , J) =0.0 
DO  122  L= 1 , NS 

KQ<  I  *  J)  =KQ(  I  ,  JJ  +  KHU  ,L  )*KP(L,  J  ) 

122  CONTINUE 
121  CONTINUE 

120  CONTINUE 

DO  123  1=1, NS 
DO  124  J= 1 , NS 
KR (  I , J ) =KG ( I , J ) +KC(  I , J  I 

124  CONTINUE 

123  CONTINUE 

DO  125  1=1, NS 
DO  126  J  =  1  ,  NS 
K  S ( I ,J)=n.o 
DO  127  L= 1 , NS 

KS(  I  ,J)=KS(  I  »  J ) - ( KI  ( I ,  L)*KR<L, J) ) 

127  CONTINUE 
126  CONTINUE 

125  CONTINUE 

DO  128  1=1, NS 
DO  129  J= 1 , NS 
KT ( I,J)=0*0 
DO  130  L=1 , NS 

KTI  I  ,J)  =  KT  <  I  ,  J)«-KBl  I  , L ) * K P ( L , J ) 

130  CONTINUE 
129  CONTINUE 

128  CONTINUE 

DO  131  I=»1,NS 
DO  132  J= 1 , NS 
KW ( I , J)=0.0 
DO  133  L  = 1 , N S 

KW(I,J)=Kki(I,J)  +  KC(I,L)*KS(L,J) 

133  CONTINUE 
132  CONTINUE 

131  CONTINUE 

DO  134  1=1, NS 
DO  135  J= 1 , NS 

K(  I , J )=KA<  I , J ) *KT(  I , J  )  +KW1  I , J ) 

135  CONTINUE 

134  CONTINUE 
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RETURN 

END 


SUBROUTINE  TOKYO  (  N  0  ,  » 

DIMENSION  AC (2C )  *  Art ( 20 ) , A ( 20 , 2P  ) 

NO  1  =NtJ-  1 

At  I , 1 )  =  1 . O / A (  1,1) 

DO  80  N  =  l,  NO  1 
DU  50  1  =  1,  N 
ABC  I )=n.O 
AC( I  1=0.0 
DO  50  J=  1 ,  N 

ABI  1  I  =  A  B (  1  )  +  A  (  I  ,  J ) *A ( J ,N+  1  ) 

AC(I)=AC(I)*A(N+1,J)*A(J,I> 

50  CONTINUE 
ACB=o .0 
DO  60  1=1, N 
ACB=ACB+AC (  I  )  * A (  I ,N+  1  I 
60  CONTINUE 

A( N+l, N+l )=1 .0/1 A  (N+l ,N+ 1 ) -ACb ) 

DO  70  1=1, N 

A ( N+  1  ,  I  )  =- A  < N+l ,  N+l ) *AC(  1  ) 

A( I » N+l )  =  -AB ( I )*A(N+1»N+1) 

70  CONTINUE 
DO  80  I = 1 , N 
DO  80  J= 1 , N 

A(  I , J |  =  A(  I , J ) -A  I  I ,N+1 ) *AC ( J) 

80  CONTINUE 
RETURN 
END 

SUBROUTINE  TORPAK 

COMMON  KJ (2Q» 20) ,KP ( 20,29 ) ,KS( 20*20 ) ,KX (4 ,20,20 ) ,KY( 3,20,20) ,KV( 20 

1 .20)  ,K( 20,20) ,NS,NX,NY, 1 CX, ICY ,XU( 20) ,Y0( 20) , KX ROTW < 4 , 20 , 20 ) , KYRQT 
2C(  3,20,20) *KYROTW( 3,20,20) ,H(20) ,C I (20  » ,WI (20i ,WBL (2C ) ,W(20) ,FBL( 2 
3Q)»W6I(2C’  )»E,Eto*EGC,EGW*FX(29),FY(20),XX(4)»YY(3)*CX(20)*CY(2Q)»RX 
4(4,20)  ,0X14,20,20)  , KDX ( 4 , 20 , 2C )  , T TX ( 20 , 20 ) , ROT ( 20 ) , V X ( 4 , 20 ) , T VX ( 20 
5 ), RY( 3,20) ,DY( 3,20 ,2C  > ,KDY( 3,20 ,20 ,TTY( 20,20 ) ,VY{ 3,20)  ,TVY( 20)  ,MT 
6(20)  , TVXA( 20 )  ,  TV  YA ( 2  0 )  ,XH ( 20 , YH( 20) ,F8I(20),KXROTC(4,20,20)  ,INDIC 
7X(4),INDICY(3),HX(4,20),CIX(4,2C),WlX(4,20),W8LX(4,2C),WX(4,20),Ed 
8LX  (4,20)  ,WBIX(4, 20)  ,FB.IX  (4,20)  ,HY(3,2C  )  ,CIY  (3,201  ,*mi3,20>  ,  WBLY(  3 

9.20)  ,WY( 3,20 ,FBLY( 3,20) ,WBIY( 3,20) ,FBIY( 3,20) 

COMMON  KAI 20,20)  , KB (20,20 ) ,KD( 20,29 ), Kb (20, 2C ) ,KC I  20,20) ,KF( 20, 20) 
1, KG (20, 20) ,KH( 20,20 ) ,KI (20,20) ,KK( 20,20 ) ,KL(2C, 20) , KM (20,20 1 ,KN( 20 

2.20)  ,K0( 20,20 ) ,K«( 20,20 ,KR( 20,20 ) ,KT (20,20) ,KW( 20,20) ,TX( 4, 20, 20) 
3,TY( 3, 20,20) ,RT( 20,20) ,SX (4, 2C ) ,SY( 3,20) ,kOTXC( 4,20) ,KOTXW( 4,20) , 
4R0T YC ( 3,20  ,ROTYW( 3,20) , TTVXC( 20,20  I ,TVXC( 2C ,20) ,TTVYC( 20,20 ) ,TVYC 
5( 2n, 20) ,CMXb( 4,20) ,CMXT( 4,20) ,WMXB( 4,20) , WMXT (4,20) ,FBMX( 4,20) , WBM 
bXC( 4 ,20 ) , WBMXw(4 ,20) ,CMYB ( 3,2C ) ,CMYT( 3,20 ) , WMYB ( 3,20 ) ,WMYT( 3,20 ) ,F 
7BMY( 3,26) , WBMYC ( 3 , 2f ) , WBMYW ( 3 , 20 ) ,C SX ( 4 , 20 ) , WSX ( 4 , 20 ) , FBSX ( 4 , 20 ), W 
8BSX(4,20) ,CSY( 3,20) ,HSY( 3 ,20) , FbSY( 3,20) ,WBSY( 3,20) ,CFX(4,20) , WFX( 

94.20)  ,FBFX(4,20  )  ,WBFX(4,20)  ,CF  Y(  3, 20)  ,h’FY(  3,20)  ,FBFY(  3,  20) 

COMMON  WBFY(3,20),DEX(4,2"I),DEY(3,2T),INSYM,ITERX,ITERY,IX,IY 
REAL  KJ ,KP  ,KS ,KX,KY, KV ,K , KXROTC , KXPUT W , KYROTC , K YROTW , KDX.KDY ,MT 
REAL  KA  ,KB,KD,KE  ,KC  ,KF,KG,KH,KI  , KK , KL , KM , KN , KU , KU , KR , KT , KW 
DIMENSION  GKT <2n ) 

READ! 5,201  )(FX(  I  )  ,1=1 , NS) 

READ( 5,201) (FY( I ) , 1  =  1  ,NS) 

RE AD( 5, 201 ) ( XH( I ) , 1  =  1 ,NS  ) 

REA 0(5, 201 )(YH< I  )  , I  =  1 , NS ) 

201  FORMAT (  10F8.2) 

RE  AD ( 5 , 2 1 2 ) ITERX,  I TERY 
212  FURMAT(2I 10) 

READ! 5,203) ( GKT (  I ) , I  =  1 , NS  ) 

2n3  FORMAT! 5 F16.6) 

DO  210  1=1, NS 
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DU  211  J=1,NS 
KVl  I  » J1  =0.0 
211  C  UNT  1  NUT 
2  in  CONTINUE 
I  NS-l 

O'  1  202  I  =1  ,L 

N  --  1  ♦  1 

K.  V  (  1  ,  1  >  =  +  (  IbKT  <  1  )  /  H  (  1  )  )  ♦  (  GKt  (N)/HtN>  )  ) 

KVIN,  I  )  -NKT  (  Ml  /H  (  N  ) 

KVl I  ,N  )  -  -GKT  <  N )  /  H  (  N  > 

2°  2  CONTINUE 

K  V  I  N  S  *  N  S  ) =  +  I  G  R  T  I NS) /  H  I  N  3  )  ) 

Kt ADI  5, 2C1  ) ( XXI  I  )  (1=1  ,  N  X  ) 

READ!  5,201  XYYI  I  )  ,  1  =  1  ,NY> 
kLACI 5, 2C 1  )  (CXI  1  ) , I  =  1 , NS ) 

RE  ADI  5, 201  I  I C Y I  I ) , 1  =  l ,NS ) 

CALL  WRIT!  1 

WRITE  (6 ,  1 0  n  )  (GKT  I  1  )  ,  1  =  1, NS) 

ICO  FORMAT  I iHl , (2X, • ST. VENANT  CONSTANT  =  • , T E 1 3 . 6 ) ) 

WR  I  T  E  I  6  »  1 r  1  )  (t,EW) 

K  l  FORMAT  I  1HC’/,2X,  •  MODULUS  OF  ELASTICITY  FOR  FRAME  =  •  ,  F  10 . 2  *  2  X  ,  •  MOUU 
1 L US  OF  ELASTICITY  EUR  WALL  =  ',Fin.2l 
2  2  3  CONTINUE 

IF! INSYM.tU.niGU  10  254 
CALL  SYMEAK 
254  CONT 1 NUE 

I F I INSYM.EU. 1 ) GO  TO  252 
ICX  =  0 

i  cy=o 

I  X=n 

I  Y  =0 

IF ( F  XI  1  )  . FQ. n. >GC  TO  5r  1 
D  J  502  I =1 ,NS 
XOI I )=CX( I )*1 .2 
5°2  CONT INUE 
GO  TO  5"3 
5T1  DO  5rt  3  1  =  1 , NS 
XOI 1 )=XX(NX>/2. 

503  CONTINUE 

IFIFYI I  I . t  0 . 0 . ) G  0  TO  6°1 
DO  6 02  1=1 , NS 
YOI  I  ) =CY I  I  >*1.2 
602  CONTINUE 
GO  TO  603 
6C  L  DO  60  3  1  =  1  , NS 
YOI  I  ) =  Y Y I NY1/2. 

6° 3  CONTINUE 

DO  604  1=1, NS 

IFICXl I ).LE. (0.55*XX(NX> ) )GU  TO  605 
604  CONTINUE 
GO  TO  brb 
60  5  CONTINUE 

DO  607  1=1, NS 

MT I  I  )  =  I  F  X<  I) *C . 8 *XH I  I ) )-( FYI  I  ) *D. 8*YH (  I ) ) 

6r 7  CONTINUE 
GO  TO  224 
606  CONTINUE 

DO  229  1=1, NS 

MTI  I  )  =  I F  X  <  I  )*XH(  I  )  ) -I FYI I> *YH (II) 

229  CONTINUE 

224  CONTINUE 
CALL  HXODX 

225  CONTINUE 
CALL  HYODY 

226  CONTINUE 
CALL  DXOTTX 
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22 7  CONTINUE 

CALL  OYOTTY 
220  CONTINUE 

CALL  RTOROT 

232  CONT INUE 
CALL  VXOTVX 

233  CONTINUE 
CALL  VYOTVY 

234  CONTINUE 

DO  2  35  I  - 1 f  NS 
TVX A ( I ) =  TVX( I ) 

T V Y A (  I  )  =T VY ( I  I 

235  CONTINUE 

IF< ICX.EQ.  ITERXIGO  TO  255 

213  CONTINUE 

IF ( ICY. EQ. ITERY ) GO  TO  255 

214  CONTINUE 

WRITE (6,21 5) (TVXA(I) ,1=1, NS) 

215  FORMAT (  1H0,2X,' TVXA=' ,9E13.6) 

C 

r  THIS  SECTION  TESTS  FOR  EQUILIBRIUM  OF  LOADS  AND  RESISTANCE 

C 

I  F ( F  X ( ll.EQ.OIGO  TO  237 
DO  236  1=1, NS 

IF ( ABS( FX(  I) +TVXA( I )  I  .GT . ABS(0.ni«FX( I)  ) I  GO  TO  231 

236  CONTINUE 
GO  TO  238 

237  DO  238  1=1, NS 
IF(ABS(TVXA(  m.GT.O.lOIGO  TO  231 

238  CONTINUE 
I F ( F Y ( 1 ) . EQ . 0 ) GO  TO  239 

240  CONTINUE 
DO  241  1=1. NS 

TFIABSl  FY(  II+TVYAd)  ).GT.ABS(0.01*FY(  imGO  TO  230 

241  CONTINUE 
GO  TO  252 

239  DO  242  1=1, NS 
I F { ABS ( TV  Y  A  ( I ) ) .GT .0. 10 ) GO  TO  230 

242  CONTINUE 
GO  TO  252 

THE  FOLLOWING  SECTION  CALCULATES  INFLUENCE  COEFFICIENTS 

231  CALL  ICXO 
ICX=ICX+1 

243  CONTINUE 

WR ITE(  6,304) ( XO (  I  )  , I  =  1 , NS ) 

304  FORMAT! 1HO,2X,'XO=’ ,9E13. 6) 

GO  TO  224 
23C  CALL  ICYO 
ICY=ICY  +  1, 

244  CONTINUE 

WRITE!  6,305) (  Y0(  I)  ,  1  =  1,  NS  ) 

3C5  FORMAT! 1H0,2X,'Y0=' ,9E13.6) 

GO  TO  225 

255  CONTINUE 
WRITE! 6,256) 

256  FORMAT ( 1H1 ,4X, ' WARNING,  NO. OF  ITERATION  CYCLES  INSUFFICIENT  -STRUC 
1TURE  NOT  IN  EQUILIBRIUM-') 

252  CONTINUE 
CALL  SWAY 

253  CONTINUE 
RETURN 
ENO 
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SUBROUTINE  SYMPAK 

COMMUN  KJ (20, 20)  ,KP( 20,20  I ,KS( 20,20) ,KX (4 ,20,20  I , KYI  3,20,201  ,KV( 20 
l *20) *K(  20*20)  ,  NS,  NX, NY, l CX , I CY , XO ( 20 ) ,Y0(20) , KXRO T w < 4 , 20 , 20 ) ,KVROT 
2C( 3,20, 20) *K YROTWI  3 , 20 *20 ) , H< 20 ) ,C I ( 20 ) * W l ( 20 ) , WBL ( 20 ) , W( 20 )  ,F BL ( 2 
30)  ,WBI  (20)  *E*EW,EGC,EGW,FX(20) , F Y ( 20 )  *  X  X ( 4 ) ,YY (3) ,CX(20),CY(20) ,RX 
4(4,20)  ,UX( 4*20,20 ,KDX(4, 20*20) ,TTX(20,20) ,ROT (20) ,VX(4,20)  *  T VX ( 20 
5)  ,RY( 3 ,201  ,0Y( 3 ,20 ,20) ,KDY(3,20  ,20 ) ,TTY( 20,20 ) , VY( 3,20) ,TVY(  20 ) ,MT 
6( 20) , TVXAI 20  )  ,T VYA( 20) , XH ( 20 ) , YH ( 20 ) , FB I ( 20 ) , KXROTC ( 4 f 20 , 20  I  , I  NO  I C 
T  X  (  4  )  i INOICYt  3 ) , HX ( 4 , 20 ) ,CIX(4,20) ,  W  I  X ( 4 , 20 ) *WBLX(4»20) , *IX  (  4 , 20  ) » F  8 
8L  X  (  4 , 20 )  •  W  B I  X  (  4,20)  ,FBIX(4,20)  ,HY(  3,20)  ,  C  I  Y  (  3,20)  »WlY  (  3,20)  »WBLY(  3 

9.20)  ,WY( 3,20) ,FBLY( 3,20) ,WBIY( 3,20) *  FBI Y( 3,20 ) 

COMMON  KA(20,20),KB(20,20),KU(20,20),KE(20,20),KCI20,20),KF(20,20> 
1 ,KG( 20,20) ,KH(20*20) ,KI (20,20)  ,KK(2C, 20  I ,KL( 20,20 )  ,KM( 20,20) ,KN(20 

2.20)  ,K0(  20,20)  ,KQ(  20, 20)  ,  KRI2C',  20  J  ,KT(20,20  )  ,KW(  20,20)  ,  TX(4,  20, 20) 
3,TY( 3,20,20) ,RT( 20,20) ,SX (4,2C ) ,SY( 3, 20) , RUT XC l 4 , 20 > , ROTX W( 4 , 20 ) , 
4R0  T YC ( 3,20) ,ROTYW( 3 , 2  0  >,  T  TVXC  (  20 , 20  ),  T  VXC  (  20 , 2C  ) , TTVYC ( 20 ,20 » ,  TVYC 
5(20,20) ,CMXB( 4,20 ) ,CMXT( 4 ,20 ) , WMXB ( 4 , 20 ) , WMXT ( 4 , 20 ) , FBM X ( 4 , 20 ) , Wb M 
6XC<4,20 I , WBMXW( 4,20)  ,CMYB( 3,20)  ,CMYT ( 3 ,20 ) , WMYB (3,2C ) ,WMYT (3,20) ,F 
7BMY( 3,20) , WBMYCI 3,20 ) ,WBMYW( 3,2O),CSX(4,20),^SX(4,20),FBSX(4,20),W 
8BSX( 4.20 ) ,CSY( 3,20) ,WSY( 3 ,20 , FBSY13,20 ) , WBSY( 3,20) ,CFX(4,20 ) , WFX( 

94.20)  ,FBFX(4,2C) ,WBFX(4,20) ,CFY(3,20) ,WFY ( 3,20) ,FBFY(  3,  20) 

COMMON  WBFY(3,20),DEX(4,20),DfcY(3,20),lNSYM,ITERX,ITERY,IX,iY 
REAL  KJ  »KP*KS*KX,KY,KV,K»KXROTC*KXROTW*KYROTC,KYKOTto*KDX*KDY,MT 
REAL  KA,KB ,KD , KE , KC , KF , KG , KH , K I , KK , KL , KM , KN , KO , KU , KR , KT , KW 

DO  100  1=1, NS 
DO  101  J= 1 »  NS 
TTX (  I , J  1=0.0 
DO  102  L= 1 , NX 
ADD=KX ( L, I , J ) 

TTX ( I, J)=TTX( I , J ) ♦ ADD 

102  CONTINUE 
101  CONTINUE 
100  CONTINUE 

DO  103  1=1, NS 
DO  104  J=1 , NS 
TTY ( I , J ) =0 .0 
DO  105  L= 1 , NY 
ADD=KY I L , I , J) 

TTY  (  I  ,  J)=TTY(  I  ,  J  M-AUD 

105  CONTINUE 
104  CONTINUE 

103  CONTINUE 

CALL  TOKYO ( NS , TT  X ) 

106  CONTINUE 

DO  107  1=1, NS 
SX ( 1,1 )=0.0 
DO  108  J=1 »NS 
PROD=TTX( I,J)*(-FX( J ) ) 

SX ( 1 , I )=SX( 1, I ) +PROD 

108  CONTINUE 

107  CONI INUE 

DO  109  I=£,NX 
DO  110  J=1 ,NS 
SX( I , J ) =SX ( 1 , J ) 

110  CONTINUE 

109  CONTINUE 

CALL  TOKYO ( NS , TTY ) 

112  CONTINUE 

DO  113  1=1, NS 

SY( 1, I )=0.0 

DO  114  J=1 , NS 

PROD  =  TTY ( I  ,  J ) *(-FY ( J  )  ( 

S Y  (  1,  I  )  =  SY(  1,  I  )  4-PROD 
114  CONTINUE 

113  CONTINUE 
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DO  115  I =2, NY 

on  H6  j  =  i  *  ns 

SY(  I  , J ) =SY( 1, J) 

116  CONTINUE 

115  CONTINUE 

117  CONTINUE 
RETURN 
END 

SUBROUTINE  WRITL1 

COMMON  KJ ( 20, 20  I ,KP( 2  0,20 ) ,KS( 20, 20) ,KX( 4, 20 ,2") ,KY(  3,20,20)  ,KV(2C 

1 .20)  »  K  ( 20,201 , NS, NX, NY ,1 C  X  » ICY , X  0 ( 2  0 )  , Y  0( 20 ) *KXR0TW(4, 20,20)  ,KYRUT 
2CI 3,20,20) ,KYROTW( 3,20,20 ) ,H (20 J ,CI (20) ,Wl (20) ,WBL ( 20) , W(20) ,FBL ( 2 
3O),W8I(2O),E,Eto,EGC,EGW,FX(20),FY(2r')»XX(4l*YY(3),CX(2Q)»CY(2O)*RX 
4(4,20) ,0X( 4,20,20) ,KDX( 4, 20,20 )  ,TTX( 2^ ,20) , ROT ( 20) ,VX( 4 ,20) ,TVX( 20 
5)  ,RY(  3,20)  ,()Y(  3,20,20  )  ,KDY(  3,20,2C  )  ,TTY(  20 ,20)  ,  VY(  3,20)  ,TVY(  20)  ,MT 
6(20) , TVXA( 20 ) ,TVYA( 20) ,XH(2C» , YH ( 20 > , F B I ( 2 0 ) , KXRUTC ( 4 , 20 , 20 )  ,  l NDI C 
7X(4) , INDICY(3 )  ,UX (4  ,20) ,CI X ( 4 , 2n ) *  W I  X ( 4 , 20 ) ,WBLX(4,20) , MX ( 4 , 2C ) ,FB 
8l_X(4,2C  >,  WB1XI  4,21  ),  FBI  X(  4,20  )  ,HY(  3,20)  ,CIY(3, 20)  ,W  IY(  3,20)  ,  WBLY(  3 

9.20)  ,WY(  3,2.1)  ,FBLY(  3,20)  ,WBIY<  3,2U  ,F6IY(  3,20) 

COMMON  KA(20,2C)  ,  KB  (20,20  ),KD(  20,20  )  ,KE(  20, 20,  KC(  20,20  ),KF(  20, 20) 
1 ,KG( 20 ,29 )  , KH(20 ,20 )  ,Kl ( 20 ,20 ) ,KK  <  2n, 2C  > ,KL( 20,20) ,KM( 20, 20) ,KN( 20 

2.20 )  ,K0( 20,20) ,KQ(2l  ,20  , KR(2C ,20) ,KT (20,20) ,KW(  20,20) , TX(4, 20*  20) 
3»TY(  3,20,20) *RT(2 3, 20)  ,SX(4,20),SY(3,2C)»RUTXC(4,20)»ROTXW(4,20)» 
4R0TYC ( 3,20) , ROT YW( 3,20) , TTVXCf 20 ,20 ) , TVXC ( 20 ,20  ) , T T V YC ( 20 , 20 ) ,  TVYC 
5(20,20) ,CMXB( 4,20 ,CMXT( 4,20 ,WMXB< 4,20) ,WMXT( 4,20) ,FBMX( 4*20) ,WBM 
6XC (4,20) , WBMX W ( 4,20) ,CMYB( 3,2r ) ,CMYT ( 3,20) , WMYB<3,20) ,WMYT(3,20),F 
7BMY(  3,20  ),  WBMYC!  3,20  ,  WBMYW(  3,20  ,CSX(  4,20  ,WSX(  4,20)  ,FBSX<  4,20)  ,  W 
BBS X( 4, 20) ,CSY( 3,20) ,WSY( 3,20) ,FBSY( 3,20 ) , WBSY( 3,20) ,CFX( 4,20) ,  WFX( 

94.20)  ,F8FX( 4, 2C )  ,WBFX (4,20) ,CFY( 3,20) ,WFY( 3,201 *FBFY( 3, 20) 

COMMON  W3FY ( 3,20 ,DEX( 4,20) ,DEY( 3,20) , INS YM, 1 TfcRX, ITERY, I  X, I Y 
REAL  KJ ,KP ,KS ,KX,KY, KV ,K , KXRUTC , KXR OT W , KYRU TC , K YRUT w , KDX, KDY , MT 
REAL  KA ,KB ,KO,Kfc ,KC , KF ,KG  ,KH ,KI  , KK , KL , KM , KN , KU , KU , KK , KT , KW 
WRITE(6,101)NS,NX,NY 

1C  1  FORMAT ( 1H1 ,2X,* NUMBER  OF  STORIES  =  • , I  2 , 4X , • NUMB ER  UF  X  FRAMES  =  * 
1 , 12, 4X, 'NUMBER  OF  Y  FRAMES  =  ’,12) 

WR I TE( 6, 1 16) I TERX, 1 TERY 

116  FORMAT ( 1H0// *2X, 'NUMBER  OF  ITERATIONS  IN  X  DIRECTION  LIMITED  TO  ', 
1 13, ?X, 'NUMBER  OF  ITERATIONS  IN  Y  DIRECTION  LIMITED  TO  ',131 

WRI T  E ( 6, 102 ) 

K2  FORMAT!  1H0//,  TIG, 'FLOOR  NUMBER  '  ,T  32  ,  •  APRL  I  ED  LOAD  (  K  I  P  )  •  ,  T6  1 ,  '  PU  1  NT 
1  OF  APPL  ICAT ION(  IN) ', T93,  'CENTROID  OF  ST  I FFNE SS l I N )  •  ) 
l*IRITE(6,103) 

1T3  FORM  AT ( T28 , • X  DI REC T I  ON'  , T43 , • Y  D I R EC T I  ON '  , T 60 , ' X  D I R EC T I  ON ' , T 75 ,  • 
1Y  DIRECTION* , T92 ,' X  D I RECT I  ON '  , T1C7 , ' Y  DIRECTION') 

DO  104  1=1, NS 

WRITE (6, 105)1, FX(I),FY(l),XH(I),YH(I),CX(I),CY(I) 

105  FORMAT(1HO,T16,I2,T2B,E8.2,T43,F8.2,T60,F8.2,T75,FB.2,T92,F8.2, 

IT  1C7 , F8 . 2 ) 

104  CONTINUE 

DO  1C6  1=1, NX 
WR I T E ( 6,107)  I 

107  FORMAT ( 1H1 ,2X, 'FRAME  X  NUMBER ', 1 10 , 2 X ,* MEMBER  PROPERTIES') 

WR I TE ( 6  » 1C  B ) 

108  FORMAT ( 1HC/,T8, • STORY  NUMBER' , T22 ,' STORY  HE IGHT ', T 3b, 'COLUMN  INERT 
1IA' , T52  » ' WALL  INERTIA'  ,T65,  'WALL  W I DTH ' , T 8 1 , ' BE  AM  SPAN (  IN  )  • , T 102 , 
2'BEAM  INERTIA! IN4) •  ) 

WR I TE ( 6 , 109) 

109  FORMAT ( T26 . ' ( IN) ' ,T40, *( IN4) ' ,T55,' (IN4) ' ,T68, ' ( IN) ' ,T77,'FRAME-BE 
1AM* ,T 89, 'WALL-BEAM' , T 100 , 'FRAME-BEAM* ,T1 12, 'WALL-BEAM* ) 

DO  110  J=1 , NS 

WR I TE ( 6 , 1 1 1 ) J , HX ( I ,J) , Cl  X ( I , J ) ,wIX( I , J ) ,WX(  I , J ) ,FBLX(  I , J ) ,W3LX ( I , J 

1) , FBIX( I , J ),WB1X( I , J ) 

111  FORMAT ( 1H0 ,T14,I3,T22,F8.2,T34,F12.2»T50,F12.2*T63»F8.2»T75*F10.2* 
1T85,F1G.2,T98,E10.2,T1°8,F10.2) 

110  CONTINUE 
WRITE(6»112)XX< I ) 
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112  FORMAT  (  1(6//  »  2  X »  ’THAME.  LOCATION  PROM  kE  PE  kE  NO  I-  PUlNT  =  '  ,  i  a  .  2  ) 

1  ri,  CONT  I  NO' 

no  m  i  - 1  *  n v 

Wk 1 Tt  ( ft, 1 1 4 ) 1 

ll4  FORM  AT  (  10  1  ,  -’X  ,  '  I  RAMt  Y  NUMBfc  R  '  ,  l  1  n  ,  2  X  ,  '  ME  MQt  R  PROPt  kT  IPS'  ) 

WR I T  E ( ft  ,  1  f  H  ) 

HR  I  T  E  (  ft  ,  l  r-  9  ) 

no  1 1  ;>  j.  i  ,ns 

HR  I  Tt  (  ft  ,  1  1  1  )  .1  *  OY  (  I  ,  J I  , Cl  Y  (  1  ,  J  )  ,  to  I  Y  (  l  ,  J  )  ,  *  Y  (  I  ,  J  )  ,  POL  Y  (  I  ,  J  )  ,  „HLY(  I.J 
l  )  ,  FB  I  Y (  1  ,  J  I  ,  WI31  Y  <  I  , J I 
1  1 5  CON  f INOP 

HR ITfc (6, 11 2) YY( I J 
1  1  3  CONT  1  NIJF 
RETURN 
END 

SUBROUTINE  HXUDX 

COMMON  Kj(20»2?)»KP(20»2’''),K3(20»20)»KX(4,2n»20),KY(3,20»2C)*KV(2C 

1  ,26  )  , K( 20 ,  2i  I ,NS ,  NX , NY , I CX  ,  I CY , XO ( 20  I  , Y0<  20 ) , KxRUTW ( 4 , 20*  20  I . KYRuT 
2C<  3, 20, 2C ) .KYRUT  to ( 3 , 2  0 ,2  o ) ,M ( 29  ) »C I l 2<  > , to  I  (20 ) , HBL ( 2C ) , to ( 20  I , F BL ( 2 
30 ) » W  B  I (  2? ) »E,EW»EGC,E.>w»FX(2C) » I YI20) , XX l 4 ) ,YY( 3) ,CX(20),CY(20) »RX 
4(4,20)  ,DX(  4,26,20)  ,KUX  U,20, 7C  )  ,TTX(  20  ,20  )  ,ROT  (  20  »VX<  4 ,2G)  ,  TVX  (20 
5)  ,R Y( 3, 20  I  ,uY( 3,2) ,  20  )  ,  K  D  Y  (  3 , 2  9 , 20  )  ,  T  T  Y  <  2C  ,  20  )  ,  V  Y  (  3 , 20  )  ,  TV  Y  (  20  )  ,  MT 
6( 20) , TVXA(20 ) ,TVYA(2C  )  ,XH( 2" I  , YH( 20) , FBI ( 20 ) ,KXKOTC (4 ,20, 2Q)  , INUlC 
7X( 4) , I NOlC Y( 3 ) ,HX  t 4  ,20  > , C I X( 4 , 2'  )  , WI X( 4,2( )  ,HBL  X( 4, 2C ) , WX (4 , 20 ) »FB 
8  L  X  (  4  »  2  0  )  ,WB I  X ( 4,2T)  , F  B I  X (  ,,20)  ,HY( 3,20) ,C1Y (3,20)  ,wlY( 3,20) ,WbLY( 3 

9.20)  , wY ( 3,20)  ,EHlY( 3 ,20)  , to 31 Y ( 3,2  0) , FBI Y (  3,20) 

COMMON  KA1 20,20 )  ,KU (20,7? ) ,KD( 2  ,2C I ,KE ( 20 ,20 ) ,KC( 2F  ,20  )  ,KF( 20,20 ) 

1 , KG i 20 ,20) ,KH(2r ,20 ) ,K I ( 20,20) ,KK( 20, 2C ) ,KL (20, 20 ) ,KM( 2 0. 20) ,KN (20 
2, 20  )  ,KO(20,2O)  ,KQ(21’  ,  ?  ^  ,KR(  2C  ,20  )  ,KT  (20,20)  ,KW  (  20,20  )  ,  TX(4,2C,  201 
3  ,  T  Y  (  3,20,20  ,  RT  (  2  ) , .  G )  ,  S  X  (  4  *  2C  I  ,  SY  (  3,20)  ,  RUT  XC  (  4 , 20  )  ,  RG  T  X  to  (  4 , 20  )  , 
4R0TYC  (  3 , 21' )  ,  RUT  Y  Vi  (  3 , 2  0  )  ,  T  T  VXC  (  20 , 20  )  ,  T  VXC  (  20 , 2C  )  ,  TT  V  YC  (  2C  ,  2  0  )  ,  T  V  YC 
5(20,20  j  ,  CM  X&(  4,20  )  ,CMXT  (  4,20  )  ,WMXB(4,2‘' )  ,  WMXT  (  4 , 20  )  ,  E  BMX  (  4, 2  0  )  ,  WtJM 
6XC  (4,20)  ,  WBMXW  (4,20)  ,CMYti(  3,20)  ,C  MY  T(  3,20  )  ,WMYB  (  3,2C  )  ,to  MY  T(  3,20)  ,E 
7 13  M  Y  <  3,20)  ,  WBMYC(  3,2'  I  ,  toHMYWl  3,2r  )  ,C  3X<  ^  ,20  ,  WSX(4,20)  ,FrtiX(  '♦,20)  ,  W 
8BSX(4»20  ,CSY(3,20)  ,toSY(  3»2n)  ,FliSY(  3,20)  ,VilUSY(3,20)  ,CEX14,2C  )  »  WE  X  ( 

94. 20)  , FBFX(4, 20 ) , toBFXt 4,20) ,CFY( 3,20) ,WFY ( 3,20) ,FBFY<  3, 20 ) 

COMMON  WBEY(3,20)  ,l)EX(4,??),DFY'  ’•  'S,!NSYM,ITFRX,ITEKY,IX,IY 
REAL  K .  I  ,  K  P  ,r.j  ,KA  ,AY,KV,K,  KXROT  C  ,  KXKOT  W ,  KY  RU  TC  ,  K  YROT  to  ,  KUX  ,  KO  Y  ,  MT 
RIAL  KA  ,KB ,KD,KE , KC , KF ,KG  ,KH ,KI ,KK, KL , KM , KN , KO , Kg , Kk , KT ,KW 

OH  5 "4  1=1, NX 
00  505  J=1 , NS 
RX ( I , J ) =X0( J ) -XX ( I ) 

5^5  CONT I NUF 
504  CONTINUE 

DO  50>6  1=1, NX 
00  507  J=1 , NS 
DO  508  L=1 , NS 
OX ( I , J,L)=o.n 

508  CONTINUE 
5 '7  CONTINUE 
506  CONTINUE 

00  509  l  =JL  ,  NX 
DO  510  J= 1 , NS 
DX (  I  ,  J  ,  J  )  =  RX (  I , J ) 

IF<  J.CQ.DGO  TO  510 

DX(  I.J, (J  —  1) > =-RX(  I  ,  ( J-l  )  ) 

516  CONTINUE 

509  CONTINUE 
RETURN 
END 

SUBROUTINE  DXOTTX 

COMMON  KJ (20, 20) , KP 1 20 ,20 ) ,KS(20»20)»KX(4,2Q,2G),KY(3,20»20) » KV ( 20 

1 .20)  ,Kl 20,20 ) ,NS , NX, NY, I CX, ICY, X0( 20 ) ,Y0( 20 ) , KXRU T W ( 4 , 20 , 20 ) .KYRUT 

2  C ( 3,20,2°) » K  YROTW ( 3,20,20) *H(2G) , C I ( 20 ) ,WI (20) ,WBL(2C ) , to (20) ,FBL( 2 


‘ 

« 

' 

. 


c*  '  1  Ok)  l  .  V  '  ! 


3(')  »W8!  (20)  »E*EW*fcGC  *EGW»FX(20)  »  F  Y  (  2  0  )  *  X  X  {  A  I  ,VY  (  3 )  ,  CX ( 2  0 )  ,  C  Y  (  20  )  ,RX 
4Uf20)  ♦0X(4,2O,?r  )  , KDX (4, 2 0,20)  ,TTX(2C  ,20)  ,  KUT  (  20)  ,  VX(4,20»  ,TVX(20 

5),RY(3,20),[JY(3,20*20),KDY(3,20»20),TTY(2C,20),\/Y(3»20),TVY(20),M1 

6  (  2n)  ,  T  VXA<  20  I  *  T  V  Y  A  (  2  C  )  ,  XH  (  20)  ,  YH  (  2C  )  ,  F  B  l  (  2  0  )  ,  K  XROTC  (  A  ,  2C ,  20  )  ,  I  NUIC 

7X(4),INOICY(3),HXI4,2Q),ClX(4,20),WlX(4,2C),WBLX(4,20),WX(4,20l,FB 

BLX{4,20),WBIX(4,2D),FHIX(4,20),HY(3,2C),CIY(3,2C),WIY(1,20),WBLV(3 

9.20)  ,  WY  (  3,20)  ,  FBL Y I  3,  20)  ,  *BI  Y  <  3,201  ,  F  B  I  Y  (  3,20) 

COMMON  KAl 20 ,2r ) , KB<  20,20 ) ,KD(20,2C ) ,KE (2C ,20 ) ,KC ( 20,20  >  *KF ( 20,20 ) 
l»KC ( 20, 20>,KH( 20, 20) ,KI (20,20)  ,KK( 2 0,20 , XL (20,20 >, KM (20,20) ,KN (20 

2.20)  , KOI  20,20) ,KQ( 2'  ,20)  , KR(2C,2n ) ,KT ( 20,20  )  , KW ( 20 , 20 ) , TX(4*20»20 ) 
3,TY(  3,  2  0,20) ,  R  T  (  2  0  »  2  C  )  ,  SX  (4 ,20  ,  SY  (  3 , 20 ) , RUT XC ( 4 , 20 ) , ROT  X M ( 4 , 20 ) , 
4R0TYC(  3,20  )  ,ROTYW(  3,2r'  )  ,  T  I  VXC  (  20 , 20  ),  T  V  XC  (  20 , 20  ),  T  T  VYC  (  20 , 20  ),  TVYC 
5(20,20)  , C  M  X  B ( 4 , 2  C )  ,  CM  XT ( 4 ,20 )  ,  WMXB ( 4, 2C )  ,  WMXT ( 4 , 20 )  ,FBMX(4,20) ,  wBM 
6XC ( 4 , 20 ) , W  BM  X  W ( 4 ,20 ) , CMYB ( 3 , 2C  )  ,CMYT(  3,20) *  WMYB ( 3 , 2C ) ,WMYT(3,20),F 
7BMY( 3,20) , WBMYC ( 3,2C ) , WB  M Y to( 3  »  20 )  ,  C  SX ( 4 , 2C )  , WS  X ( 4 , 20 ) , F  BSX ( 4 , 20 ) ,w 
8BSX(4,20) ,CSY( 3,20)  ,WSY( 3,20) ,FBSY( 3,20) ,  wB  SY ( 3,20) ,CFX(4,20) ,MFX( 

94.20)  ,  F  BF  X  (  4 , 20  ,  WBF  X  (  4 , 2  0  )  ,  CF  Y  (  3 , 20  )  ,  WF  Y  (  3,20  ,FBF  Y(  3,  20  I 
COMMON  WBFY  (3,20)  ,DLX(4*2G)  ,L)EY(  3 ,20  )  ,  INSYM,  I  TERX,  l  TERY,  l  X*  I  Y 
REAL  KJ,KP ,KS ,KX,KY,KV,K, K XROTC  ,  KXR OT W , KY KO TC , K YKOT W , KU X , KUY , MT 
REAL  KA ,KB ,KD, KE  ,KC , KF ,KG  ,KH , K I  , KK , KL ,KM, KN,KO ,KU,KR , KT ,  KW 

DO  512  1=1, NX 

DO  513  J=1 ,NS 

00  514  M=  1  ,  NS 

SUM=O.P 

DO  515  L= 1 , NS 

PROD=KX( I ,J,L)*DX( I ,L,M) 

515  SUM=SUM+PRCD 
KDX ( I , J,M)=SUM 

514  CONTINUE 
513  CONTINUE 
512  CONTINUE 

DO  516  1=1, NX 
DO  517  J=1 ,NS 
DO  518  L= 1 , NS 
T  X (  I , J,L)  =  XX(  I  ) *  KDX (  l , J, L  ) 

5 1 B  CONTINUE 
517  CONTINUE 

516  CONTINUE 

DO  519  1=1, NS 
DO  5 2D  J=1 ,NS 
T  T  X (  I , J ) =0.0 
DO  521  L=1 , NX 
A  D  D  =  T  X ( L,  I  ,  J) 

TT  X ( I, J)=TTX(  I , J ) + ADD 
521  CONTINUE 
52F  CONTINUE 
519  CONTINUE 
RETURN 
END 

SUBROUTINE  VXOTVX 

COMMON  KJ(2C,2C)  , KP (20,20 ) ,KS ( 20,20) ,KX(4,2C,20 ) , KYI  3, 20, 20) ,KV(20 
1,20  ,  K(  20 ,20  )  ,NS,NX*NY,1CX,1CY»XU(20  ,Y0(20),KXR0TW(4,20,20)  ,  K  YROT 
2C(  3,20,20)  .KYRUTM  3,20,20)  ,H(  20)  ,C  1(20  ,Wl  (20)  ,WBL(  20  ,W(  20)  ,FBL(  2 
30)*WbI(2P)»E»EW»EGC»EGW»FX(20),FY(20)»XX(4),YY(3)»CX(20),CY(2Q),RX 
4(4,20  ,DX(  4,20,20  ,  KDX  (4, 20,20  ,TTX(  2C  ,20  ),  ROT  (20  ,VX(  4,20)  ,TVX(  20 
5  ) ,RY( 3,20)  ,DY(  3 ,23, 20  )  ,KOY( 3, 20 ,2r )  ,TTY( 20,20 ) , VY( 3,20) , TVY( 20) ,MT 
6(20),TVXA(20),TVYA(2C),XH(2C),YH(20),FBI(20) , KXRUTC ( 4 , 2 0, 20)  , I  NO 1C 
7X(4)  ,  INDICYI  3  )  »HX  (4,20)  ,CI  X(  4(,  2'  )  ,WIX(  4,2f  )  *  WBL  X l  4,20)  ,  WX ( 4 , 20  )  ,FB 
8LX(4,20)  ,WBIX( 4,23  I ,FBIX( 4,20) ,HY( 3,20 ) ,CI Y( 3,20) ,WI Y(3,20) , WBLY( 3 

9.20)  , W Y ( 3,20) ,FBLY( 3, 20 , rtBI Y( 3,20)  ,FBIY( 3,20  ) 

COMMON  KA( 20,20) , KB (2 0,20) ,KD( 20,20) ,KL( 2 0,20) , KC ( 20,20) ,KF ( 20,20 ) 
l  ,KG(  2  0,20)  ,KM.(  20,20  )  ,K  I  (  20,20  )  ,  KK  (  2%  20  )  ,KL  (  20, 20)  ,KM(  20,20)  ,KN<  20 

2.20 ) , KOI  20,20) ,KQ( 20,20 ) ,KR( 20 ,20 ) ,KT(2C, 20 ) ,KH (20,20  I , TX( 4, 20, 20 ) 
3, TY( 3,20,20) ,RT( 20,20) ,SX (4, 2C ) ,SY( 3, 20 ) ,ROTXC( 4,20) ,ROTXW(  4,20) , 
4R0TYC( 3,20) ,ROTYW(3,20) *T TVXC( 20 , 2C  >  »  T  VXC  <  20 , 20  »  TTVYC  (  20 , 20  )  ,  T  VYC 
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5(  20*20)  »CMXB(  4,201  ,CMXT (  4  ,20  )  ,WMXH(  4,20)  ,  WMXT  (  4,20 )  ,FBMX(  4,20)  »  WBM 
6  X  C  (  4  *  2  0  ) * WBMXW(4 , 20 ) , CMYB ( 3  *  20 ) , CMY T(  3 , 20 ) , WMYB ( 3 , 20 ) , WMYT  (  3 ♦ 20  I , F 
7BMY(  3,20  )  , W»MYC(  3 , 20  I  ,  WBMYW  (  3  *  20  )  ,  CSX  (  4 , 20  I  ,  WS  X  (  A  ,  20  )  ,  F  fiS  X(  A  ,  20  )  ,  M 
8B  SX(4, 20 ),  CSYI3, 23), WSY  (3, 20  ), F  I3SY(3,20)»WBSY  (3, 20  ),CFX(4, 20  ),WFX( 
94 ♦ 20) , FBFXI A, 20 ) ,WBFX(4,20) ,CFY( 3,20) , WFYI 3,20) ,FBFY(  3, 20) 

COMMON  WBFY( 3,20) ,DEX(4,20) ,DEY( 3,20) , INSYM, I TERX, ITERY , I  X, l Y 

MEAL  KJ ,KP,KS ,KX ,KY,KV,K , KXROTC , KXROT W , KY ROTC , K YRUT W, KDX , KOY , MT 

RE  AL  KA  *KB ,KD,KE ,KC ,KF,KG ,KH,KI , KK,KL ,  KM  ,  KN  ,  KO  ,  K(J ,  KR  ,  KT  ,  KW 

DO  522  I  =  1 »  NX 

DO  523  J= 1 , NS 

VX( I , J)=0.0 

DO  524  L= 1 , NS 

PR0D  =  K(3X  (  I  ,J,L  )*ROT(  L  ) 

vx ( i ,  j  )  =vx (  i ,  jh-prod 

52A  CONTINUE 
523  CONTINUE 
522  CONTINUE 

DO  525  J=1 ,NS 

TVXI J) =0.0 

DO  526  1=1, NX 

TVXI J)  =  TVX<  J[+VX  I  I , J ) 

526  CONTINUE 
525  CONTINUE 
RETURN 
END 

SUBROUTINE  HYUDY 

COMMON  KJI 20,20)  ,KP I  2  0 ,20  ) ,KS( 20, 20) ,KX I  4, 20, 20) ,KY<  3,20, 20) ,KV( 20 

1.20)  ,K( 20,20) ,NS,NX ,MY, I CX , I C Y , XO ( 20 ) ,YU( 20) , KXROTW ( 4 , 20, 20 ) , KYROT 
2C( 3,20,20) .KYRUTWI 3, 20,20) ,HI20) ,C I ( 20) ,WI (20) ,WBL ( 20) ,W( 20) ,FBL( 2 
30)  ,WBI  (20)  ,E,EW,EGC,ECiW,FX(20)  ,  F  Y  (  20 )  ,XX(4)  ,YY(  3)  ,CX(20),CY(20)  ,RX 
4(4,20)  ,0X14,20,20)  »KDX(4,20,2r!)  ,TTX(2n,20)  ,  ROT  (  20  )  ,  VX  (  4 , 20  >  ,  TVXI  20 
51 ,RY( 3,20) ,DY( 3,20, 20 ) ,KDY( 3,20 ,20) ,TTY( 2C,20) , VY( 3,20) ,TVY( 20) ,MT 
6(20) ,TVXA(20)  ,TVYA( 20)  ,XH(20) , YHC2C) , FBH  20) ,KXROTC<4,20,20) , INOIC 
7X(4),INDICY(3),HX(4,20),ClX(4,20),WlX(4,2f'),WBLX(4,20t,WX(4,20),Frt 
8LX(4»20)  ,WBIX(4,20)  ,FBIX(  4,20  )  ,HY(  3,20)  ,CIY(3,20)  ,WlY  (  3,20)  ,  WBLY(  3 

9.20)  ,WY( 3,20) ,FBLY( 3,20) ,WBIY(3,20) , FB I Y ( 3,20) 

COMMON  KA<  20,20) , KB ( 20, 20 ) ,K0( 20,20 ) ,KE (20,20) ,KC( 20,20  )  ,KF( 20,20) 
1,KG(20»20) ,KH( 20,20 ) ,KI (20,20) ,KM 20,20 ) ,KL( 20,20 ) , KM (20,20 ) ,KN(2C 

2.20)  , KOI  20,20) ,KU( 2C, 20) ,KR( 20, 20) ,KT( 20,20 ) ,KW( 20, 2C ) ,TX( 4, 20,20 ) 
3,TY(  3, 20,20) ,RT( 20,20)  ,SX (4,20) ,SY( 3,20  > ,R0TXC(4 ,20) ,ROTXW( 4,20) , 
4R0TYC(  3,20  )  ,ROTYWl  3,20)  ,  T  IVXC(  20 , 20  )  ,  TVXC  (  20 , 20  )  ,  T  TVYC  (  20 , 20  )  ,  T  VY I. 
5(20,20)  ,CMXB( 4, 2G) , CM XT (4,20) , WMXB( 4, 20)  ,  WMXT (4,20)  ,FBMX( 4,20) , WBM 
6XC ( 4 , 20 )  ,  WBMXW( 4,20 ) , CMYB ( 3 , 20 ) ,C MY T ( 3 , 20 ) , WMYB ( 3 , 20 ) , WMYT ( 3 , 20 ) , F 
7BMY(  3,20!  ,Wt)MYC(  3,2  0)  ,WBMYW(  3,20)  ,CSX(  4,20)  ,WSX  (4,20  )  ,FBSX(  4,20)  ,W 
8BSX(4,20) ,CSY( 3,20) ,WSY( 3,20) ,FBSY( 3,20 ) ,W6SY( 3,20) ,CFX(4»29) ,WFX( 

94.20)  ,FBFX( 4,20) ,WBFX (4,20 J ,CFY( 3,20 ,WFY( 3,20) ,FBFY( 3, 20) 

COMMON  WBFY(3,20),DEX(4,20),DEY(3,20),INSYM,ITERX,ITERY,IX,IY 
REAL  K J, KP, KS, KX, K Y, K V, K, KXROTC, KXRUTW.KY ROTC, KYROTW, KDX, KDY,MT 
REAL  KA ,KB,KD,KE,KC ,KF,KG,KH,KI ,KK,KL , KM , KN , KO , KU , KR , KT  ,KW 

DO  604  1=1, NY 
DO  605  J=1,NS 
RY(  I  ,  J  )  =  Y«0  (  J  )  -  Y  Y  (  I  ) 

605  CONTINUE 
604  CONTINUE 

DO  606  1=1, NY 
DO  607  J=1 , NS 
DO  608  L= 1 , NS 
D  Y  (  I  ,  J  ,  L  )  =  0 . 0 
608  CONTINUE 
607  CONTINUE 

606  CONTINUE 

DO  609  1=1, NY 
DO  610  J= 1 »NS 
DY (  I  ,  J, J)  =  RY( I , J ) 

IFIJ.EQ.1IG0  TO  610 
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610  CONT 1  NUF 
6r9  CONT  I  NUE 
RETURN 
F.ND 

SUBROUTINE  DYUTTY 

COMMON  KJ(20,?~) ,KP( 20,20 > ,KS( 2' ,201 »KX(4,2r ,20 ,KY( 3,20,20)  ,  KV(2< 

1 ,2r)  , K( 20 *20  »  ,NS ,NX, NY, 1 CX, I CY, XU( 2n )  , YU( 20) , KXROTWl  4 ,2C» 20)  ,KYRUT 
2C(  3  *20 ,20  I  ,KYROTM  3,  20  *20  >  ,H(  20  (  ,CI  (20  )  ,  Wl  (  2C  >  ,WBL(  2C  )  ,  R(  20)  ,FBL  (  2 
30)  ,WBI  (  20) »E*FW,EGCfE6W,FX(20) *  T  Y ( 2  0 )  *  X  X  (  4  )  , YY ( 3) ,  C  X  ( 20  )  , C Y ( 20 )  » RX 
4(4 ,2  3)  , DX<  4  *20,20)  ,KDX (4, 2  0,20  |  ,TTX( 2C ,20 )  ,RUT< 20) , VX(4 ,20 , TVX( 20 
5)  ,RY(  3,20)  ,DY(  3,2?  ,20)  ,K0Y(3,2'',2<~  )  ,TTY<  20,2C  )  ,  VY(  3,20)  ,TVY<  20)  ,MT 
6(  20)  ,  TVXA(  2'-1 )  ,TVYA(  20  )  ,XII(20)  ,  YH<  20  )  ,  FBI  (  20)  ,KXROTC  (4,20,20)  ,  INUIC 
7X(4)flNDICY(3)*HX(4,2C),ClX(4,20),wlX(4,2C),WhLX(4,2C),WX(4,2C),Fh 
8LX(4,20)  ,  WB1X(4,2?  )  »FBU(  4,20)  ,HY(3,20  )  ,CI  Y  (  3,20)  ,Wl  Y  (  3,20)  ,  *BLY(  3 
9,2C),WY(3*20),FbLY(3,20)»WBlY(3*2r),FbIY(3,20> 

COMMON  KA( 20, 20) ,KB( 20,20 ) ,KD( 20 ,20) ,KE(2G, 20) ,KC( 20,20) ,KF(2C,2G ) 

1 , KG( 20,20) ,KH( 20,20) ,KI ( 2C, 20) ,KK( 20, 20 ) ,KL( 20,20 ) ,KM( 20,20) ,KN( 20 
2,20  )  ,KG(  20,20  )  ,KCl(2C  ,20)  ,  KR(2C,  20  )  ,KT(  20,20  ),  KW  (  20 , 20  ),  TX  (  4, 20 , 20  ) 
3, TY( 3,2  0,2  0) , R  T ( 2  0,20  ,SX (4,20 ) ,SY( 3, 20) , ROTXC< 4,20) ,ROTXW( 4,20) , 
4R0TYC( 3,20) » ROT YW ( 3,20) , T  TVXC  < 20  ,20) ,TVXC(2C  ,20)  ,T  TVYC ( 20 , 20 ) , TVYC 
5(20,20  )  »CMXB(  4,20  ,CMXT(  4 ,2C  )  ,  WMXB(  4, 2C  )  , kMXT(  4,20)  ,FbMX(  4,20)  ,*lbM 
6XC(4,20)  ,WBMXW(  4,20)  ,CMYb(  3,20)  ,  C  M  Y  T  (  3 , 2  C  )  ,  HM  YB.(  3 , 20  )  ,  WMYT  (  3 , 20  I  »  F 
7BMY(  3,2'-')  , HBMYCI  3,2C),W6MY*l(3,2''),CSX(4,2r),WSX(4,20),FbSX(4,20),k 
8BSX(4,2^) ,CSY<  3,20)  ,WSY( 3,20)  ,FBSY( 3,20) , WBSY( 3,20)  ,CFX(4,20) ,WFX( 
94  *  2°  )  ,  FBFX  (  4 , 20  )  ,WBFX(4,2'1)  »CFY(  3,20)  ,WFY(  3,20)  ,FBFY(  3,  20) 

COMMON  WBFY  (  3,20  )  *DLX(4,20)  ,L)EY  (  3,20)  ,  INSYM,  I  TEFX,  I  TERY,IX,  IY 

REAL  KJ ,KP ,KS ,KX , K Y , K V , K , KXROTC , KXROT W , N YR UTC , K YR OT W , KUX , KD Y , MT 

REAL  KA,KB *KD,KE*KC,KF*KG,KH,KI ,KK,KL ,KM»KN,KU*KU»KR,KT , KM 

DO  612  1=1, NY 

DO  613  J  =  1 » N  S 

DU  614  R= 1 , NS 

SUM=n.o 

DO  615  L  =  1 , NS 

RROD=KY( I , J, L ) *UY ( I , L , M ) 

615  SUM- SUM+PR  OD 
KD Y (  I , J , M ) =  S  U  M 

614  CONTINUE 
613  CONTINUE 
612  CONTINUE 

DC  616  1=1, NY 
DO  617  J=1,NS 
DO  618  L= 1 , NS 

T Y ( I , J,L)  =  YY(  I )*KDY( 1 , J,L  ) 

618  CONTINUE 
617  CONTINUE 

616  CONTINUE 

DO  619  I =1 , NS 
DO  620  J  = 1 , N S 
TT  Y (  I , J  1=0.0 
DO  621  1=1, NY 
ADD  =T  Y  (  L,M  ,  J  ) 

TT Y (  I  , J  )  =TTY (  I  ,  J  )  + ADD 
621  CONTINUE 
620  CONTINUE 

619  CONTINUE 
RETURN 
FND 

SUBROUTINE  VYOTVY 

COMMON  KJ  (20,20)  ,  KP  (2  0,20  )  ,KS(20,20)  ,KX(4,2C,20)  »KY(  3,2  0,20)  *  KV  (  2Ci 
1,20)  ,K(  20,20)  ,NS  ,NX,NY,  I  CX,  ICY,X0(2C  )  ,Y0(  20)  , KXRUTkM  4 , 20 , 20 )  ,  KYRUT 
2  C  (  3,20,20)  ,KYR0TW(3,2C,2r')  ,11(20)  *C  I(2C  )  ,  W  I  (20  ,  WBL  (  20  )  *  ml  i  20  I  ,FBL(2 
30 ) ,wai (20)  ,E, EW, EGC  * EGW , F  X ( 20 ) , F Y ( 20 ) , XX ( 4 ) , YY ( 3 ) , CX ( 20  I , CY ( 20 ) , RX 
4( 4,20) ,DX( 4,20,20) ,KDX( 4, 2 3,20) ,TTX( 20,2^ ), ROT ( 20) ,VX (4,20) ,TVX( 20 
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5) ,RY( 3,20)  ,DY<  3,2  3, 20  ) ,KDY<  3,20,20  ), TTY ( 20,20 ), VY(  ),20)  , TVY( 20) ,MT 
6(20)  ,TVXA(29l  ,TVYA(2C)  , XH ( 20 )  ,YH(  20)  *  FB I ( 20 )  ,KXROTC(4,20*20)  ,  I  NO  I  C 
7X  <  4) , INDICY(  3 ) ,HX(4,20) ,CI X(  4*20)  ,Wl X( 4,20  > ,WHL  X( 4,20) i WXI4.20 ) ,F» 
8LX(4  *20)  , WDI X( 4,20) ,FB1X ( 4,20  > ,HY( 3,20 ) ,C IY(3,20)  ,W I Y( 3,20) , WBLYI 3 

9.20)  ,WY( 3,20)  ,FBLY( 3,20) ,WBIY( 3,20) ,FBtY<  3, 20  > 

COMMON  KA( 20,20 )  ,K» (20,20  ) ,KDI 20, 2C  )  ,KE ( 2C  ,  20) ,KC ( 20,20  I ,KF 1 20, 20 ) 
1 »KG(  20,20 ) ,KH<  20,20 ) ,K I ( 20,20) , KK  I 20,20 ) ,  KL(  20, 20) ,KM( 20, 20) ,KN( 20 
2, 201 , KOI  20,20 > ,KQ( 2C  ,20) ,KR( 20,20 > ,KT( 20,20 ) ,KW< 20 ,20) , TX( 4, 20,20 ) 
3,TY(  3,20,20)  » RT  (  20,20)  »SX(4,20)  ,SY  (  3*20  )  ,R0TXC(4,20)  ,Kl)TXW(4,20)  , 
4R0TYC ( 3,20 ) ,ROTYW<  3 ,20) ,T  TVXCI 20 ,20 ) , T VXC ( 2C , 20 ) , T  T VYC ( 20 , 20 )  , T VYC 
5(  20,20) ,CMXB( 4,20 ) ,CMXT( 4 ,20 ) , WMXB( 4, 20) , WMXT (4 ,20)  ,  FBMXI 4, 20) , WBM 
6XC (4,20) ,WBMXW( 4,20 ) ,CMYB( 3,20 ) ,t MY T( 3,20 ) ,WMYB( 3,20) ,WMYT( 3,20) ,F 
7BMY(  3,20)  ,  WBM YC (  3,20)  ,  WDM  YW(  3,20)  ,CSX(4,20)  , WSX ( 4 , 20  )  ,FBSX(4,20)  ,  W 
8BSX( 4,20 )  ,CSY<  3,20)  ,WSY( 3,20) , F BS Y ( 3,20 ) ,WBSY( 3,20) ,CFX( 4,20) ,WFX( 

94.20) ,FBFX(4,?0),WBFX(4,20),CFYt3,2O),wPY(3,20),FBFY(3,20) 

COMMON  WBFY(3,20) ,Ut X(4,20 ) ,UE Y( 3,20) , I NSYM, I TERX, I TERY, l X, I Y 
REAL  KJ,KP*KS*KX*KY»KV,K,KXROTC*KXROTW,KYPOTC,KYRGTW»KDX*KDY,MT 
REAL  K A  ,  KB  ,  KD , KE , KC »  KPiKG  ,  KH  ,  K  I  , KK , KL , KM , KN , KO , KU , KR , KT , Kw 

DO  622  1=1, NY 
DO  623  J= 1 , NS 
VY( t , J) =0.0 
DO  624  L=! , NS 
PROD=KDY  (  l  ,  J,  L  )  *Rt)T  (  L  ) 

V  Y (  I  »  J  )  =  VY ( I  ,  J)  -PROD 

624  CONTINUE 
623  CONTINUE 
622  CONTINUE 

DO  625  J=1 , NS 

T VY ( J ) =0 . C 

DO  626  1=1 , NY 

TVY I J I =T VY ( J ) +  VY ( I , J) 

626  CONTINUE 

625  CONTINUE 
RETURN 
END 

SUBROUTINE  RTOROT 

COMMON  KJ ( 20,20) ,KP (20,20 ) ,KS<  20,20 ) ,KX (4 ,20, 2C ) ,KY( 3,20,20) ,KV( 20 

1. 20 ) *K(  20, 20), NS, NX, NY, ICX, ICY, X0( 20 ),Y0( 20 ),KXR0TW(4, 20, 20)  ,  KYRU  T 
2C( 3*20,20) ,KYROTW<  3,20 ,2C ) ,H ( 20 ) ,C I ( 20 ) , W I ( 20 ) , WBL ( 20 ) , W ( 20 ) , F BL ( 2 
30) ,WBI I  20) ,E, EW, EGC ,EGw,FX (20) ,FY(20) ,XX(4) ,YY ( 3) , CXI  20) ,CY(20> ,RX 
4(4, 20) , DX( 4,20,20) , KD X ( 4 , 20 , 20 )  ,TTX<  20 ,20) ,ROTI 20) , VX(4 ,20) , TVX(20 
5) ,RY( 3,20) ,UY( 3,20,20) ,KDY( 3,20,2  0  , TTY < 2C ,20) , VY( 3,20)  ,TVY( 20) ,MT 
6(20) ,TVXA(20 ) ,TVYA(20) ,XH(20 I , YH( 20) ,FBI ( 20) , KXROTC ( 4 , 20, 20 ) , INDIC 
7X<4) , I NDICYI 3 ) ,HX (4 ,20) , Cl X( 4,20)  ,WI X<4,20)  ,WBLX<  4,20) , WX (4, 20) ,FB 
8LX ( 4,20 ) , WBIX(4,20) ,FBIX (4,20 ) , HY ( 3,20) ,CIY( 3,20) ,WlY(3,20) , WBL Y ( 3 

9.20)  ,WY( 3,20)  ,FBLY(  3,20)  , W B I Y ( 3,2C  )  ,  F b I Y ( 3,20  ) 

COMMON  KA( 20,20 )  , KB (2  0,20  I ,KD( 2°  ,20 ) , KE  <2C  ,20 ) , KC ( 20,20) ,KF ( 20,20 ) 
1 ,KG( 20 , 20) ,KH( 20 ,20)  , KI ( 20,20)  ,KK (20,20) ,KL( 20, 20 ) ,KM( 2 0,20) ,KN( 20 

2. 20 )  , K0( 20,20 ) ,KU (20 ,20)  , KK( 20, 20) , KT( 20,20 ) ,KW  <20, 20) , TX(4, 20, 20 ) 
3, TY( 3,20  ,2  0) , RT( 20,20)  ,SX( 4,20) ,SY( 3,20 ) ,ROTXC( 4,20) ,ROTXW( 4,20) , 
4R0TYC( 3 ,20 ) ,ROTYW( 3,20) , TT VXC ( 20 , 20 ) , T VXC ( 20 , 20 ) , T T VYC ( 20 ,20 ) , T VYC 
5(20,20)  ,QMX8( 4,20) , CM  XT ( 4,20)  ,wMXB( 4,20 , WMXT (4,20)  ,FBMX( 4,20)  ,WBM 
6XC ( 4 , 20 ) , WBMXW( 4 ,20) , CMY B ( 3 , 20 ) , C MY T  <  3 , 20 ) , WMYB(3, 20) , wMYT( 3,20) ,F 
7BM Y (  3,20 ) , WBM YC ( 3, 20 )  , WBM YW( 3 , 20  I  ,C SX ( 4 , 20 )  ,WSX ( 4 , 20 ) , FBSX( 4,20) , W 
8BSX(4, 20), CSY(3, 20), WSY(3, 20 ),FBSY(3, 20), WBSY (3, 20), CFX(4, 20 ),WFX( 

94.20)  ,FBFX(4,20)  ,  WBE X ( 4 , 2n ) , C F Y ( 3 , 20 ) , WF Y ( 3 , 2 0 ) , F BE Y ( 3 , 20 ) 

COMMON  WBE Y ( 3,20  )  , D E X ( 4 , 20 ) , DE Y ( 3 , 20 )  , I  NS YM , 1 T E RX , I TE R Y , I  X , I Y 
REAL  KJ  »KP  ,KS  »KX»KY,KV ,K , KXROTC , KXR UTW , KY ROTC , KYROT W , KDX , KDY , MT 
REAL  KA  ,KB  ,KD,KE , KC , KF , KG , KH , K I ,KK,KL , KM , KN , KO , KQ , KR , KT , KW 

DO  530  1=1, NS 
DO  531  J= 1 , NS 

RT(  I , J ) =TT  X<  1 »  J ) +TTY (  I  , J ) +KV(  I , J) 

531  CONTINUE 
530  CONTINUE 

CALL  TOKYO  ( N  S , R  T ) 


,  MS..  )  <  )/■:■*'  t"  .  »*  I  *  »  <  °  Ml  H»l*i 


’  ;  ,  r  .  •  I  I  v  f  '  -  ’  ■  '  •  »  •  •  1 ' 


' 


j  . 

' 


' 


A  ',[  ■■,  l 

!t 

n  A ,  «  > 

,  H.o  •  Hi  M  •  ■>,  1,  *  .  *  .  a* 

J 

a,  »  W*<fc%ll*TT*lH  t>? 


d  Of- 


532  CONTINUE 

OH  533  1=1, NS 
KOT( I ) =o.O 
DO  534  J= 1 , NS 
PROD=RT< I , J  )  M  -M  T  (  J I  ) 

ROT (  [  )  =ROT (  I  ) +PRCD 
534  CONTINUE 

533  CONTINUE 
RETURN 
END 

SUBROUTINE  SWAY 

COMMON  KJ( 20, 2C )  ,KP ( 20,20 ) ,KS I  20, 20 I , KXI4.20, 20 ,KY( 3 ,20, 20  *  K  V ( 20 

1 .20)  , K( 20 , 20 )  ,NS  ,NX ,NY , I CX , I CY , Xu ( 20  I , YU( 20  > , KXROT W ( 4 ,20,20 )  ,  K  YRU  T 
2C(  3,20,20) » K  YROTM  3,20,20) »H ( 20 ) ,CI (20) ,WI (20) ,WBL( 20) , W  ( 20 ) ,EBL( 2 
30  )  »  W6I  (20)  *E»Ew,EGC,EGW,FX(2C)  ,FY(20)  »XX(4)  ,YY  (3)  ,  CX  (  20  )  ,CY(20)  ,RX 
4(4,20)  ,DX(4,2r,2ri),KL)X(4,2  0,20)  ,TTX(2C  »2P)*KOT(20)*VX(4,20),TVX(20 
5)  ,RY(  3,20)  ,UY(  3,20,20  ,KDY(  3,20,20  ,TTY(  20,20  )  ,VY(  3,20)  ,TVY(  20)  ,MT 
6  (  20  ,  TVXAl  20  )  ,TVYA(20  )  ,XH{  20  ,  YH(  2C  )  ,  FBI  (  20  )  , KXROTC  (4, 20,20)  ,  INUIC 
7X(  4)  ,  I  NDICY(  3  )  ,HX(4 ,20  )  ,  CIX(4,20  ,WIX(4,2C  )  ,WkJLX(  4,2C  )  ,  rtX  (  4 , 20  )  ,  E  b 
8LX(4,20),WBIX(4,2^),FBIX(4,20),HY(3,20),CIY(3,20),WIY(3,20),HBLY(3 
9,20  ,NY( 3,20)  ,FBLY( 3,20)  ,WBIY<  3,20) ,FBIY( 3,20 ) 

COMMON  KA( 20,20) , KB (20,20) ,KD( 20,20 ) ,KE( 20,20 ,KC( 20,20) ,KF ( 2C, 20 ) 
1,KG( 20,20) ,KH( 20, 20 ,KI ( 20,20) ,KK( 20,20 ) ,KL( 2C, 20) ,KM( 20,20 ) ,KN( 20 

2.20 )  ,KU( 20,20)  ,KQ( 20,20 ) ,KR(2C ,20 ) ,KT( 20,20 ) ,KW (20,20  », TX( 4, 20, 20 ) 
3»TY( 3,20,29) »RT( 20,20)  ,SX(4,20)  ,SY(3,2C)  , ROTXC ( 4 , 20 ) , ROTX Wl 4 , 20 )  , 
4R0TYC( 3,20  )  , ROT Y W ( 3 , 20 ) , T  TVXC<  20,20) ,TVXC (2C , 20)  , TTVYC (  20 ,20 ) , TVYC 
5( 20 , 20 ) ,CMXB( 4,20)  , CM  XT  <  4 ,20 ), WMXB( 4, 2C ), WMXT ( 4 , 20)  , FBMX( 4,20)  , WBM 
6XC(4, 20 ), WBMXW( 4,20) ,CMYB( 3,20) ,CMYT( 3,20 ) ,WMYB( 3,20) ,WMYT( 3,20) ,F 
7BMY  (  3  »  20  )  ,  WBMYC  (  3  ,2 r  )  ,  W'3M  Yri  (  3  ,  Zr'  )  ,  C S  X  (  4 , 2C  I  ,  WSX  (  4 , 20  )  ,  F BSX (  4 , 20  )  ,  to 
8BSX(4,20) ,CSY< 3,20)  ,W3Y( 3 ,20)  ,FBSY{ 3, 20) ,taBSY(3,20> ,CFX(4,20) ,WFX( 

94.20)  ,FBFX(4,2C ), NBFX (4,20 ,CFY( 3,20) ,WFY( 3,20) ,FBFY( 3, 20) 

COMMON  WBFY ( 3,2C ) ,UfcX (4,2C ) ,DLY ( 3,20 ) , INSYM, I TERX, ITERY, IX, 1 Y 
REAL  KJ,KP,KS ,KX,KY,KV,X, KXROT C , KXROT W , KYROTC , K YROT W , KD X , KUY , MT 
REAL  KA.KB ,KD, KE , KC ,KF,KG  ,KH ,K I  , KK , KL , KM, KN,KU, K0,KR , KT,KW 

IF  (  INSYM.EQ.  1  (GO  TO'  114 
DO  ln0  1=1, NX 
DU  101  J  =  1 » N  S 
S  X  (  I , J ) =0  .  o 
DO  102  L=1 , NS 
PROD=DX ( I  ,  J , L ) *ROT ( L ) 

S X (  I  , J )=SX( I , J H-PROD 

102  CONTINUE 
101  CONTINUE 
ICO  CONTINUE 

DO  103  1=1, NY 
DO  104  J  = 1 , N  S 

sy( i , J)  =o.n 

DO  105  L  =  1  ,  NS 

PROD=DY ( I,J,L)*(-RUT(LI) 

SY(  I , J  )  =  S Y (  I , J ) +  PROD 
1C  5  CONTINUF 
1C4  CONTINUF 

103  CONTINUF 
114  CONTINUE 

DO  106  1=1, NX 
DO  107  J=1 , NS 
ROTXC ( I , J ) =0.0 
ROTXW ( I , J ) =0.0 
DC  103  L= 1 , NS 

PR0D1=KXR0TC( I ,J,L)*SX( l , L) 

RGTXC( I , J)=KCTXC ( I , J  H-PR0D1 
PR0D2  =  KXRUTN (  I  , J , L ) *S X ( I  , L ) 

ROTXw ( I , J ) =RCTXW ( I , J ) +PR0D2 
1C 8  CONTINUE 
107  CONTINUF 
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106  CONTINUF 

DO  109  1  =  1  *  N V 

DO  1 10  J  =  1  ,  N  S 
ROTYC ( 1 ,J)=n.O 
ROTYW (  I  ,  J ) =0.0 
DO  111  1=1, NS 

PROD 1=KYR0TC l  l ,  J  *  L  )  *  S  Y  (  I  ,  L ) 

ROTYC.  (  I,J)=ROTYC(  I.Jl+PRODl 
PR0D2=KYRUTW( I , J * L ) * SY (  I  , L ) 

ROTYW (  I  ,  J|=ROTYW ( I , J I +PR0D2 

111  CONTINUF 
110  CONTINUE 
109  CONTINUE 

IF!  INSYM.EQ. I )G0  TO  112 
CALL  WRITE3 

112  CONTINUE 
CALL  CALPAK 

113  CONTINUE 
RETURN 
END 

SUBROUTINE  WRITE3 

COMMON  K  1(20,20)  .  KP  (20,20  )  ,KS(  20,20  )  ,KX  1.4, 20,20)  ,KY  (3,20, 20)  ,KV  (20 

1 .20)  »  K  (  20 ,20 ) ,NS ,NX,NY, I CX, ICY, X0<  20) ,Y0( 20) , KXROT W ( 4 , 20 , 20 )  ,KYKOT 
2C( 3,20,20) ,KYROTW(3,20,20 ) ,H(20 J ,C I ( 20) ,WI (20i ,WBL( 20) , W< 20) ,FBL(2 
30) ,WBI ( 20)  ,E,EW,EGC ,EGW,FX( 201 ,FY(20)  ,XX(4)  ,YY( 3) , CXI  20) ,CY( 20) ,RX 
4(4,20) ,DX( 4,20,20) , KDX( 4, 20, 20) ,TTX( 20,20 ) , ROT ( 20 ) , VX ( 4 , 20 ) , TVX ( 20 
5)  ,RY( 3,20) ,DY( 3,20, 20 ) ,KDY( 3,20  ,20)  ,TTY(  20,20) ,VY( 3,20) ,TVY( 20) ,MT 
4(20) ,TVXA( 20 )  ,TVYA<  20) ,XH( 20) , YH( 20) ,FBI ( 20) , KXROTC ( 4 , 20 , 20 >  , INDIC 
7X14)  »INUICY(3)  ,HX(<t,20),CIX(4,20)  »  W I  X ( 4 , 20 )  ,WBLX(4,20)»WX(4,20)»FB 
8LX(4,20) ,WBIX(4,20) ,FBIX(4,20) ,HY(3,20) ,CIY(3,20) ,W1Y(3,20) ,WBLY( 3 

9.20)  ,  WY( 3,20) ,FBLY( 3,20) , WBIY( 3,20) ,FBIY l 3,20) 

COMMON  KA<  20, 20)  , KB (20,20) ,KD(20 ,20 ) ,KE(20, 20 ) , KC ( 20,20 ) ,KF (20,20 ) 
1 , KG (  20,20) ,KH (20  ,20)  ,KI ( 20,20) , KK( 20,20 ) ,KL (20,20) ,KM( 20,20) ,KN( 20 

2 .20)  , K0(2C ,20) ,KQ( 20,20) , KR120 ,20) ,KT(20,20) ,KW( 20,20) , TX(4, 20, 20 ) 
3, TY( 3,20  ,20)  ,RT( 20,20)  ,SX (4,20 ) ,SY( 3,20 ) , RQTXC ( 4 , 20 ) , ROTXW( 4 , 20 ) , 
4R0TYC ( 3,20  )  , ROT YW( 3,20) , TT VXC ( 20 , 20 ) , T VXC < 20 , 20 ) ,  T  T  VYC  (  20 , 20  ) ,TVYC 
5(20,20)  ,CMXB(4,20) , CM  XT (4, 20) ,WMXB(4,20) , WMXT ( 4 , 20 ) , F8MX ( 4, 20 ) ,WBM 
6XC ( 4 , 20 ) , WBMXW( 4,20 ) , C MY B ( 3 , 20 ) , C MY T ( 3 , 20 ) , WM YB ( 3 , 20 ) , WMYT ( 3 , 20 ) , F 
7BMY( 3,20) , WBMYC( 3,20) , WBMYWt  3,20 ),CSX(4,20),WSX(4,20),FBSX(4,20),W 
BBS X( 4, 20) ,CSY( 3,20) ,WSY( 3,20) ,FBSY( 3,20) ,WBSY( 3,20) , CF X ( 4 , 20 ) , WFX ( 

94.20)  ,FBFX( 4,20) ,WBFX(4,2U) ,CFY( 3,20) ,WFY( 3,20) ,FBFY( 3,20 » 

COMMON  WBFY (3,20),DEX(4,20),0EY(3,20) , I  NS YM , I T ERX , I TERY, I  X, I Y 
REAL  KJ  ,KP,KS ,KX,KY,KV, K, KXROTC, KXROT W , KYROTC , KYROTW , KQX , KDY , MT 
REAL  K A  » KB , KD , KE , KC  »  KF , KG , KH , K I , KK  »  KL  »  KM , KN  »  KO , KQ  » KR , KT  »KW 
WRITE! 6,252) (TVXA( I) ,1=1, NS) 

252  FORMAT (1H1  ,  (4X, • RESULTANT  SHEAR  IN  X-DIRECTION  =',E13.6/)) 

DO  260  1=1, NX 

WRITE (6,261) ( VX ( I, J),J=1,NS) 

261  FORMAT(/1HO, (4X, • INDIVIDUAL  FRAME  X  SHEAR  =',3E13.6/)» 

260  CONTINUE 

WRITE (6, 253) < TVYA( I ) ,1=1 ,NS) 

253  FORMAT  (  1H*1  ,  (4X, 'RESULTANT  SHEAR  IN  Y- DIRECTION  =',E13.6/)) 

DO  262  1=1, NY 

WRITE(6,263)(VY(I,J) ,J=1,NS» 

263  FORMAT( /1H0, (4X, • INDI VIDUAL  FRAME  Y  SHEAR  =',3E13.6/)) 

262  CONTINUE 

WRITE( 6,254) ( XO (  I  )  , I  =  1 , NS  ) 

254  F0RMAT(1H1,(4X, 'CENTRE  OF  ROTATION  IN  X-DIRECTION  =•  ,F11.2/» ) 
WRITE(6»255I (Y0(  1  I  ,1  =  1, NS) 

255  FORMAT ( 1H1 , (4X, ' CENTRE  OF  ROTATION  IN  Y-DIRECTION  =',F11.2/)) 
WRITE!  6,25  6)  ( ROT ( I  I , 1  =  1, NS) 

256  FORMAT( 1H1 , <4X, • ROTATION  IN  RADIANS  =',E13.6/») 

WR I TE ( 6 , 25  7) ICX 

257  FORMAT( 1H1,4X, 'NUMBER  OF  ITERATION  CYCLES  X  = • , 13) 

WRITE ( 6,2bBi ICY 
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?SH  FORMAT  |  IHo  t  »  '  NUMBL  R  Li)  lUHATIUN  CYCLtS  Y  =  ',|)| 

WRITE  (6.27MI  X 

,?70  FORMAT  (  IFI0,4X  ,  "NUMBER  OF  KESLI  CYMh  k  ',li) 

WRITE  (  (St  2  71  I  I  Y 

271  FURMAT(lHO,4X, "NUMBER  OF  RESET  CYCLES  Y  -  *,Is) 

259  COM  I  NUL 
RETURN 
END 

SUBROUTINE  ICXO 

COMMON  KJ  (  2r  *20)  *  KF1 1  2 0 , 20  )  » KS ( 29  »  20  I  ,  K  X  (  4 ,2".  2C  )  ,  K  Y  (  3 , 20,  20  )  ,  K  V  (  2( 
l  ,20  )  ,K(  20 ,2°  I  ,NS  *NX  ,  NY  ,  I  CX,  1  CY  ,  X(J<  20)  ,  YG(  29  )  ,  KXR0TWI4 , 20, 20)  »KYRUl 
2C(  3,20,20 J ,KYRUTw( 3,20,20)  ,Ml 2 M  ,C I ( 20) ,w I ( 20) ,*DL<  20) ,*( 20)  ,FBL ( 2 
3n  )  ,  W  B  I  (20)  ♦  E  »  E  to  »  E  GC  »  E.G  W  *  F  X  (  20  )  *  E  Y  {  2  0  1  ,  X  X  (  4  )  ,  Y  Y  (  3  )  ,  C  X  (  2  0  )  ,CY(2C  )  ,RX 
4(4 ,20)  ,UX(  4,20,20  )  ,  KUX(4,  20,20  )  ,TTX(  20 ,2C  1  ,  ROT  (  20)  ,  VX<  4 ,20)  ,  TVX(  20 
5  )  ,RY(  3,20)  ,I)Y(  3 ,20, 20)  ,KDY(  3,20 ,20)  ,T  TYt  2C  ,2C  )  ,  VY(  3,?0)  ,TVY(  20)  ,MT 
5(20)  ,TVXA(  20  I  ,TVYA(  2C  )  ,XM(  20  ),  YH(  20  ),  FB  I  (  20  )  ,KXROTC  (<♦,  2  j,  20)  ,  I  NOIL 
7X  (  4  )  ,  I  NO  I C  Y  (  3  )  ,  HX  (  4 , 2  0  )  ,CI  X(4,21')  ,WIX(4,?r)  ,WHLx(4,2r)  ,  mX(4,20)  ,FB 
rtL  X  (  4 , 2  0  )  ,  WO  I  X  (  4 , 20  )  ,  F  t»  I  X  (  4  ,20  )  ,  HY  (  3 , 20  )  ,  C  1  Y  (  3 , 20  )  ,  tj  I  y  (  '3 , 20  )  ,  WOL  Y(  3 

9.20)  *  W Y ( 3,2  0)  ,FBLY( 3,20)  , WBIYI 3,20)  , FBI Y< 3,20  ) 

COMMON  KA(  20,2^  ),  KB  I  2  vOt  20  )  ,KD(  2C,2C  )  ,KE  (2C. ,  20  )  ,KC(  2C  ,  20)  ,KF  (  20,20  ) 
1, KG (2 0,20)  ,KH(2O,20) ,K1 (20,20)  ,KK(20,2r  )  ,KL (20,20)  ,KM(2C,20)  ,KN(20 

2. 20 )  ,KU(20,20) ,KQ( 20  ,2r )  , KR(20,20 ) ,KT(2(  , 2C  )  ,KW ( 20,20) , TX(4,20, 20) 
3,  TY(  3 ,  2  o  ,  2  r  )  ,RT(2"',2C)  ,  S  X  ( 4 , 20  )  ,  S  Y  (  3 , 2C  )  ,  ROT  XC  (  4 , 20 )  ,-ROTXW  (  4 , 20  ) , 
4R0TYC ( 3,2E  >  »RUTYto(  3,23  )  ,  TTYXCt  20 ,20  )  ,  TVXU  20,2C  )  ,  TTVYCI  20 ,20)  ,  TVYC 
51  2^,20  )  ,CMXB<  4 ,2C  )  ,CMXT  (  4  ,20  >  ,  WMXB(4, 2°  >  ,  WMXT  {  4,20)  ,  FtiMXl  4, 20)  ,  HbM 
6XC  (4 ,20)  ,  WBMXW(4,20  )  tCMYi)  (  3,20)  ,CMYT  (  3,2C  )  ,  WMYB(  3,20)  ,WMYT  (  3,20  ,  F 
7BMY(  3  ,2  0  ,  WBMYCt  3,2r  )  ,  toBMYW(  3,29)  ,CSX(  4,20)  ,toSX<4,20  )  ,FBSX<  4,20)  ,  w 
8BSX(4,20 ) ,CSY( 3,23) ,W§Y( 3,2n) ,FBSY( 3,20) ,WBSY( 3,20) ,CFX(4,29) , wFX( 

94.20 )  , FBFXt  4, 2C )  , WQF  X(4,20 ) ,CFY<  3,20) ,WFY  <  3,20) ,FbFY( 3, 20) 

COMMON  i*BFY  (  3,21  )  ,0EX(4,2C  )  ,DEY(  3,2r  )  ,  INSYM,  ITERX,  I  TERY  ,  I  X,  I  Y 
REAL  KJ  ,KP  ,KS  ,  KX  ,  KY  ,  KV  ,K  ,  KXROTC  ,  KXROT W , KYROTC , k YROT W ,  KDX ,  KUY  ,  MT 
REAL  KA  ,KB  ,Kf},KE  ,KC  ,KF  ,KG  ,KH  ,  KI  ,KK,  KL  ,  KM.KN,  K(),KiJ,KR  ,  KT  ,KW 

00  ?  33  N— 1 , NS 
1)0  2  39  1  =  1  ,  NX 

K  X (  I  ,  N ) =R  X ( I  , N ) +1.0 
0X(  I  ,  N  ,  N  )  =  0  X  (  I  ,  N  ,  N  )  *  1 . 0 
IF (N.EU.NS)GO  TO  239 
()X<  I  ,  (  N+l  )  ,N)  =0 X  (  I  ,  (N+  1)  ,  N)  -  1 .0 
239  CONTINUE 

CALL  OXCTTX 
? 4 ■'  CONTINUE 

CALL  RTORCT 

242  CONTINUE 
CALL  VXOTVX 

243  CONTINUE 

DO  244  J= 1 , NS 
T  VXC ( N, J) =T  V  X ( J ) 

244  CONTINUE 

DO  245  1=1, NX 

R  X (  I »  N ) =RX ( I , N ) -1.0 

L)X(  1  ,  N  ,  N  )  =  0  X  (  I  ,N,N)  -1 .0 

IF (N.EU.NS )G0  TU  245 

DX( I , ( N+ l ) ,N)=DX( 1 , (N+l) ,N)+1.0 

245  CONTINUE 

DO  246  J=1 ,NS 

TVXC( N,J )=TVXC(N,J ) -TVXAI J ) 

246  CONTINUE 
23B  CONTINUE 

DO  24  7  N  = 1 , N S 
DO  248  J  = 1  *  N S 
TTVXC(J,N)=TVXC(N,J) 

248  CONTINUE 

247  CONTINUE 

CALL  TOKYO  (NS.TTVXC) 

249  CONTINUE 

DO  250  1=1, NS 
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DU  251  J=1*NS 

PRnn=rrvxcd  ,j)m<-fx(j)  >  -  (  t  v  x  a  <  j  > )  i 
xni i ) =  xo ( i ) +Pkou/4.o 
25  1  CUNT  I  NUF 
250  CONTINUF 

DO  252  I = 1 »  NS 

I F ( X  0 (  I  )  • G  T . 2  0 . 0  *  X X ( N X ) ) GU  10  255 

252  CONTINUE 
GO  TO  254 

253  CONTINUE 
I  X  =  I x+l 

DO  254  1=1 , NS 

XO ( I  1 =CX ( I  ) *<  1  . 2  +  (0.2*IX)) 

254  CONTINUE 
RETURN 
END 

SUBROUTINE  ICYO 

COMMON  KJ(23»20>  ,.KP(20*20  )  »KS<  20,2C  >  ,KX(4 ,20,20)  »KY(3,23,23)  ,KV(20 

1 .20)  ,K( 20*20 ) * NS*NX,NY, I CX, ICY  , XOt  2^  » ,Y0( 20  I *KXR0TW<4 ,20,20)  ,KYRUT 
2C<  3,20,20) ,KYRQTW(  3,23,20)  ,H(29 ) ,C  I  <  20  )  ,WI  (20) ,WBL ( 20  ), W( 20) ,F0L( 2 
30 )  ,WB1 ( 2^)  ,E, EW,EGC  ,EGW,FX (20) , FY (20) , XXI 4) , YY ( 3)  ,CX( 20) ,CY( 20) ,RX 
4(4,20) , DX ( 4 , 2  0 ,20  )  ,  KUX ( 4  *  2 r , 2r  )  ,TTX(2°, 20) , HOT ( 20 ) ,VX(4,20),TVX(2r 
5),RY(3,2r'),UY(3,2'',2r'),K0Y(3,2C,2r),TTY(2r,2O),VY(3,2C),TVY(2C),MT 
6(  20)  ,  TVXM  2~  )  ,TVYA(  2n  )  ,XH(  20)  ,  YH(  20  ),  FB  I  (  2'i  )  ,KXKOTC  (4,20,20  I  ,  INDI  C 
7X(  4)  ,  I  NO  1 C  Y  (  3  )  ,HX(4,2"')  ,  C  l  X  (  4 , 2  r  )  ,  W  I  X  (  4 , 2  0  )  ,WdLX(4,20)  ,  WX  (  4 , 20  )  ,EB 
8LX(4,2  0),WBlX(4,20)  ,  F  6  I  X  (  4 , 2  0  )  ,  H  Y  (  3 , 2  0  )  ,CIY(3, 20)  »  W  I  Y  (  3 , 2  0  )  ,WBLY(  3 

9.20)  ,m/Y(  3,2°)  ,  F  B  L  Y  (  3 , 2  0 )  ,WBIY(3,20)  ,  F  B  l  Y  (  3  *  2  0  ) 

COMMON  KA(  20,20)  ,  KB  (2  0,20  )  ,  ISO  (20 ,20  )  ,  KE  (  2  C  ,  20  ) ,  KC  (  2C  ,20)  ,KF  (20,20  ) 
1 , KG (20  *20) , KH (20  ,2C ) ,KI ( 23,20 ) ,KK( 20*20 ) ,KL (20,20) , KM (20,20 ) ,KN( 20 

2.20)  ,K0( 20 ,20) ,KQ (2r ,20 ,KR(2C , 2C ) ,KT (20, 20 ) ,KW(20,20) , TX(4,2C, 20) 
3, TY( 3,20*20) ,RT<  20,20)  ,SX (4,20 ,SYt  3,20) ,RUTXC( 4,20) ,ROTXW( 4,20) , 
4RUTYC( 3,20 ) ,RUTYW( 3,20) , TTVXCI  2?  ,2C ) ,TVXC( 20 ,2C ) ,TTVYC( 20,20 ) ,  TVYC 
5(20,20)  , C  M  X  B ( 4,20)  ,CMXT ( 4 , 2  C  )  ,  w  M  X  b  (  4 , 2  C  ) »  W  M  X  T (4,20)  ,FbMX(4,20)  ,  Wb  M 
6XC( 4,20 ) , WBMXW( 4 ,20 ) , CMYB ( 3,20 ) ,CMYT ( 3, 20 ) , wMYB ( 3 ,20 ) , WMYT( 3,20) ,F 
7  B  M  Y  (  3,2^)  ,  Wt3MYC(  3,2C),W3MYW(3,20),CSX(4,20),WSX(4,20),FBSX(4,20),W 
88SX(4,2'>),CSY(3,?)),WSY(  3,20  )  ,  FBS  Y  (  3 , 2C  )  ,  WB  SY  (  3 , 20  I  ,  CF  X  (  4 , 20  )  ,  WF  X  ( 

94.20)  ,FHFX( 4,20) ,WBFX( 4,20), CFY( 3,20) *WFY( 3,20) ,FBFY( 3,20) 

COMMON  wBFY  (  3, 2  0  )  ,DEX(  4 ,20)  ,  DEY(  3 ,2°> ),  I  NSYM,  I  TERX,  I  TERY,  I  X,  I  Y 
REAL  KJ,KP,KS , KX , K Y , K V , K , KXRO TC , KXROT w , KY RUTC , K YROT w , KDX , KDY , MT 
REAL  KA , KB ,K0, KL , KC ,KF ,KG  ,KH ,KI  ,  Kk  ,  KL , KM, KN , K0 , KU , Kk , KT  ,  KW 

DO  238  N= 1 , NS 
DO  239  1=1, NY 
R Y( I ,N)=kY< I ,N)  +  1 .0 
DYl  I , N , N  )  =  D  Y (  I,N,N)  +  1.0 
IF (N.EU.NS)GO  TO  239 
DY( I , ( N  +  l ) ,N) =DY( I , (N+l) ,N)-1.0 
239  CONTINUF 

CALL  OYOTTY 
24°  CONTINUE 

CALL  RTOROT 

242  CONTINUE 
CALL  VYOTVY 

243  CONTINUE 

"DO  244  J  =  1 , NS 
TVYC(N, J)=TVY(J) 

244  CONTINUE 

DO  245  1=1, NY 

RY ( I ,N) =RY (  1  ,N)-1 .0 

DYl  I , N  »  N ) =  D  Y (  I ,N,N) -1 .0 

IF (N.EU.NS)GO  TO  245 

DY ( I  ,  <  N  +  l )  ,N)  =DY(  I  ,  (  N+l  )  ,  N)  +  1.0 

245  CONTINUE 

DO  246  J= 1 , NS 

TVYCtN, J)=TVYC(N, J)-TVYA(  J) 

246  CONTINUE 
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218  CONTI  NUF 

00  247  N=1 ,NS 
DO  248  J  =  1  ,NS 
T T  V  Y  C ( J,N) - T  V  Y  C ( N,J  I 
248  CONTINUE 
247  CONTINUE 

CALL  TOKYO  (NS,TTVYC) 

240  CONTINUE 

DO  250  1=1, NS 
00  251  J= 1 , NS 

PROD=tTVYCtI  »  J  )  *  ( l-FY(J)  )  —  (TVYAIJI  )  ) 

Y0( I ) =Y0( I ) F PkOD/4. C 

251  CONTINUE 
2^0  CONTINUE 

DO  252  1=1 , NS 

IF(YQ(I  ).GT.2r.C*YY(NY))GC  TO  253 

252  CONTINUF 
GO  TO  254 

.’5  3  CONTINUE 
I  Y  =  I  Y  ♦  1 
DO  254  1=1, NS 

YOm=CYm*(i.2+(0.2*IY)l 
254  CONTINUE 
RETURN 
END 

SUBROUTINE  CALPAK 

COMMON  KJ(  20 ,29)  ,KP ( 2^,20 ) ,KS<  20,20) ,KX(4 ,20,20 ) ,KY( 3,2U, 20)  ,KV( 20 
1 ,20 ,K(  20 ,2^ )  ,NS ,NX ,NY, ICX, ICY , XO ( 2T )  , Y0( 20) , KXROTrt ( 4 ,20,20  I  ,KYKUT 
2C(  3 ,20,20)  ,KYROTW(  3,20  ,20  )  ,H(  20  )  ,C1  I  20  )  ,V.  I  (20)  ,  WBL  <  2C  )  ,  w  <  20  )  ,  F  BL  (  2 
3F  )  ,  U  B I  (  2  0 )  »E*EW»EGC,Eow,FX(20)  ,  F-  Y  (20)  ,XX14)  ,  Y  Y  (  3  )  ,CX(20),CY(20)  ,RX 
4(4,20 , DX(4, 2C,2C ) , KDX( 4,20,2C ) ,TTX(2C,2C ) , ROT (20) , VX(4,20> , TVX(2C 
5) ,RY( 3,20)  ,DY( 3,20 ,20 )  ,KDY( 3,29  ,20) ,TTY<  20,20 ) ,VY<  3,20 ,TVY<  20) ,MT 
6(20  , TVXA( 201 ,TYYA( 20 )  ,XH (20 , YH( 20) , FBI ( 20) ,KXROTC (4,20,20) , INDIC 
7X(4),INDICY(3)»EIX(4, 20  »CIX(4, 20, WIX(4, 20), WBLX(4,2C),wX{4, 20  ),Fb 
3LX<4,20> , WBIX ( 4,20 ) ,FUI X ( 4,20) , HY( 3 ,2C ) ,C IY ( 3,20) ,Wl Y ( 3,20) , WBLY( 3 
9, 2r  )  ,  WY<  3,20)  ,FBLY(  3, 20  ,  WBI  Y  (  3,20  ,FBI  Y(  3, 20  ) 

COMMON  KA( 20,20 , KB (20,20 ) ,KD( 20  ,20 ) ,KE(2C ,20 ) ,KC ( 20,20)  ,KF( 20,20) 
1 , KG ( 20,20 ,KH(2C ,20 ), KI ( 20 ,20 ), KK < 20 , 2C ), KL ( 20, 20 ), KM ( 20 , 20 ), KNI 20 
2 , 2r  ) , KU( 2C ,2C ) *  KQ  (  2'' ,201 ,KR(2C,20) ,KT (20, 2n) ,KW(20»2G) , TX(4,20,20) 
3, TY(  3,20 ,20) ,kT(  20,2  0  ,SX (4,20) ,SY( 3,20 ) , RUT XC( 4,20) , ROTX w( 4 , 20 ) , 
4R0TY  C ( 3, 2C ) ,ROTYW( 3,20) ,TTVXC( 20,20 ) ,TVXC (20,20  ,TTVYC( 20 ,20) ,TVYC 
5 ( 20 ,20) ,CMXB ( 4 ,20 ) ,LMXT ( 4 ,20 ) , WMXB(4, 2C ) , WMXT ( 4,20) ,FBMX( 4,20) , WBM 
6XC  (.4,20  )  ,  NBMXWt  4,2  0  >,CMYB(3,2rl,CMYT(3,20),WMYB(3,2G),WMYT(3,20),F 
7BMY(  3,2'"')  , WBMY  C (  3,2'  )  »  WHM  Y „ (  3,2^  )  ,  C  SX  (  4 , 20  )  , w  S  X ( 4 , 20  )  ,FBSX(4,20)  ,W 
8BS  X  (  4  »  2  ’.’ )  ,CSY(  3,20)  ,WSY(  3 ,2C  )  ,FBSY(  3,20)  ,WBSY(  3,20)  ,  CF  X  (  4 , 20  )  ,  WF  X  ( 
94,20  ),  fiBFX  (  4,20  ),  HBFX  (4,20)  ,CFY(  3,20)  ,WFY  (3,20  )  ,FBFY(  3, 20) 

COMMON  WBFY  (3,2C  )  ,D«EX(4,2C)  ,DEY(  3,20  )  ,  INSYM,  ITERX,  I  TfcRY,  I  X,  I  Y 
REAL  KJ ,KP ,KS , KX , KY , KV ,K , KXRUTC, KXKGTW , KYRUTC , KYR OT w , K D X , KD V , MT 
real  ka,kb  ,KD,KE  ,KC,KF,ieG, KH ,K1  ,KK,KL*KM,kN,KO»KU,KR»KT ,KW 

THIS  SECTION  KEADS  IN  JOINT  RUT  AT  ION, STORY  SrtAY, MEMBER  PROPERTIES 
AND  THEN  CALCULATES  ALL  THE  MEMBER  FORCES 

DO  101  L= 1 , NX 
CALL  CALCKL) 

Ol  CONTINUF 

DO  104  L= 1 , NY 
CALL  CALC2 ( L  ) 

104  CONTINUE 

CALL  WRITE4 

RETURN 

END 
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SUBROUTINE  CALCULI  . 

COMMON  KJ<  20,20)  ,KP  (  20,20  )  ,KS  <  20,20  )  ,KX(  4, 20, 2C  )  ,KY(  3,20,20)  ,K\ll  20 

1 .20)  ,K<  20,20) ,NS  ,NX,NY, I CX, 1CY,XO(20 ) ,Y0( 20) , KXROTW ( 4 , 20, 20  I  ,KYR01 
2C( 3,20,20) ,KYRUTw( 3,20,20 ) ,H ( 20 ) , C I ( 20 ) , W I ( 20 ) , WBL ( 2C I , W ( 20 ) , F BL ( 2 
30I,WBI(20),E,EW,FGC,EGw,FX(20J,FY<20),XX<4),YY(3),CX(20),CY(20),RX 
4(4,20),UX(4,20,20),KDX(4,20,2P),TTX(20,20),RUT(20),VX(4,20),TVX(20 
5)  ,RY(  3, 201  ,L)Y(  3,20,20)  ,KOY<  3,20  ,20)  ,TTY(  20,20)  ,  VY<  3,20)  ,TVY(  201  ,MT 
6(20) ,TVXA(20) ,TVYA(20I  ,XH ( 20 ) , YH ( 20  I , F B I ( 20 ) , KXROTC ( 4 , 20 , 20 1  ,INDIC 
7X  (  4  I  , TNDICY( 3) ,HX(4,20) ,CI X(4,20)  , W  I  X ( 4 , 2  C )  ,WBLX( 4,201 ,WX(4,20) ,FB 
8LX(4,20), WBIX( 4,20) ,FU1X ( 4,20) ,HY( 3,20) ,CIY ( 3,20) ,WIY( 3,20) , WBL Y ( 3 
9 »20 ) ,WY( 3,20) ,FBLY( 3,20) ,WBIY( 3,20) , FBI Y( 3, 20 ) 

COMMON  KA<  20, 20) , KB (20,20) ,KD< 20,20 ) ,KE (20,20 ) ,KC(  20,20) ,KF (20,20) 

1 , KG ( 20,20) ,KH (20,20)  ,KI ( 20,20)  ,KK(20,20)  , KL ( 2  0  *  20 ) , KM ( 20 , 20 )  ,KN( 20 

2 .20)  ,K0( 20, 20  I  ,KU( 20 ,201  , KR (20 , 20 ) ,KT (20,20 ) ,KW ( 20 ,20 ) ,  rx(4,20, 20  I 
3,TY(3,20,20),RT(20,20),SX(4,20),SY(3,20),R0TXC(4,20),ROTXrt(',,20J, 
4R0TYC(3,20),R0TYW(3,20), TTVXCI 20,20 ) , TVXC ( 20, 20 ) , TTVYC ( 20 ,20 ) , TVYC 
6(20,20) ,CMXB(4,20) ,CMXT(4,20) ,WMXB(4,20) *  WM XT (4,20) ,FBMX(4»20) ,  WBM 
6XC(4 ,20 ) , WBMXW(4  ,20 ) , CMYB ( 3 , 20 > , CMYT ( 3 , 20 ) ,WMYB( 3,2C),WMYT(3,20),F 
7BMY( 3,20) , WBMYC (3, 20), WBMYw(3, 20), CSX(4, 20 ),WSX(4,2G),FBSX(4, 20), W 
RBSXI 4,20) ,CSY( 3,20)  ,WSY( 3,20) ,FBSY(  3, 20) ,  W  B  S  Y  I  3,20)  »  C  F  X  (  4 , 2  0  ) »  W  F  X  ( 

94.20 )  , FBFX( 4,20 ) , WBFX (4,20) ,CFY(3,20),WFY(3,20) ,FBFY( 3, 20) 

COMMON  WBFY ( 3,20 ) ,OtX( 4,20) ,DEY( 3,20) , I  NS YM , I TE RX , I TERY, I  X, I Y 
DIMENSION  SWX ( 4 , 20  I 

REAL  KJ,KP,KS  »KX»KY,KV*K,KXROTC*KXROTW*KYRGTC*KYRQTW,KDX,KDY*MT 
REAL  KA  ,KB ,KD,KE , KC , KF  ,KG  ,  KH , K I  ,KK,KL ,  KM , KN  ,  KU , KU , KR , KT , KW 
DO  106  1=1, NS 
H(  I ) =HX ( L,  I ) 

Cl (  I ) =C  I  X ( L  »  I  ) 

W  I  (  I  ) -W I  X ( L , I > 

WB  L  (  I  )  =  Vm  B  L  X  (  L  ,  I  ) 

W(  I  • =  WX ( L , I ) 

FRL(  I  l=FBLX<  L , I ) 

WBK  I  >  =  WB  I  X  (  L  ,  I  ) 

FBI (  I )  =  FB I  X ( L  ,  I  ) 

SWX ( L  ,  I  )  =  S X ( L  ,  I  ) 

F  BMX ( L , I  ) = ( E  *F  B I (  I  l/FBLI  I  ) ) * ( 6 . 0*R0T XC ( L  ,  I  )  I 

FBSX(  L, I )  =  (E*FBI  (I  )/( FBL (  1 )*FBL ( I ) ) )*( 12.0*R0TXC<L, II) 

WBMXC ( L  ,  I  )  =  ( E  *Wb I (  I  ) / WBL ( I  ))*(( 4 . 0*R0TXC ( L ,  I  )  )  ♦  ( 2 . 0*R0TXW ( L , I  ) )♦(  ( 
16.0/WBL(  I ) )*( W(  I  ) *ROTXW( L , I  1/2.0) ) ) 

WBMXW( L,  I  )  =  ( E*WBI  ( I  l/WBL (I)  ) *( ( 2.C*R0TXC( L, I  ) )  +  ( 4.0*R0T XW ( L , I ) )♦(  ( 
16.0/WBL ( I ) >*( W( I ) *R0TXW( L , 1 >/2.0) ) ) 
WBSX(L,I)=1E*W8I(II/(WBL(I)*WBL(I)))*( ( 6.P*R0TXC( L, I ) ) + ( 6 . 0*R0 TXW ( 
1L , I  ))  +  ( (  12 .O/wBL (  I  )  )*( w( I  1 *R0TXW( L, 1 ) /2 .C ) )  ) 

106  CONT I NUE 

DO  100  I  =2, NS 
N  =  I  -  1 

CMXB(L,I)=(E*CI(I)/HU()*(  (  4.0*  RO  T  XC  (  L  ♦  N )  )  M  2 . 0*ROT  XC  (  L  ,  I  )  )  -  (  6 . 0*S 
1WX ( L , I )/H( I) ) ) 

CMXT(L»I  |  =  (E*CI(  I)/H( I  I I  *  ( (2.C*R0TXC( L,N»  )  +  ( 4  .C *R<J TXC ( L  ,  I  I  ) - (6.0*S 
1WX(L,I)/H(  11)1 

CSX(L,I  )  =  <  (E*CHIIX(H(It*H(I  »)))*((  b  .  0*RCIT  XC<  L  ,  N  )  J  ♦<  b.O*RUTXC(  L,  I  ) 
1 )  -  (  1  2 . 0  *SVi  X  (  L  ,  !  )  XH(  Dll 

100  CONTINUE 

C  M  X  B ( I  ,1)=(E*CI(1I/H( I «»*<<2.0*R0TXC<L, 1$ I -< 6 .0*SW X< L , 1  )/H(  1 ) )  ) 
CMXT(L,  1)=(E*CIU1  /HI  111  *U4.C*ROTXC<L,  1)  I  -  (  6  .  C  *SW  XI  L  ,  1  ) /H  (  1  ))) 

C  SX(  L  ,  1  )  =  (  (  E  *CIf  l  1 1 1  )  )  *U<».C*ROTXC(  L,  1  1 1  -  <  12.  Q*SWX(  L  ,  1 )  / 

rn  m  i ) 

DO  101  1=2, NS 
N=  1-1 

WMX  B ( L ,  I  )=  ( EW*Wl U I /M< I  It *«(  4. 3*RGf <W( L, N ) ) ♦<  2.C*«OT  XW(L, I  I  I -(6.0* 
1 SWX( L  ,  I  ) /H<  l  )  )  ) 

WMXT(L,  I)  =  (bW*Wl  U)/H(  1)1  *<C2.0*R0TXW(L,N)  )  +  J  4.0*«0T  XW  (  L  ,  I  )  )  -  (  6.0* 
1SWX<  L, I ) /H( I )  ) ) 

WSX(  L,I  )  =  (  (EW*WI  U)  /(H(I  )*HU  )  )  t )*( (6 .0*R0TXU(L ,N I  )  *  (  6 . 0*ROT  XW  ( L  ,  I 
1 ) ) -(  12. °*SWX( l , I  I /H (I)  )  ) 
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101  CONTINUE: 

WMXB(L,ll  =  (F.W*Wl  I  ll/HEl)  ) * (  (  2 . 0*RUTXW ( L , 1 ) ) - <  6 . 0*SWX ( L , 1 ) /H (  1  )  )  | 
WMXT I L, 1 )  =  (EW*W l  ( 1 ) /H ( 1 ) ) *( ( 4. 0*ROTXW( L , l ) ) -( 6.0*SWX( L ,  l ) /H(  1) ) ) 
WSX ( L  *  1 )  =  (  (  EW*W I  ( 11/ (H ( l ) *H( 1 ) I  )  )  *(  (6.0*RQTXW(L»l l|-(  12.0*SWX(L,1) 

i  /  h  <  i  m 

WBFX(L,NSI=WSX(L,NS> 

WFX( L,NS)=WBSX(L,NS> 

M=NS-1 

00  102  1=1, M 

N=  t  +  1 

WBFX(L, I )  =  WSX I L , I  I -WSX (L,N) 

WFX( L  » ( NS- I ) )=WFX( L , ( NS+l-I )  )+WBSX( L, (NS- l )  ) 

102  CONTINUE 

103  CONTINUE 

FBFX(L,NS)=CSX( L , NS ) +WBF X < L , NS ) 

CFXl L,NS)=FRSX(I  ,NS ) - W BS X ( L , NS ) 

M=NS-i 

DO  104  1=1, M 
N=  I  + 1 

FBFXIL, I )=<CSX( L, I I-CSXl L,N) I+WBFXIL, I ) 

CFXIL, (NS- I) )=CFX( L  * (NS+l-I ) )+FBSX(L, (NS-I) )-WBSX(L, (NS- I ) ) 

104  CONTINUE 

OE  X ( L , 1  >  =  SWX( L ,1 ) 

DO  105  1=2, NS 

DEX ( L ,  I  )=DFX ( L , (  I  - 1  >  )  +  SW  X  (  L  ,  I  ) 

105  CONTINUE 
RETURN 
END 

SUBROUTINE  CALC2(L) 

COMMON  KJ( 20,201  , KP  <  2  0.20 1 ,KS(  2". 20 ) ,  KX(4 ,20, 20 ) ,KY<  3,20, 20) ,KV( 20 

1.20) ,K(29,20),NS,NX,NY,ICX,ICY,XO(20),YU(20I,KXRGTW(4,20,20) ,K YRUT 
2C( 3,20,201  ,K YROTWI 3,20,20)  ,H ( 2 0 ) ,C I ( 20 > , W I  ( 20 ) , WBL ( 20  I , W ( 20 »  , FBL I  2 
30)  ,WBI ( 20)  ,E»EW,EGC»EGw,FX(2C) ,FY(20I , X  X ( 4  I ,YY(3),CX(20I,CY(20I ,RX 
4(4,20)  ,DX( 4,20,20)  ,KDX(4, 20, 20) ,TTX(20,2C ) , ROT (20) , VX ( 4 , 20 ) , T V X ( 20 
5) ,RY( 3,20) ,DY(3,20 ,20) ,KDY(3,20 ,20 ) ,TTY( 20,20) , VY( 3,20) ,TVY( 20) ,MT 
6(20)  ,  T  v  X  A  (  20  )  »  T  V  Y  A  (  20)  «XH  (  2C  )  »YH(  20)  ,FBI  (20)  ,KXROTC(4,20,20)  ,  I  ND I  C 
7X  ( 4 ) , INDICYI3 ) ,HX(4,20) , CIX(4,20)  ,  W I  X ( 4 , 20 )  ,WBLX(4,20) , wX ( 4 , 20 ) ,FB 
8L  X ( 4 , 20 ) ,WBIX( 4,20) ,FB1X( 4,20) ,HY ( 3,20) ,CIY( 3,20) ,WIY(3,20) , WBLY(3 
9 ,20  ,WY( 3,20)  ,FBLY( 3,20) , WBI Y( 3,20 ) ,FBIY<  3,20 ) 

COMMON  KA( 20,20)  , KB (2  0,20 )  ,KD( 20, 20) ,KE (20, 20 ) ,KC<  20,20) ,KF( 20,20 ) 
1,KG( 20,20) ,KH (20,20 ) ,K  I  ( 20, 2C ) ,KK( 20,20 ) ,KL( 20,20 )  , KM ( 20, 201  ,KN( 20 

2.20)  ,KOl 20, 20) ,KQ( 20,20 ) , KR(2C ,20 ) , XT (20,20 ) ,KW( 20,20), TX( 4. 20, 20) 
3,T  Y(  3, 2n,20(  ,RT(  20,2  0  ,SX  (4,20)  ,SY(  3, 2C  )  ,R0TXC(4,20)  ,R0  rxWU,20l  , 
4R0TY  C( 3,20 ) ,ROTYW( 3,20) , TTVXC( 20,20 ) ,TVXC( 20,20 ) , TTVYCI 20 ,20) , TVYC 
5(20,20) ,CMXB( 4,20) ,CMXT( 4,20) ,WMXB< 4, 20) ,WMXT(4,20) ,FBMX(4,20) , WBM 
6XC( 4,20  >, WBMXW (4 ,20) ,CMYB( 3,20 )  ,C MY T( 3,20) ,WMYB( 3,20) ,WMYTI 3,20) ,F 
7BMY( 3,20) , WBMYC( 3 , 20 ) , WBMYW( 3 , 20 ) , CSX ( 4 , 2C ) , WSX ( 4 , 20 ) , FBSXl 4  ,20 ) f  W 
8BSX( 4,20) ,CSY( 3,20) ,WSY( 3,20) , FBSY( 3,20) ,WBSY( 3.20) ,CFX( 4,*0) ,WFX( 

94.20)  , FBFX(4,20) ,WBFX(4,20) ,CFY(3, 20 ,WFY( 3,20) ,FBFY( 3, 20) 

COMMON  WBFY (3,20) ,DEX(4,20) ,DEY( 3,20  > , INSYM, ITERX, ITERY,IX, IY 
DIMENSION  SWY ( 3  »  20 ) 

REAL  K  J  ,  M5  ,  KS  ,  KX  ,  K  Y  ,  K  V  ,  K- ,  KXROTC  ,  KXROT  W  ,  KYROTC  ,  K  YRUT  w  ,  KDX,  KD  Y  ,  MT 
REAL  KA,KB,KD,KE,KC,KF,KG,KH,KI ,KK,KL,KM,KN,KO,KU,KR,KT,KW 
DO  106  1=1, NS 
H(  I  »=HY( L, I ) 

Cl (  I  ) =C  I  Y ( L ,  I  ) 

W  I  (  I  )  =W I Y ( L , I ) 

WB  L (  I  )=WBLY(L, I ) 

W( I ) =WY ( L  ,  I  ) 

FBL (  I  I =FBLY( L , I ) 

WB 1 (  I ) =WB  I Y ( L  ,  1  > 

FBI ( I )=F8IY(L, I ) 

S  WY ( L , I ) =SY ( L  ,  I  ) 

FBMV( L , I )=(E*FBI ( I  )/FBL(  I  ) )*( 6.0*R0  I YC ( L , I  )  ) 

Fft SYIL, 1 )=( t*FBI  ( I  )/(FBL (  I  )*FBL (!))!*(  12.C*R0TYC (L .11) 
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1 r' 6  CUNT  I  NUf 

DO  1 0‘T  Ml,  MY 
DO  110  J - 1  ,  N  S 
K  0  T  Y  C.  (  1  ,JI=0.0 
KO  T Y  W (  1  ,  J  )  =0.0 
DO  111  L=  1 , NS 

PRUD1=KYR0TC  (  I  *  J  »  L ) *SY (  I  ,L) 

RUTYC (  I ,  J  )=RCT YC (  I , J  )  ♦PKOU1 

PRUt)2  =  K  YROTw  (  I »  J  »  L ) *SY ( I  , L ) 

ROTY  W<  I  ,  J  )  =RCTYw  (  l  ,  J  H-PR002 

111  CONTINUE 
1 10  CONTINUE 
109  CONTINUE 

l  F  (  I NSYM . E  U • 1 )  GO  TO  112 
CALL  WRITE3 

112  CONTINUE 
CALL  CALPAK 

113  CONTINUE 
RETURN 
END 

SUBROUTINE  WRITE3 

COMMON  KJ (20,20)  ,KP (2  0,20 ) ,KS( 20, 20 ) ,KX (4 ,20, 20  » ,KY( 3,20,201  ,KV( 20 
l,2r>,K(20,2"),NS,NX,NY,tCX,ICY,XO(2"),YU(20>,KXK0Tw(4,2n,20)  , K YRUT 
2C  (  3 ,2°  ,2")»KYRQTW(3»20,?Q),H(20)»CI(20)*WI(20),WBL(20(,W<20)*FBL(2 
30) »  WB  I (20) »E,EW»EUC,EGW,FX( 2C  ) , FY ( 20 ) , XX ( 4 ) ,YY( 3) »CX(20),CY(20) ,RX 
4(4, 20),  DX  (4, 20, 20,  KDX  (4,2  0, 20)  ,  T  TX  (  20 , 20  )  ,  ROT  (  20  I  ,  VX  (  4 , 20  )  ,  T  V  X  (  2C 
5)  ,RY( 3,20) ♦  0 Y ( 3  *  20  »  20 )  ,KDY(  3*20,2" )  ,TTY( 20,20) ,YY(3,20) ,TVY(20) ,MT 
6(20) , TVXA< 20  )  ,TVYA( 20)  ,  XH(20) , YH( 20) ,EBl ( 20) , K XRO TC ( 4 , 2 0 , 2  0 )  ,  INUIC 
7X(4),INDlCY(3),HX(4,20>,CIX(4,2"),WIX(4,2C),W8LX(4,2C),wX(4,20),Fb 
8L  X ( 4 ,2" ) »WBIX(4»20) , F  B I  X ( 4,20) , HY ( 3 , 20 ) ,C1Y (3,20) ,  W  I Y (3,20) »  WbL  Y ( 3 

9.20)  ,WY<  3,20)  ,FBLY(  3,  2")  »  W 13 1  Y  (  3,20  ,  FBI  Y(  3,20) 

COMMON  KA(  20, 2f  ),  XB  (20,20)  ,KD(2'»  ,20  )  ,KE(  20,20  )  ,KC  (  20,20  )  ,KF  (20,20  ) 
1  ,  KG  (  20,20)  »  KH ( 2"  ,20 )  »  K  /  ( 2 " *  20 )  ,  KK ( 2" *  2"  )  ,KL(20,20)  ,KM(20,20)  ,KN(20 

2.20)  ,K0(2C  ,20  )  ,K(j(2r  ,  20)  ,  KR(20 ,20)  ,KT(2D,20)  ,KW  (  20,20)  ,  TX(4,20, 20  ) 
3, TY( 3,2"  ,20)  ,RT(20,2C)  ,SX (4,20  > ,SY(3,2C) , ROT XC( 4,20) , ROTXw ( 4 , 20 ) , 
4R0T YC (  3,2" ) , ROT YW( 3,20) , TTVXC ( 2C , 20 ) , T VXC ( 20 , 20 ) , TTVYCt  20 ,2C )  ,TVYC 
5 (2", 20) ,CMXb( 4,2")  , CM  XT (4,20) ,WMXB( 4,20  ,  WMXT (4,20)  ,FBMX(4, 20) ,WBM 
6XC (4 ,2" ) , WBMXW( 4 ,20 )  , C MY B ( 3 , 20  I ,CMYT ( 3 , 2C ) , WM YB  < 3 , 20  > , WMYT ( 3 , 20 ) , F 
7bMY< 3,20 > , WBMYC ( 3,2U , WJMYw( 3, 20 ) ,CSX(4,2" ) ,WSX( 4,20) ,FbSX( 4,20 ) , W 
80SX(4,2r)  ,CSY(  3,20)  ,  ln(  S  Y  (  3,2f  )  ,FbSY(  3,2"  )  ,  wbSY  (  3,20)  ,CFX(4,20)  ,  WF  X  l 
94, 2C ) , FBFX( 4,20) , WBF X ( 4 , 20 ) , C F Y ( 3 , 20 )  , WF Y ( 3 , 20 ) ,F8FY<  3,20) 

COMMON  w’BFY  (  3,20  )  ,DEX(4,20  )  ,OEY(  3,2"  )  ,  I  NSYM,  I  TERX,  I  TERY.I  X,  I  Y 
REAL  KJ  ,KP  ,KS  ,KX  ,KY  ,  KV  ,K  ,  K.XRUTC  ,  KXROTW  ,  KYRUTC  ,K YROTW  ,  KOX.KDY  ,MT 
REAL  KA,KB»KD*KE»KC»KF»KG»KH,KI,KK,KL,KM,KN*KQ«KU*KR,KT»Kw 
WRITE (6, 2 52)  ( TV  X A ( I  I  ,1=1, NS) 

252  FORMAT ( 1H1  ,  (4X, • RESULTANT  SHEAR  IN  X-DIRECTION  =',E13.6/)I 

DO  260  1=1, NX 

WRITE(6,261)  ( V  X (  I, J),J  =  1,MS) 

261  FORMAT( /1HC, (4X, • INDIVIDUAL  FRAME  X  SHEAR  =',3E13.6/)> 

26"  CONTINUE 

WRITE l 6, 2 53) (TVYA( I  )  ,1=1 ,NS) 

253  FORMAT  (  lEfl  ,  (4X, 'RESULTANT  SHEAR  IN  Y-D  I  R  EC  T  i  ON  -  '  ,  E  1  3 . 6/ )  ) 

DO  262  I =1 ,NY 

WRITE (6, 263) (VY(  I  , J)  ,  J  =  1 ,NS) 

263  FORMAT!  /1H0,  (4X,  1  IN  01  VI  DUAL  t-RAME  Y  SHEAR  =  «  ,  3E  1 3 . 6/ )  ) 

262  CONTINUF 

WR ITE ( 6,254)  ( X  0 (  I),  1  =  1, NS) 

254  FORMAT(  1H1  ,  (4X,  '  CENTRE  OF  ROTATION  IN  X-DIRECTION  =',FU.2/)I 
WR  ITE(  6 ,256)  ( YEN  I  ),  1  =  1, MS) 

255  FORMAT ( 1H1  ,  (4X, ' CENTRE  CF  RUT  AT  ION  IN  Y-D  I  RFC  T I  ON  =',F11.2/)( 
WRITER  6',  25  6)  (  ROT  (I  )  ,  I  =  1  ,  NS  ) 

256  FORMAT (1 HI , (4X,'RJTAT  ION  IN  RADIANS  =',E13.6/)) 

WR  I  TE ( 6 , 25  7 ) ICX 

257  FORMAT ( 1H1 ,4X ,' NUMBER  OF  ITERATION  CYCLES  X  =',13) 

WR  I  TE (6,260) ICY 
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RLX(4,20) , WBIX( 4,23) ,FBI X( 4  ,20) ,HY( 3,20 ) ,Cl Y ( 3,20  > ,WIY( 3,20) , WBLY( 3 

9.20)  ,WY( 3,20) ,FBLY( 3,20) ,WBIY( 3,20) ,FBIY( 3, 20 ) 

COMMON  KA(  20,20)  ,KB  (2'0,20  )  ,KD(  2C,20  )  ,KE(2C,20  )  ,KC  (  20,20  )  ,KF  I  20,20) 

1 ,KO( 20,20) , KM (20 ,20)  ,KI (20,20) , KK (20,20 ) ,KL ( 20,20) ,KM<  20, 20) ,KN( 20 

2.20)  ,K0<  20,20 ) ,K0( 20 ,20) ,KR(20, 20) ,KT( 20, 20) ,KWI20,20> , TXI 4,20. 20 ) 
3, TY( 3,20,20) ,RT( 20,20 ,SX (4,20 ) ,SY( 3,20 ) ,ROTXC (4, 20) ,ROTXW( 4 ,20) » 
4R0 1 Y  C (  3,20) , ROT  Y  W ( 3,20) , T TVXC( 2C  ,20)  ,  TVXC ( 20 ,20 ) , T  T VYC ( 20 ,20) ,TVYC 
5(20,20) ,CMXB< 4,20) ,CMXT( 4,20) ,wMXB( 4,20) , WM XT (4,20) ,FBMX< 4,20) ,  WBM 
6XC(4 ,20  », WBMXW( 4,20 ) ,CMYB( 3,20) ,CMYT ( 3,20 ) ,WMYB( 3,20) ,WMYT( 3,20) ,F 
7BMYI 3,20) ,  W  B  M  Y  C  (  3,20)  ,  WB  M  Y  W  ( 3,20)  ,CSX(4,2C )  ,  W  S  X  (  4  *  2  0  ) *  E  8  S  X  (  4 ,20) , m 
BBS  X ( 4, 20 )  ,CSY< 3,201 ,WSY( 3,20) ,FBSY( 3,20) ,WBSY( 3,20) ,CFX( 4,20) , wFX< 

94. 20 ) ,FBFX(4, 20), WBFX(4, 20), CFY(3, 20 ),WFY(3, 20), FBFY(3, 20) 

COMMON  WBFY ( 3,20) ,DEX(4,20),DEY( 3,20) , INSYM, ITERX, ITfcRY,lX, IY 
REAL  KJ,KP,KS,K.X,KY,KV,K,KXROTC, KXP OT W , K Y RU T C , K YR OT W , KU X , KOY , MT 
RE  AL  KA  ,KB  ,  KD  ,KE  ,  KC  ,KF  ,  KG  ,KH,KI  ,  KK  ,  KL  ,  KM ,  KN  ,  KO  ,  KQ  ,  KK  ,  K.T  ,  KW 

DO  10O  1=1, NX 
WR I TE ( 6 , 10  l  )  1 

101  FORMAT (  1H1 ,2X ,' FRAME  X  NUMBER •  , T 2 3 , 1  2  ) 

WRITE(6,102) 

102  FORMAT ( IHO// ,4X, 'COLUMN  DATA') 

WRITE (6, 103) 

103  FORMAT (  1H0/,T10, • STORY  NUMBER ' , T 26 , • JU I  NT  R CT A T I  ON ( R AD )  '  , T49 , 

1 • BOTTOM  MOMENT! KI N) '  , T  72 ,  ' TUP  MOMENT ( M N )  '  , T9 3 , ' SHE  AR (KIP)', Till, 
2' AXIAL  LOAD ( K  IP)  '  ) 

DO  104  J  =  1 , NS 

WRITE(6,105) (J,ROTXC( I ,J) ,CMXB ( I , J ) , C MXT ( I , J ) , C SX ( I , J I , CF X ( 1 , J ) ) 

105  F0RMAT(1HC,T14,I2,T28,E13.6,T51,E13.6,T73,E13.6,T92,E13.6,TU2, 
1E13.6) 

104  CONTINUE 

WKl  r  e  < 6,106) 

106  FORMAT! 1HC //,4X, 'WAl L  DATA') 

WRITE(6, 1C3) 

DO  107  J  =  1 , N $ 

WRITE(6,105)(J,R0TXW(I,J), WMXB (I,J),WMXT(I,J),WSX(I,J),WFXtI,J)i 

107  CONTINUE 
WRITE(6,1C8) 

108  FORMAT ( 1H0//.4X, ' REAM  DATA') 

WR I TE ( 6,109) 

109  FORMAT ( 1HC /, T 1 0 ,' FLOOR  NUMBER'  ,T26 ,' FRAME-BEAM ', T43 ,' WALL-BEAM  MOM 
lENT(KIN)* ,T79  ,' SHEAR (KIP)  *,T1C6, 'AXIAL  LOAD(KIP)'  ) 

WR I TE ( 6 , 1 10 ) 

110  F0RMAT(T26,'M0MENT(KIN)' ,T40, 'COLUMN-END* ,T55,'WALL-END',T72,'FRAM 
1E-BEAM'  ,  T8  7 , 'WALL-BEAM', T100, 'FRAME-BEAM' ,T115* 'WALL-BEAM') 

DO  111  J=  1 , NS 

WRITE(6,112)  (J,FBMX( I ,J)  ,WBMXC(  I ,J)  ,WBMXW( I  ,J) ,FBSX( I , J ) ,WBSX( I ,J) 
1 , FBF  X(  I ,J)  , WBF  X (  I , J ) ) 

112  F0RMAT(1HC,T14,I2,T24,E13.6,T38,E13.6,T53,E13.6,T69,E13.6,T84,E13. 
16, T100,E13.6, T1  1  5 , E 1 3. 6) 

111  CONTINUE 
WRITE! 6,113) 

113  FORMAT ( 1HC//,4X, ' FLUOR  NUMBER ' * T2 5 , ' DEFLEC T ION ( IN ) • ) 

DO  114  J= 1 , NS 

WRITE(6,1J5) ( J,OEX( 1 , J ) ) 

115  FORMAT! 1H0,T8,I2,T25,E13.6) 

114  CONT INUE 
ICO  CONT I NUF 

DO  116  1=1, NY 
WRITE(6,117) I 

117  FORMAT (  1H1 ,2X ,' FRAME  Y  NUMBER '  , T2 3 , I  2 ) 

WR  ITE( 6 , 102 ) 

WRITE (6, 103) 

DO  118  J=1,NS 

WRITE(6,105)(J,ROTYC(  !,J)  ,CMYb<  I,J),CMYT(I,J),CSY(I,J),CFY(I,J)) 

118  CONTINUE 
WRITE( 6,106) 

WR  ITE<  6, 103) 
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on  no  j=i,ns 

WRITE  (  6  ,  1 0  5  )  (  J  ♦  R  0  T  Y  W  (  I.,  J  )  ,  WMY  H  (  1  ,  J  )  ,^y  r  (  I  ,  J  )  ,WSY  (  1  ,  J  )  ,  wF  Y  (  I  ,  J  )  ) 

1 1  R  CONTI NUF 

WR I TE ( 6 , 108) 

WR ITE ( 6, IOR) 

WRITF(hflln) 

DO  120  J=1,NS 

WRITEl6,ll2MJ,FliMY(t,J),WBMYCU,J),WBMYW<I,J>,HBSYU,J),WBSY<I,J) 
l.FBFYI  I  ,  J  >  <  WBF  Y (  I  ,  J)  ) 

120  CONTINUE 

WR I TE ( 6  * 1 1 3 ) 

DO  121  J= 1  *  NS 

WR ITE (ft. 1 15) ( J.DEYI I , J ) ) 

1?1  CONTINUE 
116  CONTINUE 
RETURN 
END 
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